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Abstract: Nonperturbative effects in string theory are usually associated to D-branes. In many 
'"^ cases it can be explicitly shown that D-brane instantons control the large-order behavior of string 

■4— > perturbation theory, leading to the well-known (2^)! growth of the genus expansion. This paper 

presents a detailed treatment of nonperturbative solutions in string theory, and their relation 
^ to the large-order behavior of perturbation theory, making use of transseries and resurgent 

analysis. These are powerful techniques addressing general nonperturbative contributions within 
non-linear systems, which are developed at length herein as they apply to string theory. The 
^ cases of topological strings, the Painleve I equation describing 2d quantum gravity, and the 

CN quartic matrix model, are explicitly addressed. These results generalize to minimal strings and 

general matrix models. It is shown that, in order to completely understand string theory at a 
fully nonperturbative level, new sectors are required beyond the standard D-brane sector. 
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1. Introduction and Summary 



String theory may be defined perturbatively, as a topological genus expansion, in terms of two 
couplings, a' and Qs-, 



where F = log Z is the string free energy and Z the partition function. At fixed genus, the free 
energies Fg[ti) are themselves perturbatively expanded in a', where we recall that, for instance, 
in topological string theory, the {tj} moduli are identified with Kahler parameters in the A- 
model and with complex structure parameters in the B-model. As it turns out this a' expansion 
is the milder one, with finite convergence radius. What we address in this paper concerns the 
topological genus expansion where one is faced with the familiar string theoretic large-order 
behavior Fg ~ (2(7)! rendering the topological expansion as an asymptotic expansion [1]. How 
can one go beyond perturbation theory in gs and define nonperturbative string theory in general? 

Perhaps the most popular approach to the nonperturbative definition of string theory is 
within the context of large N duality [2]. In this framework, the partition function of some 
gauge theoretic system defines, nonperturbatively, a dual large- closed-string background, 
described by some geometry which is itself determined by the particular asymptotic limit under 
consideration. Let us focus, for example, on the rather complete picture of [3]. Here, one starts 
off on the gauge theoretic side with a matrix model with some potential, V{z), and, given a 
classical vacuum — that is, a distribution of the matrix eigenvalues across the several distinct 
critical points of the potential — , the 't Hooft large A^ limit [4] yields a closed string background 
described by the topological B-model on a specific non-compact Calabi-Yau geometry. Notice 
that different choices of classical vacua will yield different large A^ geometries and the same gauge 
theoretic system will thus allow for different large A^ asymptotic expansions, represented by these 
distinct semi-classical geometrical backgrounds. The construction of the large A^ dual in [3] is 
essentially achieved by comparing free energies. On the matrix model side, the 1/A^ 't Hooft 
expansion of the free energy starts off by a choice of a semi-classical saddle-point, described 
by a spectral curve (see, e.g., [5] for a review^). Given this spectral geometry there then exists 
a well-defined procedure to compute the large A^ expansion of the free energy [8] which puts 
the results in [3] on solid ground, with an explicit construction of the genus expansion (1.1) of 
the dual closed string geometry (which can of course be checked by explicit calculations strictly 
within the topological B-model closed string theory). Moreover, this procedure later allowed for 
very interesting extensions of the proposal in [3] to more general topological string backgrounds, 
including duals of closed strings on mirrors of toric backgrounds [10, 11]. 

It so happens that the 1/A^ expansion is a perturbative expansion — in fact an asymptotic 
expansion with zero radius of convergence — and there are nonperturbative corrections of the 
type exp {—N) that need to be taken into consideration. These are associated to instantons and 
from a dual closed string point of view, in the large A^ limit, these corrections typically enjoy a 
semi-classical description corresponding to D-brane instanton effects. So we see that, given a 
gauge theoretic system and considering one of its possible large A^ limits, one obtains a closed 
string dual as semi-classical data of the gauge theory 't Hooft limit, both at perturbative and 
nonperturbative levels. This set-up indeed allows us to move beyond the perturbative 't Hooft 

^Let us further point out that this matrix model problem of constructing the 1/A'" expansion, spelled out in 
[6] and which gained an appealing geometrical flavor in [7], has recently been exactly solved purely in terms of 
spectral geometry in [8, 9], and that these results lie at the conceptual basis of our description above. 
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expansion. The question we address in this work is to which point is the nonperturbative de- 
scription complete, and whether D~branes account for the full semi-classical nonperturbative 
data in such a complete description. That is, if the finite N gauge theoretic partition function 
is the correct nonperturbative definition for closed strings in certain backgrounds, one must also 
understand how this finite system encodes, from a dual spacetime point of view, all semi- 
classical nonperturbative contributions. We shall describe how the resurgent analysis framework 
allows for the construction of exact nonperturbative solutions to rather general non-linear prob- 
lems in terms of transseries solutions (first introduced in the string theoretic context in [12]), 
and we shall further fully develop this framework as it applies to string theory. Furthermore, 
we shall find that these transseries solutions, encoding the complete nonperturbative content of 
the large N description, have sectors which cannot be associated to D-branes, at least not in a 
straightforward fashion (as first anticipated in [13, 14]). The resurgent nonperturbative solutions 
have a holographic flavor, in the sense that although one starts from the gauge theoretic (matrix 
model) side, these solutions may be understood in terms of dual large N data. Setting up a 
nonperturbative large N duality framework is of obvious relevance to many diverse issues. For 
example, a particularly interesting question is whether going beyond perturbation theory around 
some classical vacuum of the matrix model will allow, in the holographic dual, to "see" other 
closed string backgrounds (which are naturally included in the finite N gauge theoretic system). 

The present research program started in [10, 15], which proposed to generalize many of the 
nonperturbative results previously obtained within minimal strings to the realm of matrix models 
off-criticality and topological strings (intimately related via [3], as aforementioned). Indeed, the 
double-scaled instantons uncovered in [1], and studied from the matrix model point of view in 
[16, 17], were instrumental for, e.g., the discovery of D-branes in critical string theory sparked in 
[18]^. The approach of [15] used saddle-point techniques to extend results such as [16, 17, 21, 22] 
away from criticality. This is an approach that relies very much on the matrix model spectral 
curve, identifying instantons with B-cycles in the spectral geometry, and which can also be 
extended to the study of multi-instanton corrections — which was done in [23, 14], albeit not very 
explicitly on what concerns general multi-cut saddle-point configurations. These results were 
later extended in [24] to further include instantons associated to A-cycles, which play a relevant 
role in many topological string theories (in the so-called c = 1 class). 

A complimentary approach was introduced in [12], this time around making use of orthogonal 
polynomial techniques [25], where transseries were first introduced to deal with string theoretic 
problems. One of the results in the present paper is to fully generalize these ideas to obtain 
complete nonperturbative solutions to matrix models. In some sense, and as we shall make much 
more precise as we go along in this work, the transseries approach amounts to summing over 
all possible backgrounds, i.e., all possible distributions of matrix eigenvalues across the many 
cuts, which correspond to all possible large A'^ saddle geometries. In particular, multi-instanton 
corrections within multi-cut geometries [23, 14] amount to the exchange of matrix eigenvalues 
along the different cuts, which is effectively interpreted as a change of semi-classical background. 
This naturally led to the construction of a grand-canonical, manifestly background independent, 
partition function in [26] (building upon results in [27, 28]) which was further proved to be 
both holomorphic and modular covariant. Summing over all possible backgrounds or over all 
possible nonperturbative instanton corrections amounts to the same effect. This grand-canonical 
partition function is built by making use of theta-functions, implying, in particular, that there 
will be regions in the gauge theory phase diagram where there are no large expansions {i.e., 

^Of course these instantons also played a decisive role in many nonperturbative questions addressed within 
minimal string theory [5, 19] and were later precisely identified as D-brane configurations in [20]. 
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there is no expansion due to the oscillatory nature of the theta-functions). This rather 
important idea was later explored, from a large N duality point of view, in [29]. Finally, most 
of the transseries results extend beyond the context of matrix models. All they require is the 
existence of a string equation [25], typically a finite difference equation in the context of off- 
critical matrix models, or a differential equation in the context of double-scaled minimal strings, 
which is known to also exist in other examples of topological strings without a very clean matrix 
model relation, e.g., [15, 30]. There may well be larger classes of examples where this is the case. 

A very important role in all this analysis was played by the relation of instantons to the large- 
order behavior of the string perturbation theory [10, 15], i.e., to the fact that these instanton 
effects are testable via their connection to the large-order behavior of the 1/A^ asymptotic ex- 
pansion [31]. Rather impressive agreement was found for many of the calculations in the previous 
references and this will also be a very important point in the present paper: the resurgent frame- 
work we uncover, from an analytical approach, is extensively — and extremely rigorously — tested 
by exploring the connections between the asymptotics of multi-instanton sectors as dictated by 
resurgence. As we shall explain, resurgence demands for a very tight web of relations in between 
all these distinct nonperturbative sectors, which is translated into their large-order behavior. 
These relations may be very thoroughly checked, and to very high precision, making use of 
numerical tests, a fact which will clearly establish the construction we shall purpose. 

It might be fair to say that the first truly unexpected result along this line of research 
appeared in [13], which addressed the large-order asymptotic behavior of mti/ii-instanton sectors, 
rather than just focusing on the usual large-order behavior of perturbation theory. In particular, 
that work addressed the large-order behavior of the 2-instantons sector in the Painleve I system 
(the (2, 3) minimal string) and found that new nonperturbative sectors, besides the usual multi- 
instanton contributions, were required in order to properly describe the full asymptotic behavior 
of this sector. This was done at leading order, in the resurgent framework, and was in fact the 
main motivation beyond the full construction we embrace in our present work: to understand 
the complete set of nonperturbative contributions demanded by resurgence, within the minimal 
string context, and further extend it to general matrix models and topological strings. At this 
stage the reader might complain that we have mentioned the word "resurgence" a lot but have 
been a bit vague about the nature of this framework. This is due to the fact that this formalism, 
a rather general framework introduced in [32] to address general solutions of non-linear systems 
in terms of multi-instanton data, is a bit involved. In here, we wanted to motivate the need for 
more general approaches to nonperturbative issues within large duality from a purely string 
theoretic point of view. We have set up section 2 to introduce this formalism, as well as to show 
how it may be used in string theory, and we recommend the reader to regard this section as an 
enlarged introduction to the ideas that are then explored at length in the rest of the paper. 

This paper is organized as follows. As just mentioned, we begin in section 2 with an introduc- 
tion to resurgence. Asymptotic expansions, with zero radius of convergence, need to be resummed 
if one is to extract any information out of them. There are, of course, many different possible 
ressumation techniques (see, e.g., [33]) but since our models deal with asymptotic series which 
diverge factorially, the natural procedure to use in this case turns out to be the Borel resumma- 
tion framework. This naturally leads to the resurgent framework of Ecalle which we introduce in 
a string theoretic context in this section. We also discuss the relation to the Stokes phenomenon, 
which was also previously discussed in, e.g., [24, 29]. Then, in section 3, we apply some of the 
ideas of resurgence to topological string theory in the Gopakumar-Vafa representation [34, 35]. 
This is, essentially, an extension of the work developed in the context of topological strings on 
the resolved conifold in [24] . Section 4 starts developing the resurgent framework to more general 
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string theoretic systems, in such a way that we can later apply it to minimal strings and matrix 
models. This is where we develop the main structure of our nonperturbative solutions, which will 
later materialize with explicit results in the following sections. In section 5 we discuss one of our 
main examples, the (2, 3) minimal string theory, which describes pure gravity in two dimensions. 
In this section we shall construct the full two-parameters transseries solution to the Painleve I 
equation, generalizing the work of [13]. Do notice that, for the Painleve I perturbative solution, 
leading asymptotic checks have been carried out in, e.g., [36, 37, 38, 15]. A partial transseries 
analysis was done in [39]. As for its multi-instantons solutions, as mentioned above, leading 
asymptotic checks have been carried out in [13]. Our present analysis extends all these partial 
results to a full general solution. Furthermore, by analysis of the resulting resurgent structure 
we show that this solution has complete nonperturbative information about the minimal model. 
More importantly, in this complete set of nonperturbative data, and besides the standard in- 
stanton or D-brane sector, we find new nonperturbative sectors with a "generalized" instanton 
structure. We perform high-precision numerical tests of all nonperturbative sectors, including 
the new "generalized" instanton sectors, which clearly show the need for all these contributions 
in the full exact solution. We also compute many, previously unknown. Stokes constants of the 
Painleve I equation and of the (2, 3) minimal string theory. In section 6 we analyze the full 
fledged quartic matrix model, starting around the one-cut saddle-point geometry. In a similar 
fashion to what we previously did for the Painleve I equation, we construct the transseries solu- 
tion which yields the complete nonperturbative solution to this matrix model. We further show 
how this solution relates back to the Painleve I transseries solution in the double-scaling limit. 
This section includes a discussion of the new nonperturbative sectors of the quartic matrix model, 
alongside with extensive numerical checks which use the resurgent relations to prove the validity 
of these sectors. We also show that the transseries of the quartic model may be set-up in such a 
way that the Stokes constants of this problem are essentially given by the Stokes constants of the 
(2, 3) model. We close in section 7 with a discussion and some ideas for future work. Do notice 
that our analysis generated a rather large amount of data which, for reasons of space, cannot be 
all presented in the body of this paper. Mathematica files with the relevant data are available 
from the authors upon request. We do however present some partial data in a few appendices, 
to indicate how the full set-up was constructed. 

2. Borel Analysis, Resurgence and Asymptotics 

One framework to address nonperturbative completions of rather general non-linear systems 
is the resurgent formalism of Ecalle [32], building upon results of Borel analysis and Stokes 
phenomena, and we shall briefly review it in this section'^. In short, it amounts to a procedure 
which constructs solutions to non-linear problems by addressing all possible multi-instanton 
sectors, i.e., all possible saddle-point configurations in the path integral. Notice that this means 
that one constructs the full solution perturbatively in the string coupling and also perturbatively 
in instanton number — although each instanton contribution is itself nonperturbative. Besides 
allowing for the construction of nonperturbative solutions, the multi-instanton sectors also allow 
for a quantitative understanding of the large-order behavior of the corresponding perturbative 
expansions around a given, fixed multi-instanton sector (the large-order behavior of the zero- 
instanton sector being the simplest case to analyze), a subject with a long tradition in quantum 
mechanics and field theory, e.g., [40, 41, 31]. Some ideas of resurgence have also been partially 

■^The reason for the term "resurgent" — roughly meaning "reappearing" — will also be explained in what follows. 
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addressed recently within the matrix model context, see, e.g., [10, 15, 12, 26, 23, 24, 13, 14]. At 
least in principle, the multi-instanton information could provide for a reconstruction of the exact 
free energy, or partition function, in any region of the coupling-constant complex plane. 

Let us leave the string theoretic framework on the side for a moment, and consider the free 
energy in the zero-instanton sector of any given model, F{z), given as an asymptotic perturbative 
expansion^, 

^(-)-E^- (2-1) 

9=0 

In many examples one may infer that, at large g, the coefficients behave as Fg ~ g\, rendering the 
series divergent. This essentially means that there will be nonperturbative instanton corrections, 
which behave as e~"'^^, with A denoting the instanton action and n the instanton number [31]. 
By further performing a perturbative expansion around the (nonperturbative) contribution at a 
given fixed instanton number, one finds that the full n~instanton contribution is of the form (see, 
e.g., [15, 12, 42, 13] for discussions in the context of matrix models, and topological and minimal 
strings) 

+00 (n) 

(z) ~ z^""^ e""^^ Yl ^ • (2-2) 

9=1 

(n) 

Here /3 is a characteristic exponent and Fg is the ^f-loop contribution around the n-instanton 

(n) 

configuration. In many examples it is also the case that these coefficients behave as i^g ~ (7!, 
further rendering all multi-instanton contributions divergent. As we shall see, it is possible to 
precisely understand the asymptotics, in g, of the multi-instanton contributions Fg^\ in terms 

of the coefficients Fg^\ with n' close to n. This means that the asymptotic expansions are 
resurgent [42], and we shall dwell on this in the following. 

As an approximation to the exact solution these asymptotic, divergent formal power series 
must be truncated and one is consequently faced with the problem that the perturbative expan- 
sion has zero convergence radius. In particular, if we do not know the exact function, F{z), but 
only its asymptotic series expansion, how do we associate a value to the divergent sum? One 
framework to address issues related to asymptotic series is Borel analysis. Introduce the Borel 
transform as the simple map^ 

^ -\t {s) = —^——, (2.3) 

so that the Borel transform of the asymptotic series (2.1) is 

+00 

m{s) = Y.^s3, (2.4) 

9=0 ^' 

which removes the divergent part of the coefficients Fg and renders B[F] (s) with finite convergence 
radius around the origin in C. In general, however, ;B[-F](s) will have singularities and it is crucial 
to locate them in the complex plane. Indeed, if B[F] [s) has no singularities along a given direction 
in the complex s-plane, say args = 6, one may analytically continue this function on e'^M"^ and 



In the following we shall do perturbation theory around z ~ 00, rather than ~ as usual. 
^Notice that the Borel transform is not defined for a = —1, i.e., for a constant term. Thus, in order to Borel 
transform an asymptotic power series with constant term (denoted the residual coefHcient), one first drops this 
constant term and then performs the Borel transform by the rule presented above. 
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thus define the inverse Borel transform — or Borel resummation of F{z) along 9 — by means of a 
Laplace transform with a rotated contour as^ 

SeF{z)= dsB[F]{s)e-''. (2.5) 

Jo 

The function SgF{z) has, by construction, the same asymptotic expansion as F{z) and may 
provide a solution to our original question; it associates a value to the divergent sum (2.1). In 
the following we shall further define the lateral Borel resummations along 9, Sq±F{z), as the 
Borel resummations Sg±eF(z) for e ~ 0"*^. 

However, if the function 0[F](s) has poles or branch cuts along its integration contour above, 
from to e^^oo, things get a bit more subtle: in order to perform the integral (2.5) one needs 
to choose a contour which avoids such singularities. This choice of contour naturally introduces 
an ambiguity (a nonperturbative ambiguity) in the reconstruction of the original function, which 
renders F{z) non-Borel summable. As it turns out, different integration paths produce func- 
tions with the same asymptotic behavior, but differing by exponentially suppressed terms. It 
is precisely when there are such obstructions to Borel resummation along some direction 9 that 
the lateral Borel resummations become relevant: for instance, in the presence of a simple pole 
singularity at a distance A from the origin, along some direction 9 in C, one may define the Borel 
resummation on contours Cg±, either avoiding the singularity via the left (as moving towards 
infinity), and leading to Sq+F{z), or from the right, and leading to Sq-F{z) (see figure 1). One 
finds that these two functions differ by a nonperturbative term [31] 

Sq+F(z)-Sq-F{z)(x I ds- — joce-^^ (2.6) 

J (A) S-A 

Further nonperturbative ambiguities arise as one reconstructs the original function along differ- 
ent directions (with singularities) in the complex s-plane. As such, different integration paths 
produce functions with the same asymptotic behavior, but differing by exponentially suppressed 
terms. To be fully precise about these, we shall need to dwell into resurgence [43, 44, 42]. 

2.1 Alien Calculus and the Stokes Automorphism 

Let us return to our formal power series (2.1). This asymptotic expansion is said to be a simple 
resurgent function if its Borel transform, ;S[-F](s), only has simple poles or logarithmic branch 
cuts as singularities, i.e., near each singular point uj 

m^s) = " + ^{s - u) l^iil^ + $(, - a;), (2.7) 
zvri (s — u) zvri 

where a € C and $ are analytic around the origin. As it turns out, simple resurgent functions, 
under mild extra conditions, allow for the resummation of formal power series along any direction 
in the complex s-plane, thus leading to a family of sectorial analytic functions {SoF{z)}. For 
more details and definitions, we refer the reader to [43, 44, 42]. 

As should be clear — and up to nonperturbative ambiguities — in different sectors one obtains 
different resummations and one needs to fully understand Borel singularities in order to "connect" 
these sectorial solutions together. The next step in order to analyze these Borel singularities 

®If the original asymptotic series one started ofT with had a constant term, dropped in the Borel transform, one 
may now define the Borel resummation as shown, plus the addition of this constant term. 
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Figure 1: The first image shows singularities along some direction 9 in the Borel complex plane, and the 
contours corresponding to the left and right lateral Borel resummations along such direction. The second 
and third images show how to cross this singular direction, or Stokes line, via the Stokes automorphism: 
the left Borel resummation equals the right Borel resummation plus the discontinuity of the singular 
direction (given by the sum over Hankel contours around all singular points). 



in greater detail is to introduce Ecalle's alien calculus [32, 43]. At its basis lies an operator 
which, when acting on resurgent functions, behaves as a derivative; the alien derivative A;^. 
Bypassing the general definition of alien derivative (again, the interested reader may consult, 
e.g., [43, 44, 42]) let us focus on the case of simple resurgent functions: first, if u is not a singular 
point (simple pole or logarithmic cut), Ai^F{z) = 0. Second, at a singular point lo, and slightly 
rewriting^ the Borel transform as 

B[F] (oj + s) = + B\G] (s) + holomorphic, (2.8) 

27ri S ZTTl 

the alien derivative satisfies 

SA^F{z) = a + S,,g^G{z). (2.9) 

Alien derivatives thus encode the whole singular behavior of the Borel transform (they encode 
how much B[F]{s) "jumps" at a singularity) and allow for the aforementioned "connection" of 
sectorial solutions. Indeed, let us consider a singular direction 6, i.e., a direction along which there 
are singularities in the Borel complex plane. In the original complex z~plane such a direction is 
known as a Stokes line (more on this later). Understanding how to connect the distinct sectorial 
solutions on both sides of such direction necessarily entails understanding their "jump" across 
this direction, and this is accomplished via the Stokes automorphism, &g, satisfying [43] 

Sg+ = Se- o = So- o (1 - Disc^-) , (2.10) 

in such a way that the action of &g on resurgent functions immediately translates into the 
required connection of distinct sectorial solutions, across a singular direction 6. In particular, 

Sg+ - Sg- = -Sg- o Disc g-, (2-11) 

such that Discfii precisely encodes the full discontinuity^ of the resurgent function across 6. 
Geometrically, one may think of Disc g- as the sum over all Hankel contours which encircle each 

^Of course in practice it might be hard to find the function G{z) of which ^^(s) is its Borel transform. 

*A function (j> satisfying Ggfj) = "/"i or, equivalently, Discecp = 0, has no Borel singularities along the ^-direction 
and is known as a resurgence constant along this direction. In particular, in this region its Borel transform is 
analytic and (f) is thus given by a convergent power series. 
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singular point in the 0-direction, on the left, and part off to infinity, on the right (see Figure 1). 
The main point now is that, as it turns out [43, 44], one finds 



6, = exp e--«^A^ J , (2.12) 

\M I 

where {ujb} denote aU singular points along the ^-direction. Exphcitly, for singularities along 
the ^-direction in an ordered sequence, one can write [44] 

Se^F{z) = Se-F{z) + J] 1 ^-(.„^+.„^+....„,). (^^^^ • • • A^„,.F(^)) . (2.13) 

r>l;{ni>l} 

One concludes that, given all possible alien derivatives, this result provides the necessary con- 
nection, and thus allows for a full construction of the exact nonperturbative solution alongside 
with its Riemann surface domain. 

It is rather instructive to explicitly write the Stokes automorphism in our multi-instanton 
setting. Consider the positive real axis, where = 0, and where the Borel singularities are located 
at the multi-instanton points nA with n E N*. In this case: 

©0 = exp (2 e-"^^A„^ ) = 1 + e"^^ + e'^^^ (^2A + \^\] +■■■ ■ (2.14) 



\n=l 



This expression will be rather important in what follows. 

Computing alien derivatives straight from their definition is a hard task. Fortunately, as we 
shall see, there are much simpler ways to compute alien derivatives. In fact, it turns out that 
things will greatly simplify by introducing the pointed alien derivative 

A^ = e-"^A^, (2.15) 

as this operator commutes with the usual derivative [44], 

d 



0. (2.16) 

We shall now turn to explicit computations of alien derivatives in different settings. 
2.2 Transseries and the Bridge Equations 

Having understood the central role that alien derivatives play in the construction of nonperturba- 
tive solutions, the question remains: how to compute them in a — preferably simple — systematic 
fashion? The answer arises in the construction of the bridge equations, constructing a "bridge" 
between ordinary and alien calculus. 

Focusing on our familiar multi-instanton setting, with instanton action A (one may allow 
A to be complex, in which case we shall be addressing the aigA direction in the Borel complex 
plane), let us consider a transseries ansatz for our resurgent function, 

+0O 

F(z,a) = 5^cT-FW(z), (2.17) 

n=0 
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where F^'^\z) is the formal asymptotic power series (2.1), and where 



n > 1, 



(2.18) 



are the n-instanton contributions (2.2), as discussed before — the $n(-2) being further formal 
asymptotic power series. In here, a is the nonperturbative ambiguity or transseries parameter, 
selecting, in specified wedges of the complex z-plane, distinct nonperturbative completions to 
our problem. In the resurgence framework, where transseries also go by the name of resurgent 
symbols along some wedge of the complex plane [43], this is the most general solution to a given 
non-linear system. In this work we shall only be concerned with so-called log-free height-one 
transseries: this means that the ansatz above is a formal sum of trans-monomials z"e'^(^\ with 
a G M and, in our case, S{z) a particularly simple convergent series. More general transseries 
may be constructed, with S{z) a transseries itself possibly further involving compositions with 
exponentials or with logarithms, but we refer to [45] for a complete discussion. 

It is also important to realize that the transseries formalism is a rather powerful technology: 
when inserted in the non-linear equation satisfied by F(z) {e.g., a finite-difference equation in 
the case of matrix models, or an ordinary differential equation in the case of minimal strings, 
as we shall see later), it will yield back the non-linear equation for F^^\z) — which is now to be 
solved perturbatively — ; it will yield a linear and homogeneous equation for F^^\z); and it will 
yield linear but inhomogeneous equations for F^"'\z), n > 2. It is thus feasible to solve for all 
members of this hierarchy of equations and fully compute the transseries solution. Indeed, in 

(n) 

all examples of interest to us, it will turn out that all perturbative coefficients Fg appearing 
in the infinite hierarchy of formal asymptotic power series $n(-2) can be computed by means 
of (non-linear) recursions. It will be further the case that the asymptotics of these transseries 
coefficients Fg"'^ will be exactly determined in terms of neighboring coefficients Fg^ ^ , with n' 
close to n, and in terms of a finite number of Stokes constants (defined in the following). 

Finally, notice that we have assumed the transseries ansatz to depend on a single parameter, 
0", assuming that the resurgent function arises as a solution to some problem depending on a 
single "boundary condition" . More complicated problems could lead to more general transseries 
ansdtze, and we shall see some such examples further down the line, but for the moment we just 
consider the simple case where we may power series expand the transseries ansatz in a single 
parameter, i.e., the transseries is an expansion in C[[2;~^, a e""^^]]. For simplicity, we shall further 
assume that the ^n{z) asymptotic series are simple resurgent functions. 

Given the pointed alien derivative Aiy^ = e~^'^^A^^, £ E N*, one may now compute 



The key point is the following: suppose the transseries F(z, a) is an ansatz for the solution to 
some differential equation, in the variable z. Because the pointed alien derivative commutes 
with the usual derivative, it is straightforward to obtain the (linear) differential equation which 




(2.19) 



n=0 




But, clearly, this will be the exact same differential equation as the one 



that 




(2.20) 
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satisfies — simply because also this derivative commutes with the usual derivative^. Assuming for 
simplicity that the difi'erential equation is of first order, it must thus be the case that 

dF 

AeAF{z, a) = Si{a) —{z, a)- (2.21) 

a relation known as Ecalle's bridge equation [44], relating alien derivatives to familiar ones! 
Recalling that alien derivatives encode the singular behavior of the Borel transform, the bridge 
equation tell us that, in some sense, at these singularities we find back the original asymptotic 
power series we started off with — hence its name as a "resurgent" function. Let us explore the 
implications of (2.21). Spelling it out as formal power series, and given that 7^ 0, V^, this 
immediately implies 

5^(0-) = 0, l>l ^ AiAFiz,a) = 0, 1>1. (2.22) 

While one may generically expect that the proportionality factor Si{a) has a formal power series 
expansion as Se{a) = X^^^ -S*^^^ o"'^, homogeneity in a of the bridge equation (2.21) demands 
k = 1 — I. One may quickly realize this by introducing a notion of degree such that 

deg(cT'^e'"^") =n + m. (2.23) 

In this case degF{z,a) = immediately yields deg5'£(cj) = 1 — £, i.e., 

Se{a) = Sea^-^, i<l. (2.24) 

Plugging this back into the power series expansion of the bridge equation one finally obtains a 
clearer expression for the bridge, or resurgence equations^'' 

where we have used conventions in which vanishes if n is less than zero. This expression 
yields all alien derivatives, in terms of a (possibly) infinite sequence of unknowns Si £ C, i £ 
{1, —1, —2, • • • }, the so-called analytic invariants of the differential equation we started off with. 
Knowledge of these analytic invariants allows for a full nonperturbative reconstruction of the 
original function F{z), the problem we first set out to solve. However, generically, the analytic 
invariants are transcendental functions of the initial data (say, the differential equation one 
started off with) and quite hard to compute. 

For completeness, it is interesting to notice that the above resurgence equations (2.25) may 
be translated back to the structure of the Borel transform, at least near each singularity iA, by 
making use of the definition of alien derivative (2.9) for a simple resurgent function. Indeed, with 

+00 (n) +00 (n) 

^'^(-)-E^ mn]{s) = Y.^^,^^^yS^^-^^\ (2.26) 



9=1 



'in full generality this is slightly more subtle: indeed, it will be often the case that the differential equation one 
is considering will depend on some other ("initial data") functions. In this case, for the above reasoning to hold, 
these functions must either be entire functions, or their Borel transforms cannot have singularities at the points 
£A (of course these functions will also have no dependence on a whatsoever; the transseries expression (2.17) is 
simply an ansatz for the solution, introducing a new parameter). 

^"One also obtains a clearer explanation for the name "resurgent": via the bridge equations the alien derivatives, 
encoding the singular behavior of the Borel transform, are given in terms of the original asymptotic power series 
one started off with (multiplied by suitable Stokes' constants). 
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it simply follows 



B[^n]{s + eA) = Se{n + i)B['^n+i]{s)^^, e<l. (2.27) 



Going back to the connection formulae (2.13) or to Stokes' automorphism (2.14) makes clear 
how important (2.25) is: it is telling us that the somewhat initial multi-instanton data is enough 
for a full reconstruction of the nonperturbative solution. Consider first the positive real axis, 
where ^ = 0, and where the Stokes automorphism is 

©0 = exp c-^^^Am^ = 1 + e-^^ Aa + e'^^^ (^A^a + ^ A^) + • • • . (2.28) 

Given the transseries ansatz, the action of ©q -^(-2, cr) is entirely encoded by the action 
of ©0 on the several $„(z), and this can now be completely determined by the use of the 
bridge equations (2.25). But because these vanish whenever £ > 1, the Stokes automorphism 
immediately simplifies to 

©0 = exp (c-^^Aa) = 1 + e-^^ A^ + ^e'^^^ A\ + ^e'^^^ A^ + • • • , (2.29) 

where 

AT 

A^^n = iSif ■l[{n + i)-'^n+N. (2.30) 

i=l 

One may now simply compute 

©o^'^ = E f"" ^ (2-31) 



The interesting fact about the bridge equations (2.25) is that they contain much more information 
than just that concerning the positive real axis. Indeed, consider instead the negative real axis, 
where 6 = tt, and where the Stokes automorphism becomes 

6. = exp e^^^A„M^ = 1 + e^^ A_a + e^^^ (^A_2A + ^ A^ ^) + • • • . (2.32) 

The action of on F{z, a) is again entirely encoded by the action of &^ on the several ^^(z), 
and is determined by the use of the bridge equations (2.25). All one needs are formulae for 
multiple alien derivatives, which follow straightforwardly as^^ 

AT N N / i \ 

n A_ = n uh-E^A- '^n-Ef-, ■ (2.33) 

i=l 1=1 i=l \ j=l j 

Notice that the ordering of the alien derivatives in the left-hand side of the expression above is 
rather fundamental, as alien derivatives computed at different singular points do not commute. 
For example, it is simple to check that [A_„^,A__m^] oc (n — m). Also, because the alien 



'^^Of course this expression holds as long as n — X]i=i 7^ 0- soon as this term vanishes, so does the multiple 
alien derivative, and consequently so will all subsequent ones. 
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derivatives vanish as soon as one considers A_„^$„ = 0, this apparent series actually truncates 
to a finite sum, at each stage. One may simply compute 

n-i ^ ^ \ ^ ^ ( \ I 

6.<^n = <l>n + E«'^^E^ E n^-vIIh-E^^lh"-^ (2.34) 

i=\ k=i ' £i,..A>i ti=i i=i \ 

n-l £ / k . 

= cl.„ + ^e^^^ n("-7^)^-d7. |^n-^ (2.35) 

1=1 k=l ' 0=7o<7i<---<7fc=^ \i=l 

In the last line, the sum over all possible partitions ii > 1 was replaced by a sum over their 
consecutive sums 7^ = and we defined d'yj = 7^ — 'yj-i {i.e., the partitions). A few 

examples of the Stokes automorphism at = vr are given below: 

Gn^o = $0, (2.36) 

= ^1, (2.37) 

6^«>2 = $2 + 5-1 e^^«>i, (2.38) 

e^^3 = $3 + 2S-1 e^" $2 + (5-2 + 52 1) e^^' $1. (2.39) 

Finally, making use of the Stokes automorphism (2.12), one may directly apply the bridge 
equation (2.21) in order to find, e.g., 

Se+F{z,a) =Se-exp(^A^'^ F{z,a) =Se~F (^z,a{l + u;S^a-'^)^'^ . (2.40) 

In particular, when io = 1, argo; = 0, this is 

S+F{z,a)=S-F{z,a + Si), (2.41) 

in such a way that Si acts as a Stokes constant for the transseries expression. For this reason, 
we shall generally refer to the analytic invariants Se as "Stokes constants". Of course this exact 
same expression could be obtained by applying the Stokes automorphism at = 0, (2.31), to the 
transseries (2.17) (trying the same at 9 = tt, via (2.34), would be much more complicated). In 
the original complex z-plane this Borel-plane singular-direction corresponds to a Stokes line and 
what the expression above describes is precisely the Stokes phenomena of classical asymptotics — 
here fully and naturally incorporated in the resurgence analysis. At a Stokes line, subleading 
exponentials start contributing to the asymptotics and this is accomplished in here by the "jump" 
of a, the coefficient associated to the transseries formal sum over (multi-instanton) solutions. 
In other words, the "connection" expression (2.13) yields a relation between the coefficient(s) in 
the transseries solution, in different parts of its domain, or, on different sides of the Stokes line. 

2.3 Stokes Constants and Asymptotics 

One may wonder why the long detour into resurgence and alien calculus. As it turns out, 
understanding the full asymptotic behavior of all multi-instanton sectors — which is to say, fully 
understanding the nonperturbative structure of the problem at hand — demands for this complete 
formalism. Let us first recall the standard large-order dispersion relation that follows from 
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Cauchy's theorem [31]: if a function F{z) has a branch-cut along some direction, 9, in the 
complex plane, and is analytic elsewhere, it follows 



F(.) = 2, r'-~d„5i!£iZW / (2.42) 



2vri Jfj w — z /(oo) 27ri w — z 

In certain situations [40, 41] it is possible to show by scaling arguments that the integral around 
infinity does not contribute. In such cases Cauchy's theorem provides a remarkable connection 
between perturbative and nonperturbative expansions. Let us first consider our familiar pertur- 
bative expansion (2.1), within the transseries set up (2.17), where F^^\z) = ^o{z). In this case, 
the bridge equations (2.25) tell us, via the Stokes automorphisms (2.31) and (2.34), that F^^\z) 
has the following discontinuities: 

Disc a>o = - 5^ 5f e-'^' (2.43) 

Disc^^>o = 0, (2.44) 

i.e., fW(z) has a single branch cut along the Stokes direction corresponding to the positive real 
axis in the Borel complex plane. From the perturbative expansion (2.1) and using (2.42) above, 
it immediately follows 

Tpjo) r (g - fc/3) ^ r (g - fc/3 - + 1) (fc) h-1 .2 45) 

^ 2ni (kAf-'^^ ^ Tig- fc/3) ^'^ ^'^^ ' ^'■^^> 

where we have used the asymptotic expansions for the multi-instanton contributions (2.2). It is 
instructive to explicitly write down the first terms in this double-series, 

^(0) ^ Si_ r(g-/3) / (1) A (1) \ 

^3 ~ 2^i + 5 _ /3 - 1 2 + J + 

Sf Tjg- 2/3) / (2) 2A (2) \ 

53 r (g - 3/3) / (3) 3 A (3) \ 

+ 2^7^Ij^r^ ^i^W^l"^^ +---J+---. (2.46) 

This is the multi-instanton generalization of a well-known result, also from previous work within 
the matrix model and topological string theory contexts, e.g., [15, 12]. In particular, it relates the 
coefficients of the perturbative expansion around the zero-instanton sector with a sum over the 
coefficients of the perturbative expansions around all multi-instanton sectors, in an asymptotic 
expansion which holds for large g (and positive real part of the instanton action). In particular, 
the computation of the one-loop one-instanton partition function determines the leading order 
of the asymptotic expansion for the perturbative coefficients of the zero-instanton partition 
function, up to the Stokes factor Si. Higher loop contributions then yield the successive ^ 
corrections. Furthermore, multi-instanton contributions with action nA will yield corrections to 
the asymptotics of the Fg^^ coefficients which are exponentially suppressed as n~^. 

The novelty here arises due to the use of alien calculus, which allows for a straightforward 
incorporation of all multi-instanton sectors in the asymptotic formulae, as well as a generalization 



-14- 



of this procedure to all multi-instanton sectors! Indeed, in terms of asymptotics of instanton 
series, we shall find that the bridge equations (2.25) essentially tell us that, given a fixed instanton 
sector, its leading asymptotics are determined by both the next and the previous instanton 
contributions — at least in examples where a transseries ansatz depending on a single parameter 
is enough (we shall later see examples where things get more complicated). In particular, at 
the level of Borel transforms, the singularities closest to the origin, of i3[<I>„](s), are located at 
s = zizA (if n = 0, 1 there is a single closest-to-the-origin singularity located at s = A), and these 
singularities will necessarily control the large-order behavior of the multi-instanton sectors. 

Let us now address the full n-instanton sector. Consider our perturbative expansion (2.2) 
within the transseries set-up (2.17), where F^'^\z) = e~""^^$„(z). In this case, the bridge 
equations (2.25) tell us, via the Stokes automorphisms (2.31) and (2.34), that F^^\z) has the 
following discontinuities: 

Disco^n = - E f''^^) 5f e-^^^$„+,, (2.47) 

n-l e / k \ 

Disc ^<l>n = -Yl e'^' E ^ E n - ^^■) ^-d7. ^n-i, (2.48) 

e=i k=i ' o=7o<7i<---<7fe=^ \i=i / 

i.e., F^"\z) has branch cuts along the Stokes directions corresponding to both positive and 
negative real axes in the Borel complex plane. The contribution from the discontinuity at 9 = tt 
can also be rewritten as a sum over Young diagrams 7^ G r(A;, i) : < 71 < • • • < 7^ = £ of 
length ^(r) = k, and with maximum number of boxes for each part 7i (also called the length of 
the transposed Young diagram) being £(r'^) = £. This sum is only completely well-defined if we 
also set 5*0,70 = 0, in which case one finally obtains 

n-l e / k \ 

Disc^$„, = -Ee'''^EM E n(^-^i)'5-d,, U„-£. (2.49) 

£=1 k=i ■ 7,er(A:/) \j=l / 

For example, consider once again the case n = 3 in the notation above, 

Disc.ci>3 = _e^^ E ((3-7i)5-d,,)^2-e2^^E^ E ( 11 " 7.) S-d,, ) ^i- 
7i6r(i,i) fc=i ■ 7,er(fc,2) \i=i / 

(2.50) 

Expanding the sums, there will be only one Young diagram corresponding to r(l, 1), □, for which 
71 = 1. For r(l, 2) one can only find m (71 = 2), while for r(2, 2) there are two possible Young 
diagrams: P (71 = li 72 = 2) and g (71 = 72 = 2; but this will not contribute because 50 = 0). 
The expected result arising from (2.39) then simply follows. 

From the perturbative expansion (2.2) and the dispersion relation (2.42), which now needs 
to account for both branch cuts, it finally follows 

^'V" f^^^\ ^ r (g - kf3) ^ T{g-k/3-h) in+k) ,...h 



k=\ 

n— 1 I , k 



fc=l I m=l 7ier{m,fc) \j=l 
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This expression relates the coefficients of the perturbative expansion around the n-instanton 
sector with sums over the coefficients of the perturbative expansions around all other multi- 
instanton sectors, in an asymptotic expansion which holds for large g. All Stokes factors are now 
needed for the general asymptotic problem, and this analysis has essentially boiled down the 
asymptotic problem to a problem of precisely computing these Stokes factors. These numbers 
are transcendental invariants of the problem one is addressing and generically hard to compute — 
although, as we shall see, a matrix model computation partially solves this issue. Notice that 
for instanton numbers n = 0, 1 the combinatorial factor associated to the Stokes factors Se at 
negative i vanishes. As such, the contributions arising from the second and third lines in the 
expression above can only be seen at instanton number n = 2 and above. Finally, let us note that 
the explicit treatment of the leading contribution to this type of asymptotics was ffist presented, 
to the best of our knowledge, in [42]. 



3. Topological Strings in the Gopakumar— Vafa Representation 

The first concrete example we shall explore deals with topological string theory, where the free 
energy admits an integral Gopakumar- Vafa (GV) representation, see [34, 35]. Consider the free 
energy of the A-model, on a Calabi-Yau (CY) threefold X, with complexified Kahler parameters 
{ti}. As a string theory it satisfies the standard topological genus expansion (1.1) where, at 
genus g, for large values of the Kahler parameters (the large-radius phase), one finds [46] 

+ 00 

Fg{ti)=J2Ng,d{X)e-^-\ (3.1) 

di = l 

where the sum is over Kahler classes^^ and where the coefficients Ng^^i'^) are the Gromov-Witten 
(GW) invariants of X, counting world-sheet instantons, i.e., the number of curves of genus g 
and degree d in X. As we mentioned before, this a' expansion is the milder one with finite 
convergence radius tc, where the conifold singularity is reached, which may be estimated from 
the asymptotic behavior of GW invariants at large degree (here 7 is a critical exponent; see, e.g., 
[30]) [46] 

Ng,d ~ 6-^*% d ^ +00. (3.2) 

What we shall be interested in next is instead the asymptotic genus expansion. In this case, 
and as throughly investigated for the resolved conifold in [24] , the GV integral representation for 
the free energy may be interpreted as a Borel resummation formula, immediately yielding, as we 
will see in the following, the "leading" part of the topological string resurgent data. 

3.1 Topological String Free Energy and Borel Resummation 

Let us thus consider the GV integral representation for the all-genus topological string free energy 
on a CY threefold X, including nonperturbative M-theory corrections via the type IIA -H- M/S^ 
duality [35] (see also [48, 34, 49] and [24] for a discussion in the present context which further 



^^d = (di, . . . , di,^f^x)) denotes the expansion of the two-homology class d on a basis of H2 {X, Z); see, e.g., [47]. 



highlights the Schwinger-like nature of this result), 



^^(5.)-EE^^''^WE / -(2sm^)''" exp(-^(d.t + im) . (3.3) 

In here, the integers ri^''\x) are the GV invariants of the threefold Af, depending both on the 
Kahler class di and on a spin label r, and the combination Z = d-t + \m represents the central 
charge of certain four-dimensional BPS states [35]. To be completely precise, notice that in 
order to obtain the full topological string free energy one still has to add to this expression the 
(alternating) constant map contribution [50, 51] 

^-°-45(2,-2)(2,-2)!^('^^' (^-^^ 

where x{'^) = 2 {h}'^ — h?'^^ is the Euler characteristic of X. This term can also be written as 
a Borel-like resummation, where the result is^^ [34, 24] 

3=0 m6Z-^° \ 9s J 

Apart from the overall multiplicative factor of the Euler characteristic, the universal constant 
map contribution has already been fully addressed in [24] and we shall thus leave it aside for the 
moment. Let us focus on the GV contribution (3.3) instead. By rewriting the sum in m S Z as 
a sum over delta-functions it is simple to obtain the GV formula for the topological string free 
energy as [35] 

+00 +00 +00 

r=0 di=\ n=l 

This makes it quite clear how the input data for a given CY threefold is simply its set of GV integer 
invariants (and its Euler number if one is also to write down the constant map contribution). 
Expressed as a topological genus expansion one finds, at genus (7, (see, e.5., [50, 47, 52] for partial 
expressions) 



^.('.) = E \ 4^-' + B-l)"-'-^ 4*' [ "3-a» (e— ) , (3.7) 




(n) 

where the coefficients am are obtained from the generating function 



^n(2;) - -j-j^ ^^^^ 1^ ^^^^ - E "L+i^™" (3-8) 



m=0 

by power series expansion, and where Lip(x) is the polylogarithm of order defined as 

+00 



Lip(-) = E^- (3-9) 



n=l 



^■^This may also be obtained directly from the GV representation by simply setting di = and r = in (3.3) and 
properly identifying the "degree zero" and "spin zero" GV invariant with the Euler number of the CY threefold. 
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Two things to notice are the foUowing: at fixed genus g, only GV invariants with h < g 



contribute to the free energy [35]; in particular the "highest" GV invariant at genus g has h = g 
and appears with coefficient one in (3.7) as (^2g_2y, — 1; ^9- Furthermore, ag is only non- 
vanishing when g = 1, implying that n^*^'"* only contributes to the genus one free energy. 

What we want to understand in here is how or when the GV representation (3.3) may 
be understood as a nonperturbative completion of the free energy genus expansion (3.7), in 
the sense of resurgent analysis. Furthermore, one would like to understand how to relate this 
nonperturbative completion to the large-order behavior of the genus expansion (3.7) via the use 
of Stokes' automorphism. Following the approach in [24], we shall interpret the GV integral 
representation for the free energy (3.3) as a Borel resummation formula (2.5), for SgFxigs)^ in 
such a way that, after a simple change of variables, one obtains 

+00 +00 ^ / ^ 2r-2 

g[^.viM = EE"^^-'E^(^^"' 4.(,.:+i„) ) • 

with the GV representation (3.3) now amounting to the statement that 

SeFx{gs) = / dse[F;t](s) e"^. (3.11) 
Jo 

This rewriting, of course, required changing the integration with the (in general) infinite sums 
over GV invariants, a procedure which is only valid if there is uniform convergence of the partial 
sums in (3.10). As we have seen before, the sum in m is the milder one. Furthermore, at fixed 
degree, the sum in r will truncate, i.e., given a fixed two-homology class {dj}, there is r* such 
that ri^''' = for ah r > [53]. The real issue concerning uniform convergence of the GV Borel 
transform thus arises when we fix genus and sum over degree. In this case one finds that the 
asymptotic behavior of, for example, the genus zero GV invariants at large degree is [54] 



(d) ^ exp(27rt2(l) • d) 
'° ^ (P{\ogd) 



2 , d^+oo, (3.12) 



where 27rt2(l) is a critical exponent (for instance, in the example of local this would be 
27rt2(l) — 2.90759... [54]). This is an exponential growth and, as such, in strict validity, the 
results that follow only hold for threefolds with a finite number of GV invariants, i.e., without 
compact four-cycles. This is also in line with the general expectations briefly discussed in [24]. 

In this context, the only singularities of the GV Borel transform (3.10), with s 7^ 0, appear 
when r = as the zeroes of the sine (located at ujn = (2vr)^ n (d • t + im), n G Z*). In this case 
one will only find pole singularities and the Borel transform (3.10) may be written as 

B[F;,]{u:n + ^) = ^ (^^^2 - + holomorphic, (3.13) 

near each singular point The (multiple) instanton action, a;„ = nA, is further obtained as 

Am{ti) = {2Tif{d-t + im). (3.14) 

In the following we shall make use of this information in order to explore, from a resurgent point 
of view, when does the Borel interpretation of the GV integral representation (3.10) provide for 
a nonperturbative completion of the topological string free energy. 
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3.2 Simple Resurgence in Topological String Theory 

The first step in understanding the resurgence of topological strings is to compute alien deriva- 
tives. At first, this could seem non-trivial as the GV Borel transform (3.13) is not quite a simple 
resurgent function due to the second order pole. However, explicitly evaluating the difference of 
lateral Borel resummations as in (2.6), one notices that the contribution from this second order 
pole is simple to include in the alien derivative, which now becomes, for n £ Z* , 

(d^) / (27r)^ jd-t + im) Qs 



^nAFx = n'^'^ + ^ A„, (3.15) 

2-Kgs \ n n^J 

with an added unusual dependence on the coupling constant. In spite of this, the right-hand side 
above is in fact a resurgent constant, in such a way that all multiple alien derivatives vanish. In 
this case, it is trivial to compute Stokes' automorphism, (2.12). Denoting^^ hy 9 = argA, this is 

&j,Fx = Fx^Y.Y.Y.^^-^^'^\ ^ )^ (3.16) 

n=ldi=lmeZ ^ 9s J 

leading to the discontinuity 

^^'^^ - 5^ E E '4*' E m^^^ . I) e--^. ,3.17) 

71=1 di=\ mel ^ ' 

Finally, making use of the dispersion relation (2.42), where one further assumes that the contribu- 
tion around infinity may be neglected, one may now compute all coefficients in the perturbative 
asymptotic expansion of Fx^ which has the usual genus expansion form (1.1). Focusing on the 
discontinuity naturally induced by the GV integral representation (3.3), namely args = 0, and 
following a calculation very similar to the one in [24] for the case of the resolved conifold, it 
follows 



T-OO +00 I „ I 

Fxi^9.) - Y.^r' E 4'^^ Li3-,, (e-— ) . (3.18) 

Some comments are in order concerning this result. The first obvious one is that this does 
not fully match against the GV result (3.7), as it only captures the leading, dominant Bernoulli 
growth of the free energy. While this is certainly the correct expectation for an asymptotic 
formula in the case of a finite number of GV invariants, one may also ask if it is possible to do 
any better. Of course, if one is to start with the GV Borel transform (3.10), its singular part 
(3.13) will not include any GV invariant n^r^^ with r 7^ and, as such, will never be able to yield 
the subleading contributions in (3.7) unless these should arise from the singularity at infinity in 
the Cauchy dispersion relation (2.42). While this is a possibility, it is also a notoriously difficult 
case to handle — the singularity at infinity is an essential singularity, leading us far from the 
realm of simple resurgent functions — further departing from the conventional set-up of resurgent 
asymptotics. At the end of the day this "loss" of GV invariants ri'r'^ with r 7^ arises from 
the exchange of integration and infinite sums in (3.3) to obtain (3.10) and all it says is that 
another procedure will be required in order to look beyond the Bernoulli growth in (3.7), i.e., 
to study the full nonperturbative information of topological string theory. In other words, while 



'Notice that in the original integration variable of (3.3), therein denoted s, this would correspond to S = 0. 
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the GV integral representation is extremely useful in order to solve topological string theory at 
the perturbative level, (3.7), one needs extra work if one wants, in general, to obtain a closed 
form expression for the topological string Borel transform — possibly in terms of GV invariants. 

At this point it might be useful to make a bridge to the case of matrix models with polynomial 
potentials (a subject we shall study in detail later in this paper). For these, the Gaussian 
component of the polynomial potential will induce a contribution to the free energy which also 
leads to Bernoulli growth [24], rather similar to the one above arising from genus zero GV 
invariants. From a spectral curve point of view, both these contributions are associated to A- 
cycle instantons [24]. This type of instantons is always very simple to handle. As described 
above, the alien derivative is essentially trivial (it equals a resurgent constant) and at the end of 
the day the asymptotics is somewhat universal — and certainly much simpler than the discussion 
in the previous section. All multi-instanton sectors have no non-trivial large-order behavior 
(their alien derivatives vanish, their series truncate and their structure is thus rather different 
from the one in (2.2)) and the perturbative sector is essentially dominated by Bernoulli growth. 
For matrix models with polynomial potentials the truly non-trivial resurgent structure will then 
be associated to higher monomials in the potential which will induce different contributions to 
the free energy, this time around associated to B-cycle instantons [15]. More realistic examples 
of this non-trivial resurgent structure associated to B-cycle instantons will be discussed next, as 
we move to the realm of minimal strings and matrix models in the following sections. For the 
moment, let us just notice that, in general, we still expect topological strings to display full non- 
trivial resurgence: if one wants to see beyond the Bernoulli growth in (3.7) one will certainly need 
to find a proper Borel transform, leading to non-trivial alien derivatives and asymptotic growth 
of all multi-instanton sectors. Thus, in general, there will be both A and B-cycle instantons 
in topological string models, both contributing to the full instanton action, and controlling (in 
turns, depending on the absolute value of their corresponding actions) the large-order behavior 
of perturbation theory at different values of the 't Hooft moduli, as recently discussed in [55]. 

4. The Resurgence of Two— Parameters Transseries 

In order to address broader string theoretic contexts, in particular those involving minimal string 
theory or matrix models, as we shall study later in this work, we now need to generalize the for- 
malism introduced in section 2 in order to include transseries depending on multiple parameters. 
Let us start off with some words on the general transseries set-up (see, e.g., [56] for a recent 
review, or, e.g., [57, 58] for more technical accounts). 

A rank-n system of non-linear ordinary differential equations. 



may always be written, via a suitable change of variables, in the so-called prepared form [56]: 




(4.1) 





Denoting by {oi} 



the eigenvalues of the linearized system. 




- i,j=l---n 



(4.3) 



-20- 



then, in the expression above, A = diag (ai, . . . , a„) and B = diag . . . , /3„) are diagonal 
matrices and one further insures that = O z^^ii) . It is also convenient to 

choose variables such that ai > 0. Most cases addressed in the literature deal with the non- 
resonant case, where the eigenvalues {ajjj^^^...^ are Z-linearly independent, in many cases with 
argaj 7^ argOj. This will not be the case in the present work, as the string theoretic systems we 
address resonate. In the above set-up, a formal transseries solution to our system of differential 
equations (4.1) is given by [56] 

ti(z,cr) = wW(z)+ cr"z""-^e~"-"^it(")(z), (4.4) 

neN"\{0} 

where cr = (ai, . . . , a"„) are the transseries parameters, and where both the perturbative con- 
tribution, u^^\z), as well as instanton and multi-instantou"^^ contributions, u^'^\z), are formal 
asymptotic power series of the form 

+00 [n) 

9=0 

The fact that the systems we shall address in the following resonate now translates to 

3ny^n' \ n ■ a = n' ■ (X. (4.6) 

Furthermore, one often deals with proper transseries, where only exponentially suppressed contri- 
butions appear: the eigenvalues cx are such that, for some chosen direction in the complex z-plane, 
all contributions along this direction with ai ^ are exponentially suppressed; Me (n ■ a z) > 0. 
Again, as first pointed out in [13], if one wishes to fully address the instanton series in a string the- 
oretic context one will also have to allow for less studied non-proper transseries. We thus see that 
resurgence in string theory is more intricate than usual, with resonant non-proper transseries. 

As we have reviewed in section 2, asymptotic series need to be Borel resummed in order 
to extract sensible information from them. Naturally, this will also be a required step in the 
construction of a transseries solution to the non-linear differential equation (4.1), and it follows 
that [57, 58] 

Sg±u{z,a±) =Se±u^^\z)+ Yl o-± ^"""^ e""'"" cSe±it(")(z), (4.7) 

neN"\{0} 

is a good solution to our problem along a proper direction (at least for sufficiently large \z\). 
Many of the concepts introduced in section 2 now have a straightforward generalization, for 
instance a simple extension of Stokes' automorphism (2.13) where this time around one may 
write 

So+u{z,a) =Sg-u{z,a + S) (4.8) 

for the crossing of a Stokes line, with S the associated Stokes constants. 

We shall now construct the resurgent formalism for the specific case of two-parameters 
transseries, which will turn out to be the required framework to address the instanton series in 
2d quantum gravity (as first uncovered in [13] for the case of the Painleve I equation) as well as 
the instanton series in the quartic matrix model, as we shall discuss in this work. 

^''Linear systems have no multi-instanton sectors. 
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4.1 The Bridge Equations Revisited 

We have seen in section 2 how the bridge equations allow for a simple evaluation of alien deriva- 
tives (up to the determination of the Stokes invariants), (2.25), and how this result then allows for 
an exact evaluation of the Stokes automorphism along a singular direction in the Borel complex 
plane, (2.31) and (2.34). We have further seen in section 2 how the discontinuities associated to 
these singular directions end up determining the full multi-instanton asymptotics (2.51) and, in 
essence, solve the nonperturbative problem via the use of transseries solutions. 

In general one requires multi-parameter transseries in order to set up full nonperturbative 
solutions which completely encode the multi-instanton asymptotics. For the main examples we 
shall study in this work, the quartic matrix model and its double-scaling limit, the Painleve I 
equation, it turns out that a two-parameters transseries is required, as we shall see later and as 
discussed in [13]. We shall now derive the bridge equations in this situation. 

In particular, we consider the special case of two-parameters transseries where the prepared 
form eigenvalues are {it^}, with A the instanton action^^. The transseries ansatz is now simply 

+ 00 +00 

F(z,ai,a2) = Y.Y1 (4.9) 

n=0 m=0 

where the perturbative asymptotic series is 

+00 (0|0) 
3=0 

and where the generalized multi-instanton contributions take the form^'^ 

+00 rp{n\m) 

^(n|m)(_2) ~ e-n{+A)z g-m(-A)2 ^ _9 _ ^-in-m)Az $^^|^^^(^). (4.11) 

The characteristic exponent is often taken to be of the form Pnm = n/Si + m/32, but we shall also 
allow for more general combinations. Everything else is a a straightforward generalization of the 
standard result (2.2) and a simple application of our discussion at the beginning of this section. 

Because ^(n,ni)jt{n',m') \n — m = n' — m' this transseries describes a resonant system and it is 
not too hard to see that one can make the "instanton number" explicit by slightly reorganizing 
the previous transseries representation, obtaining 

+00 +00 +00 +00 

F{z,a,,a2) = ^<e-"^^ ^ ('Xiaa)™ $(^+n|.„)(^) + E ^2 e"^^ ^ (^1^2)" f Hm+n)(2). 

n=0 m=0 71=1 m=0 

(4.12) 

This also introduces a natural notion of degree, 

deg (cj^'a^e'^^^) =n-m + k, (4.13) 
such that the transseries F{z, cri, CJ2) has degree zero. 



'^^This will be the relevant case for both the Painleve I equation and the quartic matrix model. 

^^In here we are simplifying things a bit: as we shall later discuss in the Painleve I framework, ^(n\r„){z) is not 
always a plain formal power series in 2 but may sometimes also include logarithmic powers, of the form log'° z 
multiplied by formal power series in z. For clarity of discussion, we shall proceed under this simpler assumption. 
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Let us now consider the pointed alien derivative A^^ = e~^'^^A£^, £ G Z*, which, as we 
discussed earher, commutes with the usual derivative. The reasoning of section 2 used in deriving 
the bridge equation also holds now, albeit in the two-parameters case the space of solutions to 
the differential, or finite difference, string equation becomes two-dimensional [13] (we shall see 
this very explicitly in the examples that follow). It must then be the case that 

dF ~ dF 

AiAF{z,ai,a2) = Si{ai,a2) ^— (z, fii, 0-2) + Si{ai,a2) ^— (z, di, 0-2); (4.14) 

Oai (7(72 

the bridge equation in the two-parameters setting. Let us explore its implications. First of all, 
it is quite simple to notice that this immediately determines the degrees of the proportionality 
factors as 

degS£{ai,a2) = I - i and deg ^^(ai, 0-2) = -1 - ^. (4.15) 
Because these should be expressed as formal power series expansions, this further implies 

+00 

Se{a,,a2)= J2 si'^'-'''K',+'-'al (4.16) 

fc=max(0,-l+£) 

and 

+00 

Siia,,a2)= si'^'^'^'^a'l4+'+'. (4.17) 

fc=max(0,-l-£) 

Clearly, there are now a whole lot more Stokes constants than before. For simplicity of notation, 
and noticing that the Stokes constants depend only on the two parameters k and i, we redefine 
them as 

^ik+i-e,k) ^ gik+i^i) ^^^^ gikMi+i) ^ sik+i+i)_ 

Plugging these expressions back into the power series expansion of the bridge equation (4.14) 
one obtains, after a rather long but straightforward calculation, 

min(m,n+f— 1) 

^eA^{n\m) = X] {n-k + i) S'f "^"^^^ ^l^n-k+e\m-k) + 

fc=max(0/-l) 

min{m—£,n) 

+ ^ (7n-A;-^)5f+'+^)<^(„_fc|^_fc_,), (4.19) 

fe=max(-€-l,0) 

valid for all £ 7^ 0. Looking at the i > 1 case {i < —1 is completely analogous), one finds 

^eA^{n\m) = X] (n-k + l) Sf'^ $(„_fc+l|m-fc-£+l) + 

A:=0 
mm{m—£,n) 

+ ^ (m-k-l) 5f ^^n-k\m-k-i) , (4.20) 
fc=0 

which can be directly compared with equivalent expressions from [13]. In these expressions we 
have used conventions in which vanishes if either n or m are less than zero. As compared 

to the one-parameter case, (2.25), the increase in complexity is evident. Analyzing the bridge 
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equations in the form (4.19), it is not difficult to notice that the cases AM^(n|m) and ^-£A^{m.\n) 
with i > are intimately related. In fact, one can go from one to the other by performing the 
simple changes (5^, S^) -f-)- (5*"^, 5^^) and ^*(a|6) ^ where a, b can be any combination of 

indices. The same relation can be seen to extend to the full Stokes automorphisms — changing 
between the ©o^(ri|m) ^-nd S^<I>(m|n) cases — which we shall further discuss in the following. In 
any case, the main focus of our concern deals with the instanton series, ^"(,110)) where these 
formulae become 



^M'^(n|0) 



0, £>1, 

5|°)(n + l)cl.(„+i|0), ^ i = l, (4.21) 



This result clearly illustrates that in the present situation, unlike the one-parameter case, un- 
derstanding the asymptotics of the physical instanton series necessarily requires the use of the 
generalized multi-instanton contributions, due to the appearance of the term in which, 
upon multiple alien derivation, will make materialize the full generalized instanton sector. 

As we have further seen in section 2, the bridge equations may also be translated back to 
the structure of the Borel transform, at least near each singularity in the Borel complex plane. 
In the present case we have to consider, for f3nm = n/^i + m(32, 

+00 p{n|0) +00 p{n\0) 

9=1 3=1 

as well as^^ 



-00 



p(n|l) +00 (n|l) 

Then, from (4.21) above and the definition of alien derivative, it simply follows, e.g., 

B [$(„|o)] {s + M) = (Sf-'^ {n + i)B [^(„+f|o)] (s) + 5^ B [<f („+,+i|i)] (s)) £ < -1. 

(4.24) 

The next step is to use the ahen derivatives in order to fully construct Stokes' automorphism, 
allowing for a full reconstruction of the nonperturbative solution. Consider first the positive real 
axis, where ^ = 0, and where the Stokes automorphism is 

©0 = exp e-^^^AM^ = 1 + e-^^ + e'^^^ (a2A + ^ A^) + • • • • (4.25) 

Just like in the one-parameter case of section 2, given the transseries ansatz, the action of Sq on 
F(z,(7i,(T2) is entirely encoded by the action of ©q the several <l*(n|m)(^)) and this can now 
be completely determined by the use of the bridge equations. When focusing on the physical 
instanton series, and again akin to the one-parameter case of section 2, the bridge equations 
(4.21) vanish whenever i > 1, and when 1 = 1 both one-parameter (2.25) and two-parameters 
(4.21) cases are entirely analogous. Thus, the Stokes automorphism immediately simplifies to 

©0 = exp (e-^^A^) = 1 + e"^^ Aa + ^e'^^^ A\ + ^e^^A. + . . . , (4.26) 



'in the Painleve I case there will also be logarithmic contributions to which we ignore for the moment. 
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where 



N ^ 

^>{n\0) = • n + ^) • '^(n+^IO)- (4-27) 



One may now simply compute 

fc=0 ^ ^ 

a completely straightforward generalization of the one-parameter case (2.31). 

The novelties arise as we turn to the Borel negative real axis, where 6 = it, and where the 
Stokes automorphism becomes 

(+00 \ +00 £ ^ 

5^ e^^^ A_M = 1 + E E ^ E ^-^.^ • • • ^-^.^ = (4.29) 
e=i / e=i k=i ■£i,...,4>i 



= 1 + e^^ + e^^^ (^A_2A + ^) + • • • . (4.30) 

Things are now much more complicated than in the simple one-parameter transseries case, as 
the different terms in &^ will mix contributions arising from all <I>(„|m)- From the expression 
above for the Stokes automorphism it becomes obvious that, in order to find the final expression 
for 6^$(„|^), one first needs to focus on determining A^£^a • • • ^~eiA^{n\o)j with ij > 1. For 
k = 1,2, this calculation is pretty straightforward. Using (4.21) one can write 

A-£M^(n|o) = {n- ii) S^X''^ '^in-i,\o) + S% f („-^i+i|i), (4.31) 



A ...A 



i=l 



+ I (n - h) + I n + 1 - > : 4 1 5i\f + S^3. S% «^ (n+i-E^i ^.|i) + 



+2e\^it^(n+2-EL.^.|2)- (4.32) 

In order to go further and generalize these cases to an arbitrary product of alien derivatives, one 
first needs to determine A^f, , , 4$/' , v-fc « i \- After some effort one can find that 

m+1 / / fc+1 \ ~ \ 

X (4.33) 

where we have set S_f^ = 0, for any £j > 1, in order to simplify the final result. The general case 

for the ordered product of k alien derivatives of the form n^=i ^-ik+i-i^ = ■A-4A " " " A„£^^, 
acting on $(„|o), is then given by 



k k k I s 

n^-4+i-Aio) = En] E 

1=1 m=0 s=l 1^3=0 



\ i=i / \ i=i 1=1 / 
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xe(s-Xg, ) }5{y^qi,k-mj<i>^^_^^_^,^^,^^^y (4.34) 

In this expression 6{n,m) = 6nm is the usual Kronecker-delta, the function G(x) is the usual 
Heaviside step-function 

Q(x) = <( " - (4.35) 




and once again we set S_y = 0. A proof of this result can be found in appendix D. In the same 
manner as we have done earlier in the one-parameter case for the discontinuity at = vr (2.49), 
this result can also be rewritten using a sum over Young diagrams. To do so, let us first define 
6s = Ei=i (Is + 1, such that 0<(^i<(52<---<5fc = A;-m + landO<(5s<s + 1^^. As 
explained in section 2, the set of 6s, with s = 1, . . . ,k, form a Young diagram T{k, k — m + 1) 
of lengths £{T) = k and ^(r^) = k — m + 1, with the extra constraint that each component 
6s S r(^, k — m + 1) has a maximum number of s + 1 boxes. As such, one may finally rewrite 
the above result as 



(s + 1 - 6s) Sif =^ + 



k k ^ ( 

i=l m=0 5ser(fc,fc-m.+l) s=l ^ 



n 




For this expression to hold, one still needs to set 6q = \ and s'^q^ = 5q*^ = 5^^^ = (notice that 
some of these conditions will only be needed in the following) . 

Due to the complexity of this expression, let us pause for an example. Let us choose the 
case of /c = 2, which we have also described in (4.32) above, and see what the sum over Young 
diagrams (4.36) above yields. One finds: 



X 



A_,,^A_,,A$(„|o) = E E (2 - <5i) ((2 - 5i) + (n - ^1 + 2 - 5i) S^X"'"' 

X e (3 - 6^) ((3 - 6^) + (n - £i - £2 + 3 - 52) S^^'^^^) ^ (^.^^^j^.^^ • (4.37) 

The sum over Young diagrams in this expression is over 6s S r(2, 3 — m), with m = 0, 1, 2. For 
m = 0, one sums over all diagrams 6s G r(2,3) and there are three possible diagrams: 0^, 
EEP and EE0- But because 6\ < 2 and 61 = 3, only two will remain: [P^ (where 6\ = 1) and 

{61 = 2). For m = 1, one sums over diagrams 6s G r(2,2) and there are now two possible 
diagrams with 62 = '2: EP {^i = 1) ^tnd ^ {61 = 2). Finally, for m = 2, one sums over diagrams 
6s G r(2, 1), which corresponds to the single diagram: {61 = 62 = !)• Plugging these results 
back into the expression above, one easily finds (4.32) as expected. 

There is now enough information in order to completely determine the Stokes automorphism, 
at 6 = TT, of the instanton series <I*(„|o)- Going back to its definition (4.29) and making use of our 



'^^The reason for adding the one in the present definition of Ss is to make aU Ss strictly positive, and thus 
naturally labeled by some Young diagram. 
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formulae for multiple alien derivatives (4.36) it follows 



+00 i ^ 

e^'J>(n|o) = ^HO + E^'^^'E^ E A_4.^...A„,,A^(„,|o)= (4.38) 
t=i k=i ■£i,...,4->i 

= ^No) + E«"^E^ E E E (4-39) 

£=1 A;=l ^i,...,4->l "1=0 (Sser(fc,fc-m+l) 



n 



i=l 



e (s + 1 - <5,) ^ • $ 



(n—l+m\m) • 



Interestingly enough, if we further define 7^ = 'Ylii=i^i^ then the sum over the £j can also be 
rewritten as a sum of Young diagrams 7^ E r(A;, t) : < 71 < • • • < 7^ = as long as we set 
S^^ = Sq*^ = 0. In this case, one finally obtains the simpler expression 

+00 £ ^ 

e.^(nio) = ^(„io)+EEV E E E (4-40) 

£=1 fc=l ^^^Y{k,t) m=0 5^6r(fc,fc-m+l) 



n { [(. + 1 - 5.) + (n - 7. + s + 1 - <5.) e (. + 1 - 5,) } • <!>(„ 



-£+m\m) • 



Some comments are now in order. Comparing the Stokes automorphism of the instanton series 
at the 9 = TT discontinuity, for both the one-parameter (2.49)^*^ and the two-parameter cases 
(above), one can see the that the increased degree of complexity of the latter is translated in the 
fact that there is now a sum over two independent sets of Young diagrams (instead of summing 
over just one set of diagrams as in the one-parameter case). It is thus natural to infer that 
for a general ^-parameter transseries ansatz such sums would be substituted by sums over i 
independent sets of Young diagrams. It is also not too difficult to see that one can recover the 
one-parameter result (2.49) starting from (4.40) above, by simply setting 5^ = s + 1 for all 6s- 
This corresponds to choosing the Young diagrams 6s G T(k,k + 1) (with m = and consequently 



< n) where each row has one more box than the previous one, e.g. 



for k 



4.2 Stokes Constants and Asymptotics Revisited 

The main outcome of the above calculations are expressions for the discontinuities of the full, 
physical, multi-instanton series, encoded in the Stokes automorphism of $(n|o)) both 9 = 0,tt, 
directions. As we have seen earlier, in section 2, these discontinuities lie at the basis of un- 
derstanding the full asymptotic behavior of all multi-instanton sectors and we shall next use 
these new Stokes' discontinuities in order to generalize our results on asymptotics, from the 
one-parameter to the two-parameters case. Recall that by making use of Cauchy's theorem a 
given function F(z) with a branch-cut along some direction 9 in the complex plane (and ana- 
lytic elsewhere) can actually be fully described precisely by its discontinuity along that direction 
(2.42), at least as long as its behavior at infinity does not contribute. In the present case of inter- 
est, the multi-instanton free energies F^"'^^\z), which are the coefficients of the two-parameters 



^°Recall that 6„ = 1 - DisC( 
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transseries ansatz (4.9), have asymptotic expansions given by (4.10) for the perturbative se- 
ries and by (4.11) for the generalized multi-instanton contributions. Their discontinuities are 
essentiaUy given by the Stokes automorphisms of ^(^n\o)i^) previously calculated. 

Let us first look at the perturbative expansion (4.10). The discontinuities of F^^^'^\z) arise 
directly from the bridge equations (4.21), via the Stokes automorphisms (4.28) and (4.40), 



+ 00 



Disco$(o|o) = -E e-'=^^$(,|o), (4.41) 

k=l 

Disc^$(o|o) = -E (^-O e'^'^(oifc)- (4.42) 



k=l 

Note that F(o|o)(z) will now have two branch cuts in the Borel complex plane (instead of only one 
as in the one-parameter transseries case), along both positive and negative real axes. By using 
(4.10), (4.11) and (2.42) it is not difficult to find the asymptotic coefficients of the perturbative 
expansion to be given by 

f(o|o) _ V" ('^^ ^) ^ ~ Pkfl) ^ r (g - Pkfl -h + l) (fc|o) f,,.h-i , 
^ T{g- M V r (g - /3o,, -h + l) (o|,) 

As should be by now expected, we find that the coefficients of the pertubative expansion around 
the zero-instanton sector are given by an asymptotic double-sum expansion, valid for large values 
of g, over the coefficients of the perturbative expansions around (some of) the generalized multi- 
instanton sectors. The novelty in here, as compared to the one-parameter case of section 2, 
is that this expansion includes not only the coefficients of the physical instanton series Fg^^^\ 
associated with positive real part of the instanton action, but also the generalized coefficients 
Fg^^^\ associated with negative real part of the instanton action. In particular, the leading order 
of this zero-instanton asymptotic expansion is determined by the coefficients of the one-loop 
(generalized) one-instanton partition functions, but now up to two Stokes constants, namely 
S^'^ and S^^\. Higher loop contributions will arise as ^ corrections, while other multi-instanton 
contributions, with action itnA, will yield corrections suppressed as . 

Thus, what we have found in the present two-parameters transseries setting is that, such 
as in the one-parameter case, through the use of alien calculus and the bridge equations it 
is possible to include all multi-instanton sectors in the asymptotics of the perturbative zero- 
instanton sector. Furthermore, through essentially the same methods it is also straightforward 
to generalize this asymptotic result to all multi-instanton sectors. This is what we shall do 
next for the n-instanton sector, F'^'^'^^\z). Using the formulae for the Stokes automorphism 
in the directions = 0,7r, of $(n|o)) given in (4.28) and (4.40), we can easily find the related 
discontinuities in the said directions. As usual, F^'^l'') (z) has branch cuts in the Stokes directions 
corresponding to both positive and negative real axes in the Borel complex plane. Then, by 
means of (2.42) and (4.11), in particular the identification f("I'")(z) = e"*^"""^)^^ ^>(„|m)(2;), a 
lengthy but straighforward calculation leads to (it might be interesting for the reader to compare 
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this expression against its one-parameter counterpart, (2.51)) 




('gf^)^ ^ r (g + /3n,0 - Pn+kfi) ^ r (g + (3nfl - Pn+kfi - h) p(n+fc|0) 

5+/3n,0-/3n + fc,0 



{kAY 



El (9 + Pn,0 - Pn+fc,0 - /ij 77(™+fc|0) n A\h 



7ier(m,fc) <5jer(m,m.-£+l) 




^^)9+/3n.o-/3„+f-fe,^ ^ r (g + /3„,o - /3n+f-fc/) 



h=l 

where we have introduced 

S(n, j) = ((i + 1 - 5,) S'^^^t;^ + (n - 7,- + j + 1 - <5,) 5i%+^'^)) 6 (j + 1 - 5,) . (4.45) 

Recah that we have previously defined 5^^] = S^^ = Sq^ = 0, with £ > 0, and 70 = 0, 
Jo = 1, which are required to fully understand the formulae above. This result relates the 
coefficients of the perturbative expansion around the n-instanton sector with sums over the 
coefficients of the perturbative expansions around all other generalized multi-instanton sectors, 
in asymptotic expansions which hold for large g. All Stokes factors are needed to compute 
the general asymptotics of Fg'^^^\ whose computation is, in general, quite hard to do from first 
principles, but which may, nonetheless, be explored numerically in specific examples as shall be 
seen in great detail in the following sections. 

4.3 Resurgence of the String Genus Expansion 

The results we obtained in the previous subsections are rather general and do not take into 
account any symmetries or properties of the physical system that one might have started from. If 
we now specialize to the cases of interest in this work, models with a topological genus expansion 
such as topological strings, minimal strings or matrix models, then it is well known that the 
corresponding free energy in the zero-instanton sector will have a genus expansion as (1.1), i.e., 
an expansion in the closed string coupling g^, 

+00 

9'sF^'^'\9s;{U}) ^ ^g^^^Ff l°)(t.) = ^>(°l°)(g.;{iO)- (4.46) 

9=0 

This expansion resembles (4.10) if one sets z = 1/gs and assumes a tj dependence for the co- 
efficients Fg^^^\ti) in the asymptotic expansion (and similarly for the instanton action, A{ti)). 
These parameters, U, encode a possible dependence of the result on the 't Hooft moduli, as will 
be the case of matrix models. We also need to consider a string theoretic version of the ansatz 
(4.11) for the generalized multi-instanton free energies, this time around as an expansion in the 
open string coupling g^, 

(4.47) 



-29- 



Notice that in this expression we have further included an expansion in logarithmic powers of the 
open string coupling (up to some finite logarithmic power, knm) in order to account for resonant 
effects which will appear later in the Painleve I case and in the quartic matrix model, and which 
we have already mentioned at the beginning of this section^^ (see [13] as well, for the logarithmic 
terms). The integer /3nm will also be necessary in order to take into account possible different 
starting powers of our asymptotic expansions. For instance, in the case of the Painleve I equation 
we shall later find knm. = min(n,m) — m6nm and I3nm = (3{m + n) — [{knm + k)/2\j, where [•]: 
denotes the integer part of the argument, and where /3 = 1/2. We shall also make the assumption 
that the resonant effects do not appear in the n-instanton sector, that is knfl = = /co.m- Finally 
we will focus on the cases where both f3nm and knm are symmetrical in n, m. As we shall see 
later, all these assumptions turn out to be properties of string theoretic systems. 
Starting off with the zero-instanton sector, we have 

F^'^'\9s;{U}) = Y^9t''''°hr^^'\u). (4.48) 

If we compare this expansion with (4.46) above, one easily concludes that, in order to find a 
topological genus expansion, it must be the case that F^'^^^^^^'^^ti) = with /3^q = 0. Do notice 

that the free energy coefficients in the genus expansion (4.46) are given by Fg^^^^ = Fg^'^^''^^, 
which will naturally include both even and odd powers of the genus, g, as expected. Via Cauchy's 
theorem (2.42), now applied in the complex ^^-plane^^, one essentially recovers the result of the 
previous section for the F^^^^^^'^^ and, in particular, one finds for the asymptotics of F2i^_^l^^^ 



27ri 

(0)^ 

k 2- i-kAr^^-C h r(2£ + l-<l) ^ ^ ^ ^ 



^ 2^1 (^^)2m-/3- r(2^ + 1 - /3j,^i,) 

k 



( S\\ \ + 1 - /3i"[) +2^ r(2^ + 1 - - /3i^[) 

p{0|0)[0] 



The "genus expansion condition", that F^^^l^^^ = 0, now becomes equivalent to a set of relations 
between 1°^'°^, 1^'^'°^, 5^^ and We find, for each k and g, 

= i-iy+C l'^)™. (4.50) 

In the following sections we shall see in detail that by considering special properties of the systems 
we will address, such as 2d quantum gravity or the quartic matrix model, there are in fact more 
general relations between the i^j"''"'*''^', under exchange of n and m. Furthermore this will also 
allow us to find relations between sf^^ and S^l (and, in fact, relations between other Stokes 



^^Our discussion up to now solely focused on the "k — sector" of the logarithmic expansion. 

^^Notice that a blind application of Cauchy's theorem (2.42) in the ^s-variable leads to a large-order relation 
with an (incorrect) overall minus sign as compared to, e.g., (2.45). Instead, one should recall that the definition 
of the Stokes discontinuities in terms of the Stokes automorphism, (2.11), depends on what one means by left and 
right Borel resummations. Under a change of variables of the type x ^ 1/x these orientations change and so does 
the sign of the discontinuity — thus leading to the correct result. 
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constants) effectively reducing the number of independent Stokes constants needed to account 
for the large-order behavior of all multi-instanton sectors. 

The relation determined above can now be used to simplify the large-order behavior of the 
coefficients in the topological genus expansion (4.46), as 



ITT 
k=l 



which in fact, as just mentioned, reduced the number of Stokes constants one effectively needs 
to completely understand the asymptotics of the perturbative sector (comparing with the corre- 
sponding result in the previous subsection, (4.43), we see that this final expression is much closer 
to its one-parameter counterpart, (2.45)). Further notice that this expression coincides with the 
result in [13], at leading order in k, if one takes into account that in our case we are considering 
a genus expansion in the variable gs, instead of an expansion in z = I/qs as used in that paper. 

One can also use the string theoretic generalized multi-instanton expansion (4.47) to de- 
termine the large-order behavior of the physical n-instanton series F("|o)(z). This follows by 
applying Cauchy's theorem to the string coupling, Qs, and using the discontinuities for $(„|o)(z) 
determined in section 4.1. The novelty now is that we are further considering logarithmic power 
contributions to the asymptotic series of <I>(„|m)(2;). Thus, in order to obtain the large-order 
coefficients F^'^^^^^'^^ we shall apply Cauchy's theorem as before, but when making use of the 
expansion (4.47) new integrals will have to be addressed: 

r+oo r+oo 

Discontinuity at 6* = : / dx x"^"^ e~~ log"" x (-l)"" / dz z^^^ e~''^^ log" z, 

Jo z=\ Jo 

(4.52) 

Discontinuity at = TT : / dx x'S"^ e" log'' x ^ (-l)*" / dz z^?"^ e'^^Mog'' z. (4.53) 

Jo Jo 

The relevant quantity needed to perform these integral is the following Laplace transform 

r+oo / fi \ ^ r+cx) 



C[z9log'{z)]{s) = J dzz^e"'' log' z= (^j J dz z^ e 



'S z 





rr(5 + i)l 







S9+1 

and its analogous 9 = tt version 



6ro + e{r-l){Bs{g) + dgf ^ Bs{g) ] , (4.54) 



where^"^ 



C [z3 log^(-z)] {s) = ^-^^^ [Sro + e(r - 1) (Bs{g) + dg^^ B,{g)^ , (4.55) 



5,(o) = ^(a + l)-log(s), (4.56) 
Bs{a) = il^{a + 1) - log(-s) = Bs{a) - m. (4.57) 



■^■^In here ^p{z) — ^® digamma function; the logarithmic derivative of the gamma function. 
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Collecting all these results, one can now easily find the large-order behavior of Fg'^^^'^^'^^ (again, it 
might be interesting for the reader to compare this expression against the two-parameters case 
without logarithms, (4.44), or the one-parameter counterpart, (2.51)), 



+00 



F(»|o)[oi ^s-i" + k\ (sfV r(9 + p% - fti"! r(9 + - fii"|,,„ - h) |„, 



+00 I ^ fc ^ m / m 



, m\ 

k=l I ■m=l e=o 7ier(m,fc) Sj£r{m,m-e+l) \j=i 




E"'^^ '''' + ^i,0 /^li+e-k,t) r(g + f^l^^Q P\Xl-k,l ^) j^(n+l-k\e)[r] (uA\h 



X |(5,o + e(r - 1) (5fc^(a) + daj ' BkA{a)^ 



«=5+/3™o-/3!:|,_,„-/^-l 



(4.58) 



The quantity S(n,j) was previously defined in (4.45) as 



S(n, j) = ((j + 1 - S,) Si'^j;] + (n - 7, + J + 1 - <^,) 6 (j + 1 - <5,) . (4.59) 

One thing to notice is that, due to the logarithmic contributions appearing in the generalized 
multi-instanton expansion of <I>(„|m)(-2), the large-order behavior now includes contributions 
depending on the function Bs{a). The simplest possible contribution of this type in (4.58) is 

BkA {9 + /3i°'o - ^t:li-k,i -h-l)=ij[g + - /3M , -h)- log (kA) - ivr. (4.60) 

When g is very large {i.e., considering the large-order behavior) this expression may be expanded 
as 

BkA {9) c^^p{g)^logg-0{l/g), (4.61) 

where we made use of the asymptotic expansion of the digamma function around infinity. This 
shows that, in addition to the familiar 17! growth of the large-order coefficients, we now further 
find a large-order growth of the type g\ log g in the instanton sectors (which was also noticed in 
[13] for Painleve I) and generalizations thereof — as explicitly contained in (4.58). In particular, 
this is a leading growth when compared with gl and will be clearly visible at large order. 

As an application of the expression (4.58) above let us look at the case n = 1 and k = 2, that 
is, the 2-instantons contributions to i^j^'^^^*^', with particular focus on the ones which display 
a logarithmic behaviour. The contribution from the discontinuity at = is straightforward 
so we shall focus instead on the contributions arising from 6 = n. The sums in m and i have 
to be such that n + i — k > 0, which implies 2 > m > £ > 1. The cases with i = 1 will not 
have any logarithmic contributions as kn^i-k/ = ^0,1 = 0. Thus, the only case of interest is 
m = i = 2, which can have logarithmic contributions as long as kn+e-k/ = ^1,2 7^ 0. In this case 
7j G r(2,2) will have contributions from the Young diagrams EP and [J], and 5j € r(2, 1) will 
have only one contributing diagram, 0. Assuming that ki^2 = 1 (as will be the case of Painleve 
I) the 2-instantons contribution to i^j^l'^^t'^l becomes 

.,■10,101— .. (g.')^rte + MVM°i) ^ Fir"'h-2At 
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^^^i «[o] M] 2^ =h 7 M m V ^2^^^ + ^1.0 ~ ^1,2 - - iJ- 

2vri (_2^),+Mi-Mi ^o^^:^=l(5+M!.-/3S ^ 



2 - m 



The results obtained in this section can be extended to the generaUzed instanton series, such 
as, for example, the (n, l)~series. However, those generalizations yield extremely lengthy formu- 
lae. Consequently, we shall present those results only as they become needed in the following 
sections, and always in the specific form applicable to either of the particular cases of interest: 
the Painleve I equation and the quartic matrix model. 



5. Minimal Models and the Painleve I Equation 

We now want to apply the general theory of two-parameters resurgence developed in the previous 
section to some concrete examples appearing in string theory. The specific examples we have in 
mind are matrix models and minimal string theories, which, as is well known, are closely related: 
all minimal models can be obtained as double-scaling limits of matrix models. 

In this section we shall be mainly interested in the (2, 3) minimal string theory, which 
describes pure gravity in two dimensions, and whose free energy may be obtained from a solution 
of the Painleve I differential equation. Later, in section 6, we will turn to a similar resurgent 
treatment of the one-matrix model, where we shall see that, in the double-scaling limit, it exactly 
reproduces the minimal model results of this section. 

5.1 Minimal String Theory and the Double Scaling Limit 

Minimal models, labeled by two relatively prime integers, p and q, are among the simplest two- 
dimensional conformal field theories (CFT) and, starting with the seminal work of [59], they have 
been studied in great detail in the past (see, e.g., the excellent review [5]). 

Strictly speaking, the models we are interested in are not the minimal CFTs per se, but the 
string theories that they lead to. That is, we consider these models coupled to Liouville theory 
and ghosts and sum over all worldsheet topologies that the CFT can live on. The resulting genus 
expansion for the free energy is an asymptotic series, with the familiar large-order behavior 
~ (2(7)! [1], and it is the nonperturbative completion of this asymptotic series that we shall study. 
In particular, the simplest non-topological minimal string is the model with {p,q) = (2,3). It 
has a single primary operator, which after coupling to Liouville theory can be thought of as the 
worldsheet cosmological constant, and the central charge of the CFT is c = 0, meaning that the 
"target space" is a point: this minimal string theory describes pure gravity on the worldsheet. 

We shall discuss one-matrix models and their double-scaling limits in some detail later in 
section 6. For the moment, we only need one important result from the double-scaling analysis. 
The free energy F{z) of the minimal string theory depends on a single parameter, z, which is 
essentially the string coupling constant^"^. It is also convenient to define the function 

u{z) = -F"{z). (5.1) 

Then, from the double-scaling limit of the string equations of the matrix model one can show 
that the function u{z) satisfies a relatively simple ordinary differential equation which, for the 
(2, 3) minimal string describing two-dimensional pure gravity, is the famous Painleve I equation, 

u\z) - lu"{z) = z. (5.2) 
■^^More precisely, as we shall see in what follows, the c = closed string coupling constant equals z~^^^ . 
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One can solve this equation perturbatively in the string couphng constant, and the resulting 
asymptotic series gives the genus expansion of the (2, 3) minimal string free energy. What we 
are interested in here is to describe the full nonperturbative solution to this equation, in terms 
of a transseries. Since the differential equation is of second order, we expect such a solution to 
have two integration constants, and hence we should find a two-parameters transseries solution — 
exactly the type of transseries that we have discussed in the previous section. 

The construction of the two-parameters transseries solution to the Painleve I equation was 
started in [13], where the structure of the full instanton series and of the contribution with a 
single "generalized instanton" were found. Here, we complete this analysis by describing the 
structure of the full, general nonperturbative contributions to the solution. 

5.2 The Transseries Structure of Painleve I Solutions 

Let us now develop the transseries framework as applied to the Painleve I equation. 

Review of the One Parameter Transseries Solution 

As explained above, our aim is to solve the Painleve I equation, 

u\z) - lu"{z) = z, (5.3) 

in terms of a two-parameters transseries, where the perturbative parameter of the solution is the 
string coupling constant. As it turns out, in the minimal string, small string coupling corresponds 
to large z and hence the perturbative series in our solutions should be expansions around z = oo. 
It is well known, and one can easily check, that there is indeed an asymptotic series solution 
around z = oo in terms of the parameter z. It is given by 

^ (l - i .-'^ - ^ - .-'^ - . . ) . (..4) 

Note that, apart from the leading factor of z^^'^, this solution is a power series in z~^l'^ . This 
parameter is indeed known to be the coupling constant of the minimal string theory. However, 
^-5/2 -g j-^Q^ quite the perturbative parameter that we should choose for our transseries solution. 
The minimal string theory is a closed string theory, so indeed we expect its perturbative free 
energy to be a function of the closed string coupling constant. But nonperturbative effects 
in string theory, on the other hand, are associated to D-branes, and hence to open strings. As 
usual in string theory, the closed string coupling constant is the square of the open string coupling 
constant and, therefore, we may expect the nonperturbative contributions to the free energy to 
be expansions in 

X = z-^/\ (5.5) 

We shall later see that this is indeed the case. 

As a first step in finding a transseries solution to the Painleve I equation, one may now try 
to find a one-parameter transseries solution of the form 

+ 00 +00 

u{x) = cj^ e""^/^ x"^ x9, (5.6) 

n=0 9=0 

where x is expressed in terms of z by the relation above, and A and f3 are coefficients that still 
need to be determined. Plugging this ansatz back into the Painleve I equation (see, for example. 
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[15, 13]), one finds that a solution exists if one chooses 

A^.'-f, ,^\. (5.7, 

The same result could be obtained by writing the Painleve I equation in prepared form, (4.2), 
where one would find 

du 





n 8v^ 

, 5 . 



u{z) + ---. (5.8) 



For the "instanton action", A, there is a choice of sign. In the one-parameter transseries one 
usually chooses the positive sign, since with that choice the instanton factor exp(— A/x) is expo- 
nentially suppressed as expected. Doing this one finds, for example, the one-instanton correction 

Note that indeed we now find a series in the open string coupling x = whereas the purely 

perturbative part (5.4) of u{z) is a series in the closed string coupling = z~^l'^ . The coefficients 
in this expression can be determined recursively by plugging the ansatz (5.6) into the Painleve 
I equation. One finds that this determines all coefficients except the leading one, Its (non- 
zero) value can in fact be chosen arbitrarily without loss of generality, since we can rescale it by 
choosing the nonperturbative ambiguity ai. For now, we adopt a normalization where Uq^^ = 1. 

The Two Parameters Transseries Solution 

So far, we have only considered the positive sign choice for the instanton action A in (5.7). 
However, at the level of formal solutions, the negative sign choice is also required in order to 
obtain the most general solutions of the Painleve I equation, i.e., we should really apply the 
machinery developed in section 4 and solve the Painleve I equation using a too-parameters 
transseries. To do this, it is very convenient to change variables once again. Recall that the 
/3-parameter we found for the ansatz (5.6) equals f3 = 1/2. As we will see, the x-dependent 
prefactor in the two-parameters transseries will no longer be of the simple form x"^. It is 
therefore no longer convenient to take it outside the perturbative sum over g, as we did in (5.6). 
The analogue of x"^, on the other hand, will still be a half-integer power of x, so if we want to 
consider it as part of the perturbative series, it is more convenient to use the variable 

W = xl/2 = ^-5/8^ (5_10) 

Of course, up to a possible odd overall power in li;, we still expect all perturbative series to be 
expansions in the open string coupling constant, w'^, and we will find that this is indeed the case. 
Let us be a bit pedantic and stress this point once again, in order not to raise any confusions 
later on: the open string coupling constant is x = w"^ and we shall mostly work in the w variable. 

It is also useful for calculational purposes to scale away the overall power of z^^"^ in u{z), 
and set 

u{z) 



u{w) = -j^ . (5.11) 

Here, we slightly abuse notation; it would have been more precise to call the function on the 
left-hand side u{w), but to avoid writing too many hats we will stick to the above notation and 
simply remember whether we use the rescaled u or not by looking at the variable that we use. 
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It is now a simple exercise to rewrite the Painleve I equation in terms of the function u{w); 
one finds 

1 25 25 

u^iw) + u{w) - —w^ u'{w) - —w^ u"{w) = 1, (5.12) 

where we want to solve this equation using a two-parameters transseries ansatz, 

+ 00 +00 

u{w,ai,a2) = E E <^72"^e-("-™)^/"'' 1>(„|„)(t/;). (5.13) 

71=0 m=0 

Note that, as we have mentioned, we have not included a factor of w^""^ in the transseries 
expansion, but absorbed it in $(„|„)(?i;). This means that the leading coefficients in ^(^n\m) '^iH 
in general not multiply the constant term. Conversely, we can find back the analogue of the 
prefactor w^""' (as we will do below) by finding the first nonzero coefficient in $(„|^)(w). 

One may now be tempted to complete the ansatz above by assuming that (tt)) is a 

power series in w. However, an ansatz of this form turns out not to work, essentially since the 
Painleve I equation is a resonant equation (a property we have previously discussed in section 
4 and to which we shall come back in a moment). It turns out that, for a correct ansatz, one 
needs terms multiplying powers of log{w), a phenomenon first observed in [13]. In that paper, 
the authors calculated $(„|i)(t(;), and found that it had the general form 

+0O +CO 

9=0 9=0 

In fact, for n = 0, 1, the logarithmic terms are absent, but they are always present whenever n > 1. 
One may now wonder what the general form of $(n|m) is. From the li^-term in the Painleve I 
equation, one sees that is determined recursively in terms of products ^(^n-p\ni-q)^{p\q)- 

This means that, starting^^ at ^(4|2)) we can expect to encounter log^ w terms coming from terms 
such as '&(2|i)*^(2|i)- Extending this reasoning, we see that a natural ansatz for the general ^[n\m) 
is 

min(7i,m) +00 
k=0 9=0 

Our job now is to determine if a solution for all coefficients Ug^^'^^^^^ can be found. It is a tedious 
but straightforward exercise to plug the ansdtze (5.13) and (5.15) into the Painleve I equation 
(5.12) and, in this process, to find that the coefficients tt^"'"^^''^' must satisfy the relation 

n m g k ^ ^ 

5- <5o -^fj = E E £ E 4"'™^ ^'^^ -) [^^-^i _ (5.16) 

n=0 m=0 3=0 J=o 

+ 1 [m -n){g- 3) A^-j™)^ " ^ + 2) + 1) 4-^^"^'' " 
~{k + l){g-A) n("jf I'^+^l - ^ (5^^ - 16) {^9 - 24) u^^^^^'\ 



^We shall actually see below that, due to resonance, the log'^ w behaviour already sets in at $ 



(3|2)- 
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This relation is valid for any 4-tuple {n,m,k,g) if we assume that all non-existent coefficients — 
that is, the ones with k larger than min(n, m) and the ones with g < — are vanishing. The 
relation can be used to recursively determine y^"''"''''^' in terms of coefficients which have smaller 
n, m, k or g. A Mathematica notebook with the results is available from the authors upon request. 



The Consequences of Resonance 

In using the relation (5.16), one finds that something special happens whenever |n — m| = 1. In 
this case, the first term on the second line of the relation equals 



25^42 
96 



^HHW = _2^^H"')['=1, (5.17) 



where we inserted the explicit value (5.7) for A. However, this is not the only term multiplying 

Ug 

up to 



^{n\m)[k]_ sum in the first line of (5.16) also contains two terms with this factor, which add 



24°l°)Mu("l-)W=2n("l-)[^l, (5.18) 

where we read off the leading coefficient Wq'^''^^^'^^ = 1 from (5.4). Thus, we see that, whenever 
|n — m| = 1, the leading terms in the recursion formula cancel. This is precisely the phenomenon 
of resonance! The cancellation of the leading terms in itself is not a problem — it simply means 
that one should use our formula to determine u^^^^^''^ instead. However, it could potentially be 

a problem whenever y^'^^^t'^l does not exist — that is, when we try to determine the leading term 
in w for each perturbative series, given n, m, k. Here, two things can happen: 

1. The recursion relation may reduce to const = 0, in which case it cannot be satisfied. This 
is what happens if one does not include the correct logw terms. For example, if we would 
include no logarithmic terms at all, the recursion for n = 2, m = 1 would lead to such 
an inconsistency. Thus, resonance forces us to include the logarithmic terms. In a similar 
way, we will need log^ w terms starting at n = 3, m = 2. Note that above we have already 
argued that such terms must appear for n = 4, m = 2; now we find that we also need 
to include them in <I'(3|2), as we did in our ansatz. Only at n = m = 2 are the log^ w; 
terms absent. This pattern actually continues to higher m: ^(m|m) '^ill never contain any 
logarithmic terms; but the log'" w terms set in immediately at n = m + 1 due to resonance. 

2. The recursion relation may reduce to = 0. This is of course consistent, but it means that 
we have a leading coefficient which can be chosen arbitrarily. We already saw an example 
of this: u'^\ the leading coefficient of the one~instanton series, can have an arbitrary 
value due to the choice in the normalization of the nonperturbative ambiguity cji. In our 
two-parameters transseries terminology, this coefficient is now denoted u^"*^'*^^'*^^. The same 
thing now holds for n^'^'^'*^''', its value can be absorbed into 02- However, it turns out that 
the recursion relation allows for a whole lot more free parameters: for any m > 0, the 
coefficients ^^^^^^"^^^^^ and u^^^"^^^^^^^ are not fixed by our recursion relation. 

The second property above seems confusing at first sight. How can a two-parameters transseries, 
solving a second order differential equation, have infinitely many free parameters? The answer 
turns out to be that our ansatz still has a large degree of reparametrization symmetry. 
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Reparametrization Invariance 

Recall that our general transseries ansatz for the solution to the Painleve I equation has the form 

+00 +00 

"(^) = E E <^r^^2"'e-(-— )^/"'' $(„|„)(u;). (5.19) 

n=0 m=0 

It is important to note that the nonperturbative factor in each term only depends on the difference 
n — m. This means that when we make a (degree preserving) change of variables, 

+00 +00 
cri = Up {aia2y , = 0=2 E ^5 (^1^2)'' , (5.20) 

p=0 q=0 

with arbitrary coefficients ap, f3q, we will find a new expression with exactly the same nonper- 
turbative structure. Let us work this out in some detail. From the above change of variables, we 
get expansions of the form 

+00 +00 

^1 = ?r E (?i?2)^ , air = E (^1^2)^ , (5.21) 

where it is not too hard to find explicit formulae for the coefficients 7", 5^*, given by 

71 m 

7" = En"A. and 5T = Y.\{fi,^. (5.22) 

{A} i=l i=l 

In here, {A} and {//} are ordered partitions, where "ordered" means that, for example, we 
consider {0, 1,4} and {4, 1, 0} as different partitions of the integer 5. In the first sum, {A} runs 
over all ordered partitions of r with length n, and the analogous statement holds for the second 
sum. These formulae only hold for n, m > 0; for n = we have that 7q = 1 and all other 7^ = 0. 
The same thing of course holds for 5^. 

Inserting these results in (5.19), it follows 

+00 +00 +00 +00 

= E J^5^ar''+^ar^+^7."Ce-("-"^)^/'"' ^in\m){w). (5.23) 

n=0 771=0 r=0 s=0 

Changing the summation variables (n, m) to (n, fh) = [n + r + s,m + r + s), one obtains 

+00 +00 ro so 

= Y.Y. ^f?2^e-("--)^/-' x;E^r''"'^?"'"'^(--^i---)(^)- (5.24) 

n=0 m=0 r=0 s=0 

In this expression, ro = min(n, fh) and sq = min(n, fh) — r. In other words, r and s run over the 
triangle given by 

r > 0, s > 0, r + s < min(n, m). (5.25) 

Thus, we have found that, after reparametrization, u{'w) can be written in exactly the same form 
albeit in terms of new functions, 

^{n\m){w) =Y.ir''-'5T-'-'^^n~r~s\m-r-s){w). (5.26) 

r,s 
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Let us write out the first few of those: 

5(n|0) = "O ^HO), $(0|m) = ^o'^(Olm)- (5-27) 

Since we have already fixed the leading coefficients of ^(i\o) and <J>(o|i) to equal one, this means 
that we cannot freely choose and Pq: we have to set them equal to one as well. Using this, 
one finds for the next few 

$(n|i) = ^(n|i) + ai{n- 1) $(„_i|o), (5.28) 

$(l|m) = ^(l|m) +Pl{m- 1) $(o|m-l), (5.29) 

where n, m > 1. Thus, the ^(ri|i) are only defined up to additions of $(n-i|o)- can continue 
like this: after fixing qi and /3i it turns out that the free parameters 02 and show up for the 
first time in $(n|2) and <l'(2|m); and multiply possible additions of ^"(,1-210) ^(o|m-2)- 

This explains the fact that, in the previous subsection, we found that our recursive transseries 
solution had an infinite number of undetermined parameters. They are simply the parameters 
Op and f3q that determine the freedom in the parametrization of the coefficients ai and a2- One 
will find a unique two-parameters transseries solution to the Painleve I equation only after fixing 
these parameters by some sort of "gauge condition". 

Two Parameters Transseries: Results 

There is a rather natural condition^^ to fix the free parameters in our transseries ansatz. It turns 
out that the transseries component ^(rn+i\m) never has a constant term. We have also seen that 
^{i\o) starts at order w^, and we now know that we can use reparametrization invariance to add 
an arbitrary multiple of ^(i|o) to ^*{m+i|m)- Thus, one can tune the free parameter Um in such a 
way that the w^-teicm in $(m+i|m) vanishes. That is, one can fix half of the reparametrization 
invariance by simply setting 

^(m+l|m)[0] ^ > ^5_3Q^ 

In the exact same way, one can use the /3„-parameters to set 

^(n|n+l)[0] ^ > ^5_3^^ 

by adding the appropriate multiples of '&(o|i)- 

This fixing of the undetermined parameters is the last ingredient one needs in order to use 
the recursive formula (5.16) and solve for the entire transseries. Using a computer, this can be 
efficiently done up to n = m = 10 and 5 = 50 in a matter of minutes, and we have tabulated 
some of the $(„|^)(ti;) in appendix A. One thing the reader should note from those expressions 
is that the resulting functions are always, up to an overall factor, indeed expansions in the open 
string coupling constant x = w'^. 

The choices (5.30) and (5.31) simplify our results a lot, and sets many more of the leading 
coefficients to zero. Let us fix n, m and k, and ask ourselves what the lowest index g is for which 
^{n\m)[k] nonzero. We will call this index 2/3ifm (the factor of 2 is essentially due to the fact 
that we are now working with the w variable rather than x); it is the analogue of the f3nm in the 
general logarithm-free two-parameters transseries (4.11). Whereas in that case (3nm is usually 
of the form (n + m)/3 for a fixed /?, in the Painleve I case we find a more complicated structure. 

^^This condition is applied implicitly in the function ${2|i) reported in [13]. 
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Table 1: Values for 2/3„m (left) and 2/3„m (right). An asterisk in the second table means that there are 
no logarithmic terms in 



We tabulate 2/3j^m and 2/3^m in table 1. This table clearly has some structure and, in fact, it is 

\k] 

not too hard to find a general formula for 2f3nm- When n = m, none of the contributions have 
logarithms, and we have that 

2/3M = 2n. (5.32) 

For n 7^ m, it is easiest to write separate formulae for the cases n > m and m > n. When either 
n or m is smaller than k, we have no log*^ corrections. When n > m > k, one finds 

2/3]f], = n- A;+ (m + A; mod 2). (5.33) 

For m > n > k, the formula is the same, but with n and m interchanged. This can be summarized 
by defining, for all n and m, 



2/3W 

' nm. 



n-\- m 



(5.34) 



where ["Jj represents the integer part, and 

knm = min(n, m) - m bnm (5.35) 
is just the maximum power of the logarithm appearing in the expansion of ^(n\m)k^\ 
5.3 The String Genus Expansion Revisited 

We now have enough information to address the string genus expansion of the Painleve I solution, 
applying the general formulae previously obtained in section 4.3. Let us start by re-writing the 
asymptotic expansion for the $(n|m)) given in (4.47), as^'^ 

^nm -|-CxD fcnm -|-00 

$(nM(^)^Elog'^E^3"'"^'''^ =E^og'^ E^'^i"'^^ (5.36) 

k=0 g=0 k=0 g'=0 ^ 

i.e., as an expansion in w rather than as an expansion in x. Recall that our formulae in section 
4.3 were written in terms of the open string coupling gg = x = nP' , while in here we find it more 



^''Notice that the coefficients in the following just denote coefficients of a general transseries solution, 

in the abstract setting of section 4, and not the free energy of the (2, 3) model. We shall discuss the relation 
between the Painleve I solution and the (2, 3) free energy at the end of this section. 
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convenient to work directly in the tf-variable. Furthermore, in this expression it is understood 
that all the i^^"!"*^!'^^ with g odd vanish (in order to have an expansion in integer powers of x). 

2 

We can now directly compare with the expansion (5.15) for the Painleve I solution, and easily 
find that knm is given by (5.35) above, and^® 

^(n\m)[k] ^ 2^ F^"^^"^^^^^ ^ ^ 2*^f("I"')[''] f5 371 

9 I 2g 9 ' \ ■ I 

In particular, this implies that the ^x^"''"^''^' vanish for odd g. Moreover, the lowest index in g for 
which 'u^^"'!™')^ is non-zero, finm, can also be obtained via a comparison with the results of the 
previous section, being given by (5.34) above. 

In this way, we can rewrite the expansion of $(„|^)(x) for the Painleve I solution as 

kmn 1 h, -|-00 kjim i fc 

k=0 3=0 k=0 

with the knm and f3nm given earlier. It is now straightforward to apply the results of section 4.3 
to the current case. But, before that, let us address two important properties arising from the 
Painleve I recursion relations (5.16), i.e., from the physics of the (2,3) model, which will refine 
our results even further (also see appendix A). The first of these properties relates the coefficients 

\k] [01 

{n\m)^ at the fc-th logarithmic power, with the contribution without logarithms, as 

(n\in) ^1 I / {n—k\m—k) ^ ' 

This is a rather important relation; it amounts to saying that the logarithmic terms in (5.38) 
are actually not independent of each other, as their coefficients are all related to the coefficients 
of the logarithm-free term. In other words, these logarithmic contributions simply amount to a 
useful arrangement of the resonant transseries solution. The previous relation can be written in 
terms of the tt^"'!™^^^] by noting that /3nm = Pn^k m-k ^'^^ i^ViS 

^in\mm = i. |^ ^(^ ' ^{n-k\m-km . (5.40) 

The second property we shall be using relates the different it^^'^^t'^l under interchange of 
n ■v^ m. This relation can be found in appendix A and is given by^^ 

^{n\m)[k] ^ ^_^y+i3y^l_(j^+rn)/2 ^{rn\n)[k] ^ ^_-^y-[{k„^m+k) /2\ ^{rn\n)[k] ^ {^.41) 



■^^More precisely, the relation between and ^'"l'"'!''' is given by 

2*: p(n|m)[fc] _ (n|m)[fc] 
S g' ' 



where = 2 + and g starts at 0. To write the expansion of (x) we performed a shift on the variable 

g such that u^"'™)'''' yt^l™-)!*^! where now the expansion starts at ^(^"l'""*' a;''"™. 

^^Recall that we previously performed the change ^i^^'™"''' — >■ itjg''"''*', with 5' = 2 + f^nX^^. 
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Note that the exponent of (—1) in the above expression is always an integer. In the case where 
n = m (and consequently k = 0) we find that this relation returns 

Consequentially, the (n|n)-instanton series will always have a topological genus expansion 

+ 00 +00 

<^(n|n)(x) ^ X- ^Sj"!'^) x^^ = X- )M x'^. (5.43) 

9=0 9=0 

Looking back at the zero-instanton series from section 4.3, we have the genus expansion 

+00 

9=0 

where x = Qs and where the large-order behavior follows from 

57(0|0) _ ,,(o|o)[o] _ ('^^ ^) ^(^g - ^ki)) r(2ff -h- /ggl)) (^io)p] 

One can also write large-order formulae for the asymptotics of the Painleve I multi-instanton 
coefficients in the current language. This amounts to inserting these coefficients, written as 
(5.37), back in (4.58). The condition Ug^gm+i) = was studied in equation (4.50), which, when 
applied to the present case and by further using (5.41), yields 

ymy^Wom ^ (5i°))'n(f )M = (-l)<i (5.46) 

This immediately implies the following relation between sf'^ and 5"^^^ 

5i°) = (-l)lsi°}, (5.47) 

which coincides with a result found in [13]. 

The aforementioned properties (5.41) and (5.42) for the Painleve I coefficients can, in prin- 
ciple, allow us to find many possible relations between the Stokes coefficients sj^^ and S'^'"\ We 
will present one more such example in the following, with the study of the (n, l)-instanton series. 
First, using the same tools as in section 4.1, we can find the Stokes automorphism for the series 
^{n\i)i^)j both at 6* = 0, 

60'J>(n|l)(^) = <^>(n|l)(^) + Ern " "''^'^ ^^-^^^ 

fc=l ^ ^ 

V^(°)f<D. , I ( ^^^-^) c(0) I fc(2n + fc-l) (1) (o)\ , . 

and at = TT, 

©,$(„ii)(z) = $(„|i)(z) + EE^EE E n J') • -^(n-i-fc+^K), 

fc=lm=l ■ i=0 7i6r{m,fc) (5j6r(m,m-^+2) j=l 

(5.49) 
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where this time around we find'^'^ 

S(^) in, j) = ({j + 2-5,) + (n - 1 - 7,- + j + 2 - 5,) s'^^'^'t'^A Q{j + 2-5,). (5.50) 



With these results in hand, one can use the asymptotic expansion (5.38) and Cauchy's 
theorem to obtain the large-order behavior of the coefficients ^2^'^''^'^, with n > 1 and r = 0, 1. 
In order to simplify this calculation we shall now make use of the property (5.40), relating the 
logarithmic sectors, and thus write the expansion of <l>(„|i)(a;) as 

+ 00 ^ +00 

= ^Si)(-) + log^ • ^Hi|0)(^)- (5-52) 

At this stage, we already know the asymptotic behavior of u'^ g^j^^j j^q-^^ want to determine 

the asymptotics of Uj^'^^^*^'. Furthermore, we know the discontinuities of <l>(„|i)(x) given the 

Stokes automorphisms above. Thus, applying the Cauchy formula to the function <I>'^|-|^-|(a;), and 
making use of the relation (5.52), one obtains 




dw Disced ^n\i){w) 2(n-l) f'""^ du; Disc 



^ gt^^\Jo 27ri w-x Jo 27ri ^ w - x 

(5.53) 

The asymptotics of u^^^^^^'^ will have a contribution from each of these integrals, except in the 
case when n = 1, where only the first integral is present. In this case we have already seen that 
^(i|i) (^) "^ill have a genus expansion, as a consequence of the condition that ^*2(2rn+i) ~ ^' Solving 
this condition, using (5.41), we find more relations between the Stokes coefficients. Summarizing, 
these relations are 

5W = (-l)i5it (5.54) 
5f = 8% (5.55) 

As discussed before, requiring a genus expansion of <I)(„|,„)(2;) for n > 1, which is equivalent to 

setting u^2{2m^i^ = 0, will then yield a tower of relations between different Stokes coefficients, 
effectively reducing the number of independent coefficients needed to account for both the full 
multi-instanton asymptotics as well as any possible Stokes transition one might wish to consider. 

5.4 Resurgence of Instantons in Minimal Strings 

The recursion formula (5.16) provides us with a tool to calculate the two-parameters transseries 
solution of the Painleve I equation, to arbitrary precision. In particular, this allows us to do 
high-precision tests of the resurgent properties that were discussed in general terms in section 4 
and that were discussed in the specific Painleve I case in the preceding paragraphs. 



^"Comparing against the (n, 0) case, (4.45), the reader may want to guess a solution for the arbitrary (n, m) 
instanton series. 
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Resurgence of the Perturbative Series 

One of the main new phenomena that our resurgence analysis uncovers is the fact that the large- 
order behavior of transseries coefficients is itself subject to nonperturbative corrections. This 
phenomenon is already present in the simplest case: the large-order behavior of Ug^^^''^^\ the 
zero-instanton, perturbative expansion coefficients of the Painleve I transseries solution. 

Recall that, in our normalizations, these coefficients are only nonzero when g is a multiple 
of four. To avoid writing unnecesary factors, let us rescale 

i^-^^^'^ ^,{o|o)[o] 
r {29 - 

We can then write the large-order formula (4.51) as 

+00 -n /'o„ I, l\ 



U4g 



V^„(i|o)[o] fc r(2g-/t-^) q(o) (2|o)[o] „h-29+i Ah+h r(2g- /t- 1) 



h=o ^ ^ h=o 



The ratios of gamma functions in this expression should be thought of as perturbative 1/g 
expansions. For example, we can rewrite the ratio of gamma functions in the first sum of the 
first line above as 

r(2,-i) - l\2g-k-l- 2^' + 2'^^2 9 + ■ (5-59) 

In this way, we can define these ratios as (possibly asymptotic) series for any values of g and 
h. In particular, this allows us to work with expressions such as, for instance, the factor of 
r(2gi — h — 1) in the second sum in (5.58), even when 2g — h — 1 is a negative integer for which 
the actual gamma function would have had a pole. 

Thus, the ffist sum in (5.58) gives a purely perturbative description of the large g behavior 
of the u^g coefficients, as a series in 1/g. This perturbative large-order series has been studied 
in detail in [15, 13] and was found to give correct results up to high precision. What we see now 
is that, nevertheless, the perturbative large-order behavior is not the full story. For example, 
the second sum in (5.58) contains further corrections that come with a factor 2~^^, and therefore 
are invisible in a perturbative study. The sum in the second line of (5.58) gives 3~^^ corrections, 
and so on; one keeps finding subleading multi-instanton corrections in this way. 

The question is: can we actually see those nonperturbative corrections to the large-order 
behavior? It should be intuitively clear that in order to see an effect as small as 2~^^ at large 
g, we first need to subtract the leading perturbative series to very high order. Here one actually 
runs into a problem since the perturbative series in 1/g, the first sum appearing in (5.58), is not 
convergent — it is an asymptotic series. This should not come as a great surprise: we know that 
the presence of nonperturbative effects in a quantity is closely related to the nonconvergence of 
its perturbation series. This phenomenon pops up again in the large-order formula. 

Optimal Truncation 

The simplest way to deal with asymptotic series is to do a so-called optimal truncation: one 
simply sums the terms in the series for as long as their absolute value decreases, and cuts off the 
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Figure 2: The logj^g of the absolute value of the first 200 coefhcients in the 1 /^-expansion associated to 
the first sum appearing in (5.58), for the case where g ~ 30. 



sum at this point. As an example, let us look at the case where g = 30. In figure 2, we have 
plotted the log^g of the absolute value of the first 200 terms in the l/(7-expansion associated 
to the first sum in (5.58). The smallest term in the series occurs at order g~'^^ and equals, 
approximately, —2.8 x 10~^^. We see from the figure that, after this term, the terms in the 
asymptotic expansion start growing again. Thus, optimal truncation instructs us to cut off the 
sum after the order term. We expect that the size of the final term gives a good indication 
of the precision of the calculation. This is indeed true: one finds that 

U4.30 = 0.9978832395689425456292 . . . , (5.60) 
u1%o = 0.9978832395689425456257 . . . , (5.61) 

where "ot" stands for "optimal truncation", and, in the first line, we have calculated the exact 
value using (5.57). Thus, we get the correct result within an error of 3.5 X 10-21— indeed of the 
order of magnitude of the last term in the optimally truncated sum. 

The problem with this method is that it is only barely sufficient to distinguish the 2"^^ 
effects associated to the second sum appearing in (5.58). For our example value of g = 30, the 
leading term in this sum is 



(0) „(2|0)[0] 

2.33... X 10^2" i. (5.62) 



^T^^T^ r(59) _ ^ ,^ _^_,o. 

259 p /'119 



2 



We see that this leading term in the 2"^^ corrections is roughly of the same order of magnitude 
as the error in the optimal truncation'^^ . In other words, this term is only just within the 
"resolution" that optimal truncation allows us, and any 1/g corrections to it (let alone the 3"^^ 
corrections) will be completely washed out by the error due to optimal truncation. This is 
not just an unlucky coincidence: one can show using general arguments (see, e.g., [56]) that 
optimal truncation always leads to an error which is of the same order of magnitude as the first 
nonperturbative contribution. 



"^We will soon also explain the perhaps surprising fact that this term is imaginary. 



-45- 



Borel Pade Approximation 

Since optimal truncation is not powerful enough, we need a better method to approximate the 
asymptotic series associated to the first sum appearing in (5.58). That is, we actually need to 
resum this series. Of course we already know of a very powerful method to resum asymptotic 
series: the method of Borel resummation, discussed at length in section 2. To employ this 
method, we would in principle need to find the Borel transform (2.4) of the first sum in (5.58), 
and then do the Laplace transform (2.5) that inverts the Borel transform. The problem with 
this procedure is that we only have a recursive definition of the coefficients in the asymptotic 
series and, as a result, it seems impossible to find an exact expression for the Borel transform. 
Note that approximating the Borel transform by a Taylor series will not do: the inverse Borel 
transform will then simply give back our original divergent series. 

The solution to this problem lies in the method of Borel-Pade approximations. Let us write 
the 1/g expansion associated to the first sum of (5.58) as 

+0O 



P{g) = J2^ng-''. (5.63) 

n=0 

The Borel transform (2.4) of this asymptotic series is 

+ 00 



n! 

n=0 



One can check that the a„ grow factorially with n, so that this new series has a finite radius of 
convergence. However, we can only calculate the a„ recursively, so in numerical calculations we 
will actually have to cut off the above sum at some large order. For convenience, we choose this 
order to be an even number, 2A^, 



2N 

a. 



n=0 

Instead of directly performing the inverse Borel transform (which, as we mentioned, would give 
back the original asymptotic series), we now further approximate this function by an order N 
Pade approximan^"^ 

That is, the degree 2N polynomial in 1/g is replaced by a rational function which is the ratio 
of two degree N polynomials in 1/g. The coefficients in this approximation are chosen in such 
a way that the first 2 + 1 terms in a 1/ (/-expansion of i3[^l[p](s) reproduce i3[P](s). When 
one furthermore chooses cq = 1, to remove the invariance under homogeneous rescalings of all 
coefficients, this requirement can be shown to lead to a unique set of (6„,Cn). There exist fast 
algorithms to determine Pade approximants; for instance in Mathematica such an algorithm is 
implemented under the name PadeApproximant. 

^^More precisely, this is the order {N,N) Pade approximant. One could, in principle, choose different orders of 
for the numerator and the denominator, but in numerical approximations this so-called diagonal choice often 
leads to the best results. As we shall see, in our case it indeed leads to very precise numerics. 
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The virtue of replacing the polyomial by this rational function is that, for small 1/g, both 
functions look very similar, but for large the rational function approaches a constant and is 
therefore much better behaved. As a result, one can now calculate the inverse Borel transform, 
or Borel resummation, 

SrP{9)= / ds^[^][P](s/5)e-^ (5.67) 

Contrary to our original asymptotic series, this result will indeed converge in the limit where 
N ^ oo. Note that the subject of Borel-Pade approximations has been studied intensively from 
a mathematical point of view, and has been applied to several physical problems in the past — the 
reader can find further details, for example, in [60, 12]. 

One thing one needs to be careful about when doing a Borel-Pade approximation is that 
the rational function i3[^l[-P](s) will, in general, have poles on the positive s-axis, making the 
integral (5.67) ill-defined. This problem is precisely the same as the one we encountered earlier 
for the ordinary Borel resummation in section 2, and we now know how to solve it: instead 
of integrating (5.67) along the real s-axis, we need to integrate around the poles using a +ie 
prescription^^. As a result, the resummed approximation S^^^P{g) will no longer be purely real, 
but will have a small imaginary part. For example, using a Borel-Pade approximation for the 
first sum appearing in (5.58), in our example case of g = 30, we find the value 

u^^^^^ = 0.9978832395689425456292 ... - 2.26 . . . x 10"^° i, (5.68) 

where the (1) indicates that we only resummed the first sum in (5.58). Comparing this to (5.60) 
and (5.62), we notice two very important facts. First of all, the Borel-Pade approximation 
indeed gives better results than optimal truncation: at the precision to which we are presently 
calculating, the real part of the above expression exactly reproduces (5.60). Moreover, shedding 
light on our previous evaluation, the imaginary part of the above result is of the same order of 
magnitude as (5.62), albeit of opposite sign. That is, it is largely canceled by the leading 2~^^ 
term in (5.58) which, as we now understand, indeed needs to be imaginary. The fact that the 
cancellation is not precise is because in (5.62) we only calculated the leading term in the 2~^^ 
corrections; adding further terms will give more precise results. 

Testing the 2~^f Corrections using Richardson Transforms 

We shall see in a moment how incredibly precise these results can be made, but first we want 
to perform an additional test on the validity of our large-order formula (5.58). Traditionally 
(see [15, 12, 13] for many examples), large-order formulae are tested as follows. One finds a 
g-dependent quantity, Xg, such that the ratio 



Rg = ^ (5.69) 



approaches a certain coefficient, i?oo, at large g, and such that the corrections to this large-order 
value take, at least to a good approximation, the form of a 1/g expansion. One then calculates 
Rg for a sequence of low values of g, and finds Roo using the numerical method of Richardson 
transforms (see, e.g., [60, 15]). 



^■^This sign is a matter of convention; integrating using a — ie prescription will lead to the same large-order 
formulae, but with the imaginary Stokes constant S^"' replaced by — S'J"''. 
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Figure 3: The first 100 values of tlie sequence Rg, and tfie first tliree Ricliardson transforms of this 
sequence. The sequences accurately approach the numerical value 4, as expected. 



Let us work this out for a particular example. Since the perturbative large-order behavior 
associated to the first sum in (5.58) has already been tested extensively in [15, 13], the first 
interesting thing we can test is whether the corrections to this term scale as 2~^^ as g ^ oo. To 
this end, let us define 



U4g 



"■4q 



sf^A\T{2g-l) 
3 X 29r (2^ - 



+ 0{l/g) 



(5.70) 



where the right hand side is the result we expect if our formula (5.58) is correct. From this, we 
find the expectation that the ratio Rg should go like 



Rg = A + d/g) 



(5.71) 



as g — 7- oo. To check this expectation, we have plotted the values of Rg for the first 100 values of g 
in figure 3 (top blue line). We see that the sequence Rg indeed approaches the numerical value 4, 
albeit slowly. To increase convergence, we can remove \/g effects by calculating the Richardson 
transforms of this sequence. This method is explained in some detail in [60], for example. The 
figure shows, from top to bottom, the first three Richardson transforms of the sequence Rg. We 
see that the sequences accurately approach the value i?oo = 4. The best convergence happens 
after seven Richardson transforms, and in this way we find a limiting value of 



Rc, 



4.000000000038. 



(5.72) 



The fact that we numerically find the expected answer up to one part in 10^^ gives us a lot of 
confidence that the 2~^^ behavior in our large order formula (5.58) is correct. 

We could continue and define a new Xg which tests the prefactor of the 2~^^ corrections, 
and then the subleading terms in l/g, and so on. We will not do this here, since we shall now 
see that there are other tests that check the coefficients in our formula to even higher precision. 

Direct Numerical Evaluation 

One may jump ahead and wonder: can we also see the 3~^^ corrections in our large-order 
formulae numerically, and perhaps even go beyond those? It is clear what needs to be done for 
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this: to see 3~^^ effects, we not only need to resum the perturbative asymptotic series, but also 
the asymptotic series multiplying the factor 2"^^ in the first line of (5.58). As we shall now see, 
the method of Borel-Pade approximations is powerful enough for this to be done. 

Let us begin by again focusing on an example. To make sure the numbers that follow fit on 
a single line, let us now look at the case g = 10. We will denote the Borel-Pade resummation of 
the first n distinct sums in (5.58) by u£io- ^ numerical evaluation of the Fade approximant and 
the consecutive Laplace transform gives the following results: 

n4.io ~ 0.995695607481681532429, (5.73) 

S4.10 - S^io ~ 0.000000000000249496840 + 0.000000041490689176523 i, (5.74) 

S4.10 - uflo ~ -0.000000000000498993666 + 0.000000000000000063033 i, (5.75) 

S4.10 - Sjio ~ -0.000000000000000000043 - 0.000000000000000063033 i. (5.76) 

From these numbers, we learn the following. First of all, we see again that already the leading 
Borel-Pade approximant u^^iq gives a very good approximation to the actual value ■^4.10. It is 
off by a term of order 10~^ in the imaginary direction, and only by a term of order 10^^'^ in 
the real direction. This imaginary error is then canceled to very high precision by the order 
2~^3 terms, leaving an imaginary error of order 10~^^. Meanwhile, the real error is not further 
corrected at this level. The reason for this last fact is that the real error in both the perturbative 
terms and in the order 2"^^' terms come from 3~^^ effects, and are therefore of the same order of 
magnitude. We see even more: they are not only of the same order of magnitude, but actually 
related by a simple rational factor: the real error in the 2~^^ terms is —3 times the real error 
in the perturbative terms, thus giving the overall real error a factor of —2, as seen above. That 
these errors are so simply related could have been anticipated: both come from the 3-instantons 
series in the transseries solution to the Painleve I equation. 

The remaining real error is then canceled to order lO"^'' by the 3-instantons effects and, 
at this order, the imaginary error stays of the same magnitude, again being related by a simple 
rational factor to the imaginary error at the previous level. One can continue like this: the 
remaining imaginary error will now be canceled by 4"^^ effects, the next improvement in the real 
error will occur at order 5~^^, and so on^^. 

~{n) 

To see how well this method works, we show in figure 4 the precision of u\g for g ranging 
from 2 to 30 and n ranging from 1 to 6 (i.e., we have tested our results up to six instantons). 
To obtain these numbers, we have done the appropriate Borel-Pade resummations up to orders 
(200, 180, 160, 120, 80, 80) for the (1, 2-2^, 3-2^, 4-29, 5-29^ g-^s) corrections, respectively. Along 
the vertical axis, we have plotted the precision, which is defined as 



log 



10 



~{n) 



(5.77) 



~{n) 

that is, the number of decimal places to which u\g gives the correct result. 

^*Since we know that the precise value of ung is real, we could actually have ignored all imaginary errors. This 
reduces the number of Borel-Pade approximations that one needs to make by a factor of two. As the simplest 
example, we could take the real part of u^g — u''^^ , and find a result which is correct up to 3^^^ corrections instead 
of the expected 2~^^ corrections. In the explicit calculations, however, to make sure that our methods are correct 
in more general cases, we have not used this simplification. 
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Figure 4: The precision of u\ for g ranging from 2 to 30 and n ranging from 1 to 6. 



Thus, we see that our large-order formulae lead to extremely accurate results. We saw before 
that, for g = 30, optimal truncation of the perturbative 1/g large-order series gave the correct 
result up to approximately 20 decimal places. Now we see that using Borel-Pade approximants 
the nonperturbative n~^^ effects play a crucial role in getting higher precision, and that by 
including up to 6~^^ corrections we can get results that are correct up to 60 decimal places. 

Two final remarks about these results are in order. First of all, even though we are speaking 
of "large-order behavior", we see from figure 4 that already at (7 = 2 we get results which are 
accurate up to 10 decimal places. There is still, however, a limit to this procedure. The reason 
for this is that, in our normalization, the n-instanton series at genus g comes with a factor of 



The leading [h = 0) terms for n = 8 will therefore blow up when g = 2. This is also the reason 
why we have not included 5 = 1 in our graph: there, already the n = 4 and n = 6 contributions 
blow up. It would be interesting to know if this is indeed a fundamental problem or whether it 
is simply a matter of normalizations, and can be solved in a similar way to how we circumvented 
the analogous gamma function singularities for nonleading values of h. 

A second remark is that this test can be viewed as a much more accurate test of certain 
coefficients than the traditional tests using Richardson transforms. For example, at (7 = 30 we 
have seen that the 2"^^ effects set in at order 10"^*^. However, we have now checked formulae 
for (7 = 30 up to order 10"^*^. If the base 2 in 2~^^ would have had an error, 5, this would have 
shifted the 2~^^ effects to effects of order 



so from the fact that we get correct results up to order 10~ , we see that 6 can be no larger 
than of order 10"^". This is a huge improvement compared to the accuracy of order lO"'^^ that 
we found using Richardson transforms. 

This does not mean that the method of Richardson transforms has become useless. Note 
that, in the above tests, we have essentially "reversed the burden of proof" : we have assumed that 



T{2g-h-f) 



(5.78) 





(5.79) 
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the n~^^ corrections (coming from higher instanton coefficients in the transseries) were correct, 
and checked these against the known perturbative coefficients. In the method of Richardson 
transforms, one starts from the known perturbative coefficients, and reproduces the expected 
coefficients in the n~^^ corrections. Whereas the ffist method is more powerful as a test, in 
practical cases one is more likely to know the perturbative coefficients in a transseries than to 
know all nonperturbative coefficients, as we do in this example. Thus, in those cases, Richard- 
son transforms can be used to learn something about the nonperturbative data, starting from 
perturbative data. This approach can be useful for example when studying topological string 
theories, where detailed nonperturbative information is often unknown. 

Resurgence of the (n|m) Instanton Series 

Now that we have gained some confidence in our resurgent techniques from studying the pertur- 
bative series ^(o|o)' a-PPly these techniques to the n-instanton perturbative series ^^^n\o) 

\k] 

and, more generally, to the generalized instanton series 

A new phenomenon occurs here: the large-order behavior of the series coefficients, u^"'™'^^'^', 
no longer depends only on the single Stokes constant sf'\ and further Stokes constants will 
appear. For example, applying (4.58) to the one-instanton series one finds that, up to order 2~^, 
its large-order behavior has the following six contributions 



^(l|0)[0] 



^Als^jmm (-i)^5i°^ ^ a|i)[o] r(g-2/.-i) 

29+1 - 27ri Z^''2h+2 • h^l + 27ri Z^''4h+2 • .g-2h-i + 

h=0 ^ ^ h=0 ^ ^ 

(310)10] r (,-/.-!) (-^)-KT f^.^(2,)to] ri,-h-i) 

i-^y{sn^fsnm)^- ,,,, n^-h) 



2vri ^ 2h+i ^2Ay-h 

Several facts should be noted about this expansion: 

• The first two sums determine the perturbative large-order behavior of the one-instanton 
coefficients, as a series in 1/g. In the zero-instanton case, (5.58), we saw that this behavior 
was determined by the next instanton series — in that case, the one-instanton series. In the 
first sum above, we see this "forward resurgence" again: the large-order behavior of the 
one-instanton series is partly determined by the two-instantons series. However, we see 
from the second sum that there is also "sideways resurgence" : the large-order behavior of 
the one-instanton series also depends on the (1|1) generalized instanton coefficients. Thus, 
even though the physical interpretation of these generalized sectors is somewhat mysterious, 
they do influence the physical instanton sectors in a very important way. 

• Even though it does not happen in the above example, from the structure (4.19) of alien 
derivatives, one can easily see that, in general, also "backward resurgence" will occur. For 
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example, the large-order formulae for the two-instantons series will contain contributions 
coming from the previous, one-instanton series. Thus, already at the perturbative level 
in 1/(7, we find a very intricate pattern of relations between the different generalized in- 
stanton series. This pattern gets even more intricate at higher nonperturbative orders. 
For example, in the last four sums of the above formula, we see that, at order 2~^, the 
large-order behavior of the one-instanton series is determined by the 3-instanton series, 
by the generalized (2|l)-instanton series (including its logarithmic contributions ti2ft'+i^^)' 
and even recursively by the 1-instanton series itself^^. 

In the last sum above, two new Stokes constants appear: S^^ S^l (recall that S^l = 
iSi'\ so it is not a new constant). The new constants appear in the combination 

T = s^^l+^-S^'ls% (5.81) 

so that by matching the right-hand side of (5.80) to the left-hand side for large values of g, 
we can determine T up to corrections coming from terms. Note that to do this, we need 
to calculate Borel-Pade approximations to the infinite sums in (5.80). This is a procedure 
which takes some (computer) time, but other than that is relatively straightforward. 

Finally, recall that i?2A(5 — h) = tp {g — h + 1) — log {2A) — ivr, where is the digamma 
function. At large g the digamma function has the asymptotic expansion 

1 B 

*(.)=log(.)---E5j;,. (5.82) 

n=l 

where in here B2n stands for the Bernoulli numbers. The leading term in this expansion 
implies that, at large order, B2a{9 — /i) ~ fogs', i-^-, we find at the 2-instantons level a 
growth of type gllogg, leading as compared to g\. In calculating T, the easiest way to 
deal with this behavior is to gather all terms multiplying \ogg, divide out the log 5, do a 
Borel-Pade approximation and then multiply with log 5 again. The further terms coming 
from the above asymptotic expansion can then be treated as before, using Borel-Pade 
approximation to resum all of them directly. 

Carrying out the Borel-Pade approximations, we have found that 

T = -0.90573009110532780736 .... (5.83) 

The precision of this number can be determined as follows. Note that, for any g, we can determine 
Tg from (5.80), and we expect the result to become better as g becomes larger. In fact, the true 
value of T should be Too. We have calculated Tg for values of 5 up to 151. In figure 5, we plot 
the number of decimal places to which Tg agrees with T151. When g <C 151, this is essentially the 
number of decimal places to which Tg agrees with Too- At g ~ 151, this is no longer true, since 
we are not really comparing with Too, but with T151. One finds that for g <C 151, the precision 
increases linearly, so by extrapolating this linear behavior, we find the expected precision of T151, 
which in this case is a bit more than 20 decimal places. 



''''This behavior is a consequence of the symmetries of the problem: it is really '"2^+1'^^ that appears in the last 

(1|0)[0 
2h + l 



sum of the large-order formula, but we have used equation (A. 21) to rewrite these coefficients in terms of u''^'"""' 
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Figure 5: The precision of Tg with respect to T151, and the resulting Unear extrapolation to g = 151. 



Of course, it is really the separate values of S^^l and S^l that we want to calculate, and not 
the value of the particular linear combination T. This can be achieved, for instance, by looking 
at the large-order behavior of the series n^^'^^'^' which, at the perturbative level, only depends 
on S^l and on the known constant In exactly the same way, we can then calculate S^l 

and from that constant and T determine S^l- Applying this procedure to several generalized 
instanton series, we have calculated a series of Stokes coefficients that are tabulated in table 2. 
Note that, to calculate these numbers, we have tested the resurgence of several of the generalized 
instanton series up to three instantons. In this table we have also indicated the number of 
decimal places to which we have calculated the answer. In the case of sf^\ an analytic answer 
is known — see for example [36, 17, 21, 15] for derivations. One has 

= (5.84) 



For 5f , we have actually listed more decimal places than we have calculated; we will see in a 
moment why we are able to conjecture some further digits. For readability, we have only listed 
about 20 decimal places for each Stokes constant; the authors will of course provide further 
data to the interested reader, upon request. In the table, we also list which type of large-order 
behavior the Stokes constants determine. For example, the constant S^^ appears for the first 
time in the large-order expansion of '^'|2|o)' ^'^^'^^ multiplies the terms of order 3~^. Apart 
from sf^\ and to the best of our knowledge, the only other number in this table which has been 
calculated before is S\ [13]. This number, called S-i in equation (5.38) of that paper, was 
calculated numerically in there up to 19 decimal places, and our result agrees with this up to 
17 decimal places {i.e., the final 2 decimal places that are reported in [13] are incorrect — this is 
possibly a result of the onset of 2~^ effects that were not taken into account in that paper). 

Note that we have only listed Stokes constants S^^^ and S'^"^ with i > 0. The reason is that 
all of these are purely imaginary, but, from them, one can then easily calculate the corresponding 
set of Stokes constants with i < using relations such as (5.54) and the ones that follow it. 

More interestingly, we find that the Stokes constants with i > also satisfy several (at this 
stage, unexpected) relations amongst themselves. The first thing one notices is that it seems 
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-0.371257624642845568... i 


00 
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(0|0) 


1-9 


c(0) 


0.500000000000000000... 1 
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<p , , , 
(1|0) 
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-0.897849124725732240... i 


13 


CP , , , 
(2|0) 
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0(1) 


-4.879253817220057751... i 


81 
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9.8od8759804878dz735... 1 


19 


"^(211) 
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-22.825711248125715287... i 


36 


(2|2) 


1-9 




2.439626908610028875... i 


112 


20) 


1-9 




15.217140832083810191... i 


108 


<!> ° . 

(3|1) 


1-9 




45.334204678679729580... i 


108 


(4|2) 


1-9 



Table 2: The Stokes constants that we have calculated. The third column gives the number of decimal 
places to which the answer is explicitly calculated. The fourth column lists the generalized instanton 
series for which the Stokes constant appears for the first time, and the fifth column lists what type of 
large-order behavior this constant determines. 



extremely likely that 



5f = ^. (5.85) 



Studying table 2 some more, one also finds that 



^f) = _ls{'\ (5.86) 



S['^ = -Is?\ (5.87) 



5f = (5.88) 

sio)sW = ^s['\ (5.89) 

are satisfied, at least up to the order to which we have calculated the relevant constants. Of 
course, we conjecture these results to be exact, even though we have no clear idea on how to prove 
these relations. Proving these relations and generalizing them to arbitrary Stokes constants^^ is 
a very interesting problem whose solution will very likely give us a much deeper understanding 
of Stokes phenomena and resurgence in the Painleve I framework. 

5.5 The Nonperturbative Free Energy of the (2,3) Model 

As discussed at the beginning of this section, we know that the free energy, F{z), of the (2,3) 
minimal string theory is related to the solution, u{z), of the Painleve I equation by 

u{z) = -F"{z). (5.90) 



^^Using the limited amount of available data, one may make further bold guesses such as n S}" = 
- (n - 1) S'J""^' and nS^^^ = i""^ (Sf ^)^"". Also, it seems natural to write (5.89) as 25*^' = 3S'J^'S'^"\ 
since these two products often occur in the same alien derivatives. It is further tempting to guess that, in general, 
every Stokes constant can be expressed as a rational function of the Sj;"' alone. 
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We now want to investigate how our results for the Painleve I solution translate into results for 
this free energy. Let us begin by studying the perturbative contribution to u{z), 



^(0|0)(^) = E4r™-"^'^"^/'- (5-91) 

9=0 

Here, the reader should note that we have re~inserted the factor of ^/z that we had removed 
ear her in (5.11). Two integrations then lead to 

+ 00 „ „ +00 

i7(0|0)(^) ^ _^40|0)[0] / / d^^-(59-l)/2 ^ ^i7(0|0)^2,-2^ (5 92) 

9=0 -^-^ 9=0 

where x = z~^l^ is the string coupling constant, and where we defined the perturbative expansion 
coefficients for F^^'^^^iz) as 

pirn - ^ ^(o|o)M f5 

- (5,-3)(5<7-5)^^^ • ^^-^^^ 

Our reason for not including a "log index" [0] in the free energy coefficients pf^^^ will become 
clear in a moment. This asymptotic series is, once again, the perturbative part^''' of a transseries 
expansion for the free energy F{z). To see what form the one-instanton contribution takes, let 
us integrate the leading one-instanton term in the u{z) transseries, 

r r (1|0)[0] 

In this expression we have explicitly written the leading coefficient n^^'*^^^*^'. Recall from our 
discussion in section 5.2 that the value of this constant can be absorbed by a rescaling of a\ and, 
for this reason, we have so far worked in a convention where u^"^'*^^^*^' = 1. This was a very useful 
normalization for constructing the two-parameters transseries solution for u(z) but, to discuss 
the free energy F{z), we now actually want to change to a different convention. 

The reason for this new choice of normalization is that we would like our one-instanton 
contribution to the free energy to agree with the equivalent result that was computed in [17, 
15], for the free energy around the one-instanton configuration, straight out of a matrix model 
calculation associated to eigenvalue tunneling. That is, we want our one-instanton contribution 
to have the normalization (compare with, e.g., formula (4.35) in [15]) 



Notice that this coefficient is computed directly from the (2, 3) minimal model spectral curve 
[15]. To find this answer, one simply has to rescale 

3I/4 

cTi = -i;-^^f. (5.96) 



^^Here and in what follows, we will not explicitly include any integration constants. In principle, these lead to 
undetermined terms in F{z) which are constant and linear in z, and which cannot be fixed by using the Painleve 
I analysis alone; they must be derived from the minimal string theory directly. It turns out that naively setting 
these terms to zero actually leads to the correct string theory result. 
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In order to keep the symmetry between instantons and generalized instantons, which we discussed 
before, we shah also choose to rescale (T2 with this exact same factor. 

The reader may have noticed that the above result may be equivalently written as 

ai = 5}°^ a[. (5.97) 

The appearance of the Stokes constant sf^^ in this formula turns out to be quite natural. A 
rescaling of the variables cTj does not only rescale the transseries components, but also the Stokes 
constants. A quick calculation shows that, under a general rescaling of the Uj, these quantities 
scale as follows (recall (4.14) for example) 



0-1 


= cidi, 


(5.98) 


0-2 




(5.99) 


{n\m) 


— C-^ C2 ^(n|m)5 


(5.100) 




_ „l-k S-k-l q{k) 
— C2 0^ , 


(5.101) 


oik) 


— C2 oi . 


(5.102) 



In particular, our rescaling sets = 1. Of course, physical quantities cannot depend on 
arbitrary normalization choices, so any physical quantity must be a scale invariant combination 
of the above quantities. As we shall see in the following, we will be particularly interested in 
quantities which can be made scale invariant by multiplying with powers of the first Stokes 
contstant. When this Stokes constant equals 1, this means that the scale invariant quantity is 
numerically equal to the "bare" quantity. 

Stokes Constants for the Free Energy 

Recall from (5.58) that the large-order behavior of ^4^'*^^^*^^ has a leading term 



^(0|0)[0] ^ 25f V {2g - 1; ^(i|o)[o]_ ^^^^^3^ 
^ 27ri A^^~2 



For Fg^^^^ we can do the exact same large-order calculation and the result is very similar 

(010) ^ 2Sf^ r(^g-j) k(i|0) (5.104) 

where we denoted the leading one-instanton coefficient^^ in the free energy transseries by F^^'^\ 
We see that the only difference between the above two equations is in the argument of the gamma 
function and the power of the instanton action, A. Both of these are shifted by —2, as a result 
of the double integration involved in going from u{z) to F{z). One can see quite easily that this 
is a general property: all large-order formulae for u(z) and F{z) are the same up to these shifts. 

Notice that in (5.104) we have denoted the Stokes constant as s[^^^ . Indeed, nothing guar- 
antees that the Stokes constants for the transseries F{z) equal those for the transseries u{z) — and 



^^For the free energy, and in order to avoid fractional indices, we will use a convention where all perturbative 
series start with a coefficient Fq"^"^\ with lower index 0. 
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in fact we shall see that, in general, they are different. However, the Stokes constants for F{z) 
can easily be obtained from those for u{z). For example, comparing (5.103) and (5.104), we can 
calculate the Stokes constant s[^^^ for the free energy. First of all, note that we can rewrite 
(5.93) as 

pm = i + o -]). (5.105) 

' 25 (2<7-i)(2g-i) V KaJJ ^ ' 

Moreover, we know from (5.94) and (5.97) that 

i.(i|0) = _^^(i|0)[0]. (5.106) 



Inserting both of these in (5.104) gives 



2 



Comparing this to (5.103), we find that 



5j°^^ = 1. (5.108) 



This once again indicates why the rescaling (5.97) was a useful choice. In general, quantities such 
as Fq^^^^ or S^^^ cannot be physically meaningful quantities: only "scale invariant" quantities 
such as 

5f)^.F(^l°) (5.109) 

can carry physical information. Having chosen our present normalization in such a way that 
^(o)F _ have shifted the full physical information into Fq^^^\ We could of course just as 

well have done the opposite thing, i.e., choosing a normalization where Fq^'*^^ = 1 and absorbing 
all physical information into sf*^^ . The reason we have chosen the present normalization is 
that it agrees with the one usually chosen in the litarature. For example, in [15] the above 
normalization is chosen and the resulting physical quantity i^Q^'*^^ (called /ii in that paper) is 
calculated directly from the spectral curve of a matrix model. 

The above calculation, relating 5^*^^^ to can be repeated for any Stokes constant. One 
simply finds a term in a large-order formula in which the Stokes constant appears, calculates the 
normalization of this term for both u{z) and F{z), and then compares the two formulae. Doing 
this carefully one finds the following relations between the Stokes constants for F{z) and for u{z) 

2k+e-2 



Sf^ = ^2^f " , (5.110) 

5f)^ = ^2 (5.111) 



Note that the right-hand side in these equations consists of scale invariant quantities; the left- 
hand side consists of quantities which are implicitly scale invariant as well, due to the analogous 
powers of S[°^^ = 1. The factor of comes from taking a second derivative of the instanton 
factor ex.p{±£Az^^^) in the free energy transseries. 

We have listed the numerical values for the free energy Stokes constants in table 3. Our main 
reason for listing the free energy Stokes constants separately is that we expect those numbers 
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1.000000000000000000... 
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2.000000000000000000... i 


20 




-3.000000000000000000. . . 


13 




-1.811460182210655615... 


81 




-5.434380546631966844... i 


19 




1.168020496900498115... 


36 




0.905730091105327807... 


112 




-0.778680331266998743... 


108 


gWF 


0.319744372344502079... 


108 



Table 3: Stokes constants for the (2, 3) minimal string free energy. The third column gives the number 
of decimal places to which the answer was computed. See table 2 for the corresponding quantities for the 
Painleve I solution u{z), from which the above numbers are derived using (5.111). 



to be the ones that can eventually be calculated from minimal string theory or spectral curve 
considerations, similar to the way in which one can calculate Of course, to actually carry 

out such calculations, one needs a physical understanding of what the generalized instantons are. 

As was the case for u{z), not all of the free energy Stokes constants are independent: using 
equation (5.111), the relations (5.86-5.89) directly translate into relations between these numbers 

SP^.I^P, (5.112) 



5P = _^5P, (5.113) 
:WF _ o- cii)F 



5^^^" = 3iS\''\ (5.114) 

As in footnote 36, one can then conjecture analogous further relations for the free energy Stokes 
constants that have not been calculated yet, such as, e.g., Sn^^ = i"~^n. 

The Free Energy Transseries Coefficients 

We now want to calculate the explicit form of some of the (n|m)-instantons contributions to the 
free energy transseries. In the case where n = m, there are no logarithmic contributions to the 
u-transseries, and the double integration is easily carried out as we did for n = m = in (5.93). 
Let us therefore study the "off-diagonal" (n|m)-instantons contribution to u[z) 

+ 00 

^ ^ ("1"*)^ ^12 23+2/3^1 ' ^ ' 

g=0 

where n ^ m. It is convenient to add to this term all the logarithmic terms that are proportional 
to it by (5.40), i.e., all terms of the form 

e-(«-)^/-^ log^u;) • $!tVfc|,n+fc)(^)' (5-116) 
with A; > 0. Using (5.40), we can rewrite these terms as 

i (-i(m-n)cTia2 log u;)' ^a^"^ e-("— )^/-' $[°j|^)(u;), (5.117) 
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and summing all of them over k we find that we can incorporate all of those terms by simply 
replacing ^f^^^^ by 

Formally, we can write this as^^ 

[sum] / N 45(™-")'710"2 ai[0] / N /r- 1in\ 

Rewriting the result in terms of z = w~^^^ and reintroducing the scale factor z^/^, we get the 
following (n|m)-contribution to the free energy transseries 

+ 00 [0] 

,n-.m -(n-m)A^5/4 (n|m)[0] lOg+sp -4 ^ 4(n-^).,.2 



a^a^ e 



Y^^in\mm ^- + (5.120) 



9=0 

To integrate this part of the transseries, we use the fact that 

dz z-^^-'^^'^' = ^ g-M..V4 ^ ^^^^^ _ (^_Mz5/4) , (5.121) 

fc=l 



where 

afc(7) = —^^ ^ ^ ^ (5.122) 

is a polynomial of degree /c — 1 in 7. It is important to notice that, in the components of our 
logarithmically summed transseries (5.120), the coefficient 7 is linear in aia2 and thus a/c(7) in 
(5.121) above will be a polynomial of degree A; — 1 in cji(T2. This means that integrating the 
(n|m) transseries component in u[z) will not only contribute to the {n\m) transseries component 
in F{z), but also to all (n + a|m + a) components with a > 0. 

Using (5.121), the double integration of the u-transseries is now easily carried out in a 
computer. We find the result that the free energy has the following transseries structure 

+00 +00 

F{Z, af , ) =Y.Y. {sf^Y {^^F^^ ^-in-m)A.^'^ ^i(n-mK.r fH-)(^), (5.123) 

n=0 m=0 

where the i^("l™)(2;) are perturbative expansions'*^ in the string coupling z~^^^. The formal power 
of z should once again be interpreted as 

^i(n-mK<x2^ = exp (n - m) af logz^ , (5.124) 



^^This also illustrates in a rather clear way, and as explained before, that the logarithmic sectors do not seem 
to represent any new nonperturbative sectors. Herein, they simply encode an irrational power function. 

''"We use this term with a bit of hand-waving since these expansions contain half-integral overall powers of the 
string coupling constant and two logarithmic terms actually appear in the lowest F'^"^'^\z). 
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which can be expanded to give log z-dependent contributions exactly analogous to the ones we 
found for the n-transseries. That is, we could leave out this factor in (5.123) and instead replace 



with 



+00 +00 +00 +00 +00 

n=Om=0 n=0m=0fc=0 



(n|m) 



1 /I 



(5.125) 

(5.126) 
(5.127) 



, , , (n — m) log z I F; , n (z) 

Keeping the (5.123) transseries structure for the free energy, the first few of the f("I"^)(z) are 
4 J 



i?(3 
i?(4 



I0)(z 
|o)(z 
I0)(z 
l^)(z 
l^)(z 
l^)(z 
I2)(z 
I2)(z 

l^)(. 



15 

1 

12 
1 

'288 



-22 
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16957 _i5 

=z 4 
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-5 
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+ 
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-z 4 
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21233664\/3 
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(5.128) 
(5.129) 

(5.130) 

(5.131) 
(5.132) 

(5.133) 
(5.134) 
(5.135) 
(5.136) 



One easily checks that inserting these expansions in (5.123), and taking minus its second deriva- 
tive, reproduces the results for the u{z) transseries that we listed in appendix A. We only listed 
the F('^l™)(z) with n > m here; the ones with n < m can be obtained by the rule 

pim\n) ^ ^_-^y+ln/2],pin\m)^ n > m. (5.137) 

The starting exponent of ^("1"*) follows straightforwardly from that of u^"'™^. One has 

~ z-t^'-+i, (5.138) 
i?!"!™) ^ ^-Ia'I^ (5.139) 



where (3nm is defined in (5.34) and the second line above is valid for n ^ m. This concludes the 
nonperturbative solution to the (2, 3) minimal string. 
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6. Matrix Models with Polynomial Potentials 

While there are many examples of exactly solvable matrix models (see, e.g., a few such examples 
within the context of nonperturbative completions in [24]), it is certainly the case that in most 
situations one does not have access to anything more than perturbative techniques, most noto- 
riously those introduced a long time ago [6, 25, 7]. Enlarging these old techniques by the use of 
resurgent analysis naturally becomes of critical importance for the extraction of nonperturbative 
information out of a rather large class of string theoretic examples [12]. In here, we shall focus 
upon matrix models with polynomial potentials, mostly on the quartic one-matrix model, devel- 
oping the two-parameters resurgent framework as it applies to this example. Notice that in the 
large N limit all matrix model quantities will now depend upon 't Hooft moduli, an additional 
complication as compared to the case of minimal strings. However, we shall further see how to 
make the bridge back to Painleve I via a natural double-scaling limit of the quartic model. 

The resurgent analysis of matrix models has another added feature, as compared to minimal 
string models. Within this context, perturbative techniques construct asymptotic expansions 
which are formal power series around a given saddle-point of the partition function of the the- 
ory. In other words, one performs perturbation theory around a chosen background — where one 
expects that a full nonperturbative solution should be background independent, i.e., it should 
include all possible backgrounds [26]. This is where the full transseries framework comes into 
play: only by properly considering the correct multiple-parameters transseries (a too-parameters 
transseries in the quartic example) can we expect to construct fully nonperturbative, background 
independent solutions. In fact, while it is possible to consider a one-parameter transseries ansatz 
for the quartic matrix model, still yielding a rather interesting amount of nonperturbative in- 
formation, this is not the most general multi-instanton expansion required and, as such, cannot 
possibly see all other backgrounds [12]. In the following we shall construct the full two-parameters 
transseries solution to the quartic matrix model around the so-called one-cut large saddle- 
point. Because this is the most general solution to this problem, it is naturally applicable to the 
problem of changing of background: one can envisage starting off in the one-cut phase and, via 
Stokes transitions, reach other stable saddle-points of the quartic matrix models such as, e.g., 
its well-known two-cut phase. We hope to report on these issues in upcoming work. 

6.1 Matrix Models: Spectral Geometry and Orthogonal Polynomials 

For the purpose of completeness on what follows, let us begin with a lightening review of matrix 
models, both in the spectral geometry and orthogonal polynomial frameworks (for more complete 
accounts we refer the reader to, e.g., the excellent reviews [5, 61]). 

The one-matrix model partition function for the hermitian ensemble is 



with 't Hooft coupling t = Ngs (fixed in the 't Hooft limit). In standard diagonal gauge one has 



where A(Ai) is the Vandermonde determinant. The simplest possible saddle point for this integral 
is the one-cut solution, characterized by an eigenvalue density normalized to one, and where 
the cut is simply C = [a,b]. A rather convenient description of this saddle point is given by 
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the Riemann surface which corresponds to a double-sheet covering of the complex plane with 
precisely the above cut. This geometry is described by the corresponding spectral curve^^ 



y{z) = M{z) ^/{z-a){z-h), (6.3) 

where^^ 

X / du; V'il/w) 1 , 

M{z)=i —= =. 6.4 

7(0) 2m 1 - wz ^{l-aw){l-bw) 

For future reference, it is also useful to define the holomorphic effective potential ^h-eff(-^) ~ vi^)^ 
which appears at leading order in the large N expansion of the matrix integral as 



i=l \ i=l / 



(6.5) 



There are many ways to solve matrix models. A recursive method, sometimes denoted by the 
topological recursion, was recently introduced for computing connected correlation functions and 
genus g free energies, entirely in terms of the spectral curve [8, 9]. However, for our purposes of 
computing the genus expansion of the free energy, one of the most efficient and simple methods 
is still that of orthogonal polynomials [25], which we now briefly introduce. Considering again 
the one-matrix model partition function in diagonal gauge (6.2) it is natural to regard 

dfi{z) = e-^^(^) — (6.6) 
27r 

as a positive-definite measure on M, and it is immediate to introduce orthogonal polynomials, 
{pn{z)}, with respect to this measure as 

dfj.{z) Pn{z)pm{z) = hn6nm, n > 0, (6.7) 

where one further normalizes Pn{z) such that Pn{z) = z" + ■ ■ ■ . Further noticing that the 
Vandermonde determinant is A(Ai) = detpj_i(Aj), the one-matrix model partition function 
above may be computed as 

N-l N 



Z=llhn = h^ll r^-, (6.8) 

n=0 n=l 

where we have defined = for n > 1, and where one may exphcitly write 

ho = / d/x(z) = — / dze~^^^'\ (6.9) 



The r„ coefficients also appear in the recursion relations of the orthogonal polynomials, 

Pn+l{z) = {z + Sn)Pniz) - rnPn-l{z), (6.10) 

together with the new coefficients {sn}, which actually vanish for an even potential. 

Where the imaginary part of the spectral curve simply relates to the eigenvalue density. 
''■^This particular expression only holds for polynomial potentials. 
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The key point that fohows is that once one has a precise form of the coefficients in the 
recursion (6.10), one may then simply compute the partition function of the matrix model (and, 
in fact, all quantities in a large N topological expansion). In the example of main interest to us 
in the following, that of the quartic potential V{z) = }^z^ — ^ z^, it is simple to find that Sn = 
and [25] 



fn ^1 - {rn-i +rn + ?^n+i)^ = ngg. (6.11) 



This recursion sets up a perturbative expansion around the one-cut solution of the quartic matrix 
model which, as briefly outlined above, is described by a single cut C = [—2a, 2a] where 



a' 



1 



1 - Vl - 2Xt) (6.12) 



and the spectral curve is 



+ 2a^)^ Vz^-Aa"^. (6.13) 



Before attempting a nonperturbative transseries solution to the quartic matrix model, let us 
briefly consider its perturbative solution [25] and what it implies towards resurgence. 

6.2 Resurgence of the Euler MacLaurin Formula 

In the 't Hooft limit, where — >• +oo with t = QgN held flxed, the perturbative, large N, 
topological expansion of the free energy F = logZ of the matrix model (6.1) is precisely given 
by a standard string theoretic genus expansion (1.1). This is usually normalized against the 
Gaussian weight, where Vg{z) = \z^, thus following from (6.8) 

g=0 ^ U ^■^ n=l " 

In this expression one flrst needs to understand the large expansion of the recursion coefficients, 
{r„}. Given the Gaussian solution = ngg it is natural to change variables^^ as x = ngg, where 
X ^ [0, t] in the 't Hooft limit, and define the function 

7^(x) = r„, with 7^^(a;) = x. (6.15) 

In the example of the quartic potential, (6.11) is then rewritten as 

n{x)i^l-^{n{x-gs)+n{x)+n{x + gs))^=x. (6.16) 

Noticing that this equation is invariant under gg ^ —gs it follows that TZ[x) is an even function 
of the string coupling and thus admits an asymptotic large N expansion of the form 

n{x) = Y,gl'R2g{x), (6.17) 

9=0 



The X variable in tiiis section should not be confused with the x variable of section 5. 
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which ahows one to solve for the R2g{x) in a recursive fashion, given i?o(0) = 0. Further noticing 
that in the 't Hooft hmit, where x becomes a continuous variable, the sum in (6.14) may be 
computed by making use of the Euler-MacLaurin formula^^ 



lim — 



N 



n=l 

we finally obtain 







1=1 +00 



1 B' 



(6.18) 



log 



n{x) 



rIS \ 

+00 



(2^)! dx'^9-1 



+ 

x=o/ 9. 
{t - x) lo, 



1 /■* 

/ dxit- x) log 

n{x 



n{x) 



+ 



x=t 



x=Q 



(6.19) 



or, explicitly, using the expansion of TZ{x) in powers of the string coupling [25], e.g., 



jro(i) = f dx{t-x) 
Jo 

T^{t) = I dx{t-x) 
Jo 



log 



Roix) 



R2{x) 

Ro{x) 



1 d 
12 dx 



(t — x) \oi 



Roix) 



x=t 



+ -t\. 



(6.20) 
(6.21) 



x=0 



It is worth making some comments concerning these expressions. First notice that the Euler- 
MacLaurin formula is an asymptotic expansion, thus only capturing perturbative contributions 
to the matrix model free energy. These perturbative contributions to the free energy at genus 
g then arise from a recursive solution to the string equation, i.e., out of the coefficients R2g{x), 
computed recursively in the quartic potential example out of the large string equation (6.16), 
and similarly for different potentials. For instance, in our main example, it is simple to obtain 
out of (6.16) that 

i?o(x) = ^ (l - Vl - 2Ax) , (6.22) 

which is the one solution satisfying the initial condition Ro{0) = 0. Finally, when computing 
J^g{t) above, we have partially restricted our result to the quartic potential, as we have made use 
of the fact that in this case one has [25] 



ho n{x) 



log 



x=0 



ho{+\\\) 
ho{-\\\) 



1 



11 

48 



>?9l 



(6.23) 



We now arrive at the main point concerning the construction of perturbative solutions to 
matrix models, in the orthogonal polynomial framework, and its relation to resurgence. It can 
be shown that the asymptotic expansion (6.18), defining the Euler-MacLaurin formula, may also 
be written as a finite difference operator of Toda type [12], 



F{t + gs)- 2F{t) +F{t-gs) = log 



n{t) 



(6.24) 



*In here the B2k are the BernouUi numbers and x = t^. 
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This expression encodes the relation between TZ{t,gs) and J^{t,gs), expressed by the Euler- 
MacLaurin asymptotic formula, and it is essentially the large N version of the identity 



which is in itself an immediate consequence of (6.8). The above finite-difference equation makes 
it clear that if the recursion function TZ{t, gs) has a non-trivial resurgent structure, arising via 
a transseries solution to the string equation (6.16), then, (6.24) will immediately induce a non- 
trivial resurgent structure to the matrix model free energy J-'{t, gs), of the exact same form [12]. 
This is essentially a statement concerning the particular solution to the non-homogeneous Toda- 
type relation above. One has, however, to check the general solution to the homogeneous version 
of (6.24), i.e., check whether the Euler-MacLaurin formula induces any other new resurgent ef- 
fects before further ado! But all such homogeneous "Toda" resurgent effects have already been 
studied in [24]. Furthermore, it was shown in [62] that, essentially because the homogeneous 
Euler-MacLaurin relation (6.24) is linear with constant coefficients, it only has Borel singular- 
ities associated to A-cycle instantons, of the type discussed in [24]. In this scenario, B-cycle 
instantons, displaying fully non-trivial resurgence, originate in transseries solutions to the string 
equation (6.16). These translate to the free energy as the non-homogeneous contribution to the 
solution of (6.24) (relating back to our discussion on A and B-cycle instantons in section 3). 

The bottom line is thus that the nonperturbative resurgent analysis can be all done at the 
level of the non-linear recursion, or string equation (6.16), alone. This will capture the full non- 
trivial resurgent structure of the matrix model free energy; the addition of "Toda" or A-cycle 
instantons then being completely straightforward to implement, following the results in [24]. 

6.3 The Transseries Structure of the Quartic Matrix Model 

As just discussed, the solution TZ{x) to the string equation (6.16) completely determines the free 
energy of the one-cut solution to the quartic matrix model. In order to nonperturbatively solve 
this model, our aim is now to construct TZ{x) as a transseries solution. 

The string equation is, in this case, the finite-difference analogue of a second-order differ- 
ential equation. For this reason, one expects the full transseries solution to contain two free 
parameters, which is further consistent with the fact that the double-scaling limit of the quartic 
matrix model reproduces the (2, 3) minimal string theory. As we have seen, the free energy of 
that theory is described by the Painleve I equation, which is also solved by a transseries with 
two free parameters. The one-parameter transseries solution to the string equation (6.16) was 
first discussed in [12], building upon the perturbative results obtained in [25]. Below, we review 
those results, and then continue to describe the full two-parameters transseries solution. 

Review of the One Parameter Transseries Solution 

In [12] the one-parameter transseries solution to the string equation 



was investigated. It was found that such a solution can indeed be constructed, having the form 



Zn+iZn-1 



(6.25) 





(6.26) 




(6.27) 



+00 

= e-'^M^)l9s ^g9ji(n)^x), (6.28) 

5=0 

where we used a notation which shghtly differs from that in [12] but which is more convenient for 
our purposes. Note that, as in the Painleve I case, the nonperturbative answer is an expansion 
in the "open string coupUng constant", gs, and not in the "closed string couphng constant", g^. 

To find expressions for the R^g^\x), one simply plugs (6.28) into (6.26), and solves the 
resulting equation order by order in a and gs- For example, at order a^g^, one finds the equation 

l--]=x, (6.29) 



where we have introduced the shorthand 



do\x). (6.30) 



Solving this quadratic equation leads to the answer we have already mentioned, 

r = j(l- Vl - 2Aa;) . (6.31) 

Here the square root is defined to be positive on real and positive arguments and we chose the 
sign in front of it in such a way that r has a finite A — t- limit. 
At order a^g^, (6.26) gives the equation 

R^'^ {x){l- Ar) - ^ = 0. (6.32) 

Using (6.31), one can now solve for R^\x) in terms of x. In fact, it will turn out to be useful 
to write this answer, as well as all other answers that will follow, in terms of r. Doing this, one 
obtains 

This procedure is easily continued to order cr^gs^ , which then determines all coefficients R^g{x). 
In this way, one reproduces the perturbative results that were first obtained in [25]. Note that, 
at order a^, we are skipping all odd orders in gs since our answer should be an expansion in the 
closed string coupling constant g^. As aforementioned, since equation (6.26) is itself even in gg, 
it is indeed possible to find a perturbative solution 7^pcrt(a;) which is also even in gs- 

The next step is to calculate the one~instanton contributions, which appear at order a^. 
Expanding (6.26) at order cr^gg, one finds 

?Wf^^ ( ^+A'{x) ,-A'(x) ■ 6 



R);> {x) (^e+^ + e"^ + 4 - ^ J = 0. (6.34) 

One sees that the overall factor R^\x), which we will soon find to be nonzero, drops out. Hence, 
this equation determines the possible values for the instanton action A[x). Expressed in terms 
of the variable r, these values are 

r / 3 \ 1 

A{x) = ±- (2 - Ar) arccosh 2 T ^3(1 - Ar) (3 - Ar) + 

z y Ar J ZX 
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+^ipr(2-Ar) + Cint, (6.35) 

where the branch cuts are chosen such that both the arccosh and the square root are positive 
when Ar — )• 1~, and p € Z. Furthermore, notice that in the first hne there is only one single sign 
ambiguity: one can either choose both upper signs or both lower ones. The integration constant 
Cint and the integer ambiguity p were fixed in [12] by requiring that this expression reproduces 
the Painleve I instanton action in the corresponding double-scaling limit. It turns out that, for 
this, both constants need to vanish. The sign in the first line was also fixed in [12]; to obtain the 
positive Painleve I instanton action, one needs to choose 

r /3 — 2Ar\ 1 

A{x) = -- (2 - Ar) arccosh — + — ^(3 - 3Ar) (3- Ar). (6.36) 

2 y Ar J 2A 

In our two-parameters case, we shall eventually be interested in both choices of sign. We simply 
take the above expression as the definition of A{x) and, once we move on to the two-parameters 
transseries, one will see that both A(x) and —A(x) appear symmetrically. 

Akin to the Painleve I case, essentially the same results may be obtained by writing the 
string equation (6.26) in prepared form. Indeed, also for finite difference equations there is a 
very similar story to the one we described in section 4 for ordinary differential equations, and 
which we shall now mention very briefiy [63]. This time around one can show that, via a suitable 
change of variables, a rank-n system of non-linear finite difference equations 

n{x + l) = F{x,n{x)), (6.37) 

may always be written in prepared form as [63] 

nix + 1) = A{x)n{x) + G{x,n{x)), (6.38) 

with G{x,'R{x)) = O {\\7l\\^ ,x^'^7i}j and where 

A(x) = diag (e-'^i (l + x-i)^\...,e-"" (l + x"!)''") . (6.39) 

Within this setting formal transseries solutions to our system of non-linear finite difference 
equations essentially have the same form and properties as those discussed in section 4. 

Once we have fixed the instanton action A (to keep the notation readable, we shall many 
times suppress the x-dependence of all our functions), one can continue to higer orders in Qs- At 
order a^g\ (6.26) gives terms involving two unknown functions, E^'^ and R^i \ However, it turns 
out that the terms proportional to R^^^ actually are 



i2i^^{e+^ + e-^ + 4-^), (6.40) 



and hence vanish by (6.34). One is left with the equation 



dx \ 7^2 

This differential equation is not too hard to solve; where the multiplicative integration constant 
is once again fixed by requiring that the double-scaling limit yields the Painleve I solution [12]. 
Its solution is thus 

= TlJ^ 777 > (6-42) 

° (3-Ar)i/^(3-3Ar)i/4' ^ ' 
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where the quartic roots are defined to be positive as Ar — )■ 1 .In fact, we shall use this convention 
for any of the fractional powers that will appear in what follows. 

Proceeding in this way, one finds a similar pattern: at order a^gs both R^gli and R^^'^ appear 

as unknown functions, but rI^^ multiplies the same terms as R^^ in (6.34) and hence drops out. 
This is nothing but the phenomenon of resonance that we have also encountered in the Painleve 
I case. What is left is a linear first-order differential equation for R^gli, which can then be easily 
solved. The integration constant in this solution can be fixed by the requirement of a good 
double-scaling limit. In [12], the answers for i?f ^ and R^^^ were calculated in this way. Using a 
Mathematica script, we have calculated the one-instanton contributions up to -R3o^. The general 
structure of these solutions will be described below. 

In principle, one could now go on in the same way and calculate the higher instanton con- 

(n) 

tributions Rg , for n > 1. Instead of doing this in the one-parameter formalism, we shall now 
move on to the two-parameters case, and calculate the higher instanton contributions as part of 
this more general setting. 

The Two Parameters Transseries Solution 

In the framework of the present paper, one should not restrict to a single sign choice for the 
instanton action. Rather, we would like to find the general two-parameters transseries solution 

+00 +00 

^(^) = E E <^2"^^^"''"^(^), (6.43) 

n=0 m=0 

+00 

R(n\m)^^^ ^ ^-in-m)A{x)/g, ^ 5fi?H'")(x), (6.44) 

Q—i^nm 

to the quartic model string equation. Note that, apart from the obvious changes in this ansatz, as 
going from one parameter a to two parameters ai, a2, we have also included a "starting genus" 
/3nm) which plays the same role as the Pnm in our previous examples. The reader may also wonder 
if it is not necessary, as in the Painleve I case, to introduce log gg terms in our ansatz. As we 
shall see below, there is in fact no need for such terms in the present context^^. 

Once we have made this ansatz, solving the string equation (6.26) order by order in n, m and 
g \s a, tedious but relatively straightforward exercise. As in the one-parameter case, one simply 
inserts (6.44) into the string equation, isolates the terms multiplying a certain power of ci, (T2 
and gs, and solves the resulting equations inductively for i?g"'™^(x). 

In the case of the ordinary instanton series one finds algebraic equations for R^^^'^\x), with 
n > 1. For example, at order o"f(T2fi'^, one finds the equation 



, , , / a\ p(l|0) p(l|0) 



which, after inserting (6.42) and (6.36), is solved by 



4"°^ = # 372- (6-46) 

° 2(3-Ar)i/2(3_3^^)3/2 



''^As we will see in section 6.4, however, it may be useful to change variables in such a way that logpa terms do 
appear. This will turn out to be especially useful when we want to study the double-scaling limit. 
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Going beyond the instanton series, we can now also calculate the "generalized instanton con- 
tributions" R^J^^"^\x), with nonzero m. At order (t|(T2(?s, for example, one finds an algebraic 
equation that is solved by 

(3 - Ar)^/2 (3 - 3Ar)^/2 

Continuing to higher genus, one finds that all -Rg^'^^ with odd g vanish, so that the resulting 
perturbative series is a series in the closed string coupling constant g'^. The same holds for all 
other functions Rg^^"'\ with as many instantons as "generalized anti-instantons" . 

At generic order aia^gi one has to solve an algebraic equation to find i?^"''"^ Generically, 
i.e., when n ^ the answers also contain "open string" contributions with g odd. When 
n = mil, we again encounter the phenomenon of resonance: the terms multiplying i^^"'"*^ 

drop out, and we actually need to solve a differential equation to obtain i?^'^™^ Some of the 
integration constants that appear in the solutions to these differential equations are equivalent 
to the ambiguities we found in the Painleve I case: they parameterize the choices we have in 
rearranging ai and o"2 into new nonperturbative parameters. We fix those integration constants 
as for Painleve I, by requiring that jSnm is as large as possible. Other integration constants do 
not have this interpretation, and need to be fixed by requiring the correct double-scaling limit. 

The solutions to the differential equations for n = m ± 1 are not all of the form that we have 
encountered so far. Starting at n = 2, m = 1, we also have logarithms entering the game. For 
example, for i?Q^'^\ we find that 

(2|i) _ A\/A^ (54 - 45Ar - GA^r^ + SAV^) 3X^/>^ (6 + 3Ar - 6X^r^ + 2X^r^) 
^° 4r (3 - 3Ar)iV4 (3 _ 32r (3 - 3Ar)iV^ (3 - Arf 



with 



(6.48) 
(0.49) 



Note that, once again, we see logarithms appearing as was previously the case for the Painleve 
I equation. The big difference as compared to the aforementioned situation is that now the 
logarithmic factors are functions of x, and not of the perturbative parameter gg. As it turns out, 
all instanton corrections still take the form of open string theory perturbation series. Only in 
the double-scaling limit (where, as we shall see shortly, x becomes a function of gg) do we find 
back the logarithmic coupling constant dependence of the Painleve I solution. 

Another interesting result is that, generically, the "starting genus" /3„m in (6.44) is nonzero. 
In fact, it is usually negative: for example, one finds that the series for n = 2, m = 1, does not 
start with the above function but with 

SO that /32,i = — 1- We find that the non-logarithmic terms have a true genus expansion in gs, 
but that the expansion for the logarithmic terms actually starts at "genus —1/2". At higher 
generalized instanton numbers, the non-logarithmic terms will in general also appear with neg- 
ative powers of gs. While this may seem surprising, it is not a big problem: as we shall see in 
section 6.5 the transseries solution for the free energy of the quartic matrix model still only has 
nonnegative genus contributions. 
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Two— Parameters Transseries: Results 

We have written a Mathematica script to solve the equations for i?g"''"^(x) to high orders in n, 
m and g. In appendix B, we present some further exphcit results. Here, let us write down a 
formula for the generic structure of the answer: 

W (3 _ 3Ar)P3 (3 _ Ar^ ^3 ^^■^^> 

where the powers in the prefactor are the following functions of n, m and g, 

Pi = ^ (3n - m - 2) , (6.52) 
P2 = n + m + g — 1, (6.53) 
P3 = ^ (5n + 5m + Wg - 4) , (6.54) 

P4 = ^ (3n + 3m + 6g + 26-4), (6.55) 

with 5 = (n+m) mod 2. In general, the gs expansion starts at g = (5nm = — min(n, m), whereas n 
and m only take on nonnegative values. Finally, at each order (n, m, g) we find a finite expansion 
in logarithms, 

min(n,m.) 

pH'»)(;r)= ^ PHHW(2;) .iog'=/(x), (6.56) 

A:=0 

with f{x) the function defined in (6.49). The resulting components pj"''"^''^' are now polynomials 
in Xr, of degree {6g + n + 5m + 6 — 2) /2. 

These formulae look somewhat complicated, but the crucial point is that all the information 
about the two-parameters transseries is now contained in a set of simple polynomials. Moreover, 
up to an overall rational factor consisting of powers of some small prime factors, the coefficients 
of these polynomials are integers. Thus, we have reduced the full nonperturbative solution of the 
quartic matrix model to the determination of a list of (n + m + 63 — 5 + 2)/2 integers for every 
n, m, k and g. This result makes one wonder if these integers have any further relations between 
them, and whether they contain any geometrical information, as for example in the case for GV 
invariants we have discussed in section 3. We have no concrete suggestions in this direction, but 
it would be very interesting if such an interpretation could be found. 

The reader may have observed that both the power pi and the degree of the polynomials 
are not symmetric under the exchange of n and m. The reason for this is that we wrote (6.51) 
in such a way that, in general, when n > m, the pj"'!™''^'^! are irreducible polynomials^^. When 
n < m the structure formula is still valid but the polynomials are no longer irreducible. In fact, 
the symmetry of the string equation dictates that 



g ^ ' g 

and as a result there is a relation 



pHm) = (_i)9pHn)^ (6.57) 
p{n\n,)[k] ^ (Ar)2™-2" (6.58) 



''^There are a few low-index exceptions to this rule, for example, p^'^^M irac f^r.i^r^^ \^ p(M''-)l' 

and all p^^_^l^^l*l contain a factor of (2 — \r). 
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when n < m. After inserting this back in (6.51), the symmetry in n and m is indeed restored. 

When n = m, the polyno: 
these polynomials factorize as 



When n = m, the polynomials pj*^'"^^^"' are highly reducible. It turns out that in this case 



p(n\nm ^ (^^)P2-Pi Q(n)^ (g ^g) 



with qIP^ a polynomial of degree 2g + 2n — 1 in Ar. Thus, in these cases, one can rewrite (6.51) 
as follows 

^»"""'->= (3-3Arr(3-A,0" Q^"''"'- 

(n) 

When n = 0, Qg factorizes even further, and can be written as 

Q(^) = Xr{3-Xrf^Sg, (6.61) 

with Sg a polynomial of degree {g — 2)/2 in Ar (recall that Sg is only nonzero for g even, so that 
this degree is always an integer). Thus, we now have 

This expression is only truly valid for > 0, although formally we can use it for = as well if 
we choose 

as the "degree —1 polynomial". 

6.4 Resurgence of Instantons in Matrix Models and String Theory 

Now that we know the full structure of the one-cut two-parameters transseries solution to the 
quartic matrix model, we can test the theory of resurgence as described earlier in this paper. 
Before we do this for the full solution, let us discuss the double-scaling limit, in which the string 
equation reduces to the Painleve I equation that we studied in section 5. 

Double Scaling Limit 

It is well-known that there is a double-scaling limit in which the double-line Feynman diagrams 
of the quartic matrix model reproduce the worldsheets of the (2, 3) minimal string theory (for 
details on the physical aspects of this relation, we refer the reader to the review [5]). At the level 
of equations, it is not too hard to see directly that this limit exists. To this end, we first change 
variables from {x,gs) to 



(- 



■ 1 — 2Xx \ 

{z,gs) = 2, 9s , (6.64) 

and replace TZ{x,gs) by a function u{z,gs) using the substitution 

n{x, gs) = ^ (l - (8A2<72) ^ u{z, gs)) . (6.65) 

A little algebra then shows that, in the limit where gs ^ and z is held fixed, the string equation 
(6.26) indeed reduces to the Painleve 1 equation 

u^{z) - ^u"{z) = z. (6.66) 
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Note that this result is true for any value of A. This is a consequence of the fact that we 
have a redundancy of variables: the coupling constant A in the quartic matrix model potential 
V{M) = \M'^ — can essentially be absorbed into gs (or the 't Hooft coupling t = ggN) by 

a rescaling of M. We will encounter this redundancy of variables a few times in what follows. 

Of course, the fact that the string equation reduces to the Painleve I equation does not 
automatically imply that the same is true for the particular solutions TZ{x,gs) and u{z) that 
we have constructed. It is well known (see, e.g., [5]) that this is nevertheless the case at the 
perturbative level; 

TZpcrt{x, gs) 'Upcrt(z) (6.67) 

in the double-scaling limit. One might therefore hope that the same holds true for the full 
two-parameters transseries solutions. It turns out that this is indeed the case, but not in a 
completely straightforward way. As we shall see, the correct double-scaling limit also requires a 
subtle transformation between the nonperturbative ambiguities (iTi,o"2) for the two solutions. 

To further understand this limit, let us look at the full two-parameters transseries solution 
TZ{x). It turns out to be useful to make some shifts in the summation indices, and write the 
transseries in the form^^ 

+00 +00 +00 +00 

n=0 m=0 g=l3nm k=0 

Here, we have used the shorthand (6.49) 

(3-Ar)^(3-3Ar)^ 
/(^) = ^^^4^:4 ' (6-69) 

and split the i?^'^'™'' components into logarithmic contributions in the obvious way 

min(n,m) 

^(n|m)(^^^ E log'(/(^))-4"''"^f'k^)- (6.70) 
fc=0 

The reason for writing TZ{x) in the above form is that we may now apply the same trick as we 
did for the Painleve I solution: from (B.26) and (6.51) one easily deduces that 

{n+k\m+k)[k] _ 1 f X{n-m) Y (n|^)[0] . . 

- J., { u J ' (6.aj 

so that we can sum the full logarithmic sector in order to find 

+00 +00 +00 

^(^) = JZ5Z Yl ^^rf^re-^"-'")'^^"^/^^5f4"''"^'°k^)-(/(a^))^^"""^"'"'- (6.72) 

n=0 m=0 g=/3„r„ 

Next, we want to manipulate this expression in such a way that it gives the correct double-scaling 
limit, u{z). To this end, we note that, in this double-scaling limit^*^, one finds 



''^Recall our convention that Rg""^"^^ = if g < /3„ 



^^In here we have scaled i^'"!™'!"' with the same overaU factor of —A ^ {SX^gl) ^ that was present for R 



(0|0)[0] 

a 
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/ — > ^IMX^glz^'^, (6.74) 

- ^QMM ^ ^Pi (6.75) 
9s 

where we denoted the quartic matrix model instanton action by j4qmm and the Painleve I in- 
stanton action by A-pi. In what follows, whenever there is danger of confusion, we will label 
Painleve I quantities with a subscript PI and the analogous quartic matrix model quantities with 
a subscript QMM. When no subscript is present, we always refer to the quartic matrix model 
quantitiy. In the above formulae, we have also introduced the constant 

C = -^^. (6.76) 

Of these three double-scaling formulae, the last two can be simply derived from their definitions. 
However, since we have no closed form expression for ijj"'™^'"', we cannot derive the first — we 
shall see nonetheless that it is necessary for the double-scaling limit to work. Moreover, we have 
explicitly checked its validity on all of the (more than 100) ^"''"^["l that we have calculated. 

As an example, consider the expressions for Pg'^^^^ (the polynomial components of i?^^''^^'''^) 
in (B.12-B.15). In the double-scaling limit we find that they yield the following terms: 

1^-3/4 _ 55 ^-2 _^ 1325 ^-13/4 _ 3363653 ^,9/2 _^ 

6 576^/3 36864 53084160\/3'' 

After removing the overall normalization of y/z and substituting z = w~^^^, this reproduces 
the ti(2:)-component *I'(2|o) {■^•^)^ should be expected. More generally, for the polynomial 

components pj"'™"^ of the fi^"'™^ coefficients that we present in appendix B, one easily derives 
from (6.73) that the double-scaling limit gives^^ 

$(nM(^) = - E (-^) 23.2V2 35./2 ^^^'"^^(1)- (6-78) 

We thus conclude that the double-scaling limit works nicely at the component level. However, 
when inserting (6.73-6.75) into (6.72), we see that for the full TZ{x) the naive double-scaling 
limit has two problems: 

1. The factors of C and yfgl in (6.73) are not present in (6.72). The factors of C can be 
absorbed into a redefinition of the (Tj, but the absence of the factors of y/~g^ will make the 
(n|m) 7^ (0|0) terms blow up in the double-scaling limit. 

2. The power of / in (6.72) should reproduce the power of z in (5.120) in the double-scaling 
limit. We see from (6.74) that this is essentially what happens but that, in the present 
form, the double-scaling limit of / also has an unwanted (7s-dependence. 

Both of these problems can be solved by the following somewhat unconventional change of vari- 
ables: 

ai = ^,ai-{72\gs)--^^'^\ (6.79) 



'Recall that 5 = {n-\- m) mod 2. 
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= V9^a2-(72Ag,)+6'^^"^ (6.80) 

We have discussed in section 5.2 that one is allowed to make criO"2-dependent changes of variables 
in a two-parameters transseries. The somewhat surprising fact in here is that we now find a 
transformation which is also (^^"dependent. In the Painleve I case, the expansion parameter in 
the transseries was z. Thus, in that case, one was not allowed to make z-dependent changes of Oi 
for the simple reason that this would spoil the Painleve I equation, itself a differential equation 
in z. However, in the present quartic matrix model case, although the expansion parameter in 
the transseries is g^^ the string equation is not an equation in — it is an equation in x, for 
which gs is a parameter. For this reason, a (7s~dependent change in cjj does not spoil the string 
equation, and we are in fact allowed to make the above change of variables. 
Inserting the new variables into (6.72), we find that 



+00 +00 +00 , J., . ^ j^(n-m)o-icr2 

„=0 m=0 g=l3„„, ^ ^ 

(6.81) 

and we see from (6.73-6.75) that if we define the transseries parameters for the Painleve I equation 

as 

f^i,Pi = ^, (6.82) 
we indeed get the correct double-scaling limit, u{z), 

+00 +00 +00 , 

Ex — ^ \^ — ^ / \ A 5/4 ('nlm^rnl 10g + 5(n + m)-4 , 4(n-m)CTi pjCT2 pi 
E E <,Pi<Pie-('^-"^^^-^' 4in+l-^ ^ • (6-83) 

n=0 m=0 g=/3nm 

Notice that the only difference between this expression and (5.120) is that the present formula 
has a lower starting genus: the first term in the ^-sum is the one with u'^p"^^^'!^_^^. However, as 
we have defined all coefficients with genus smaller than the starting genus to be identically zero, 
this is not a problem (in principle, we could have started all ^f-sums at — oo). 

Choice of Resurgent Variables 

Having indentified the correct double-scaling limit of the transseries TZ{x), we can now test its 
resurgent properties. Recall that also for non-linear difference equations there exists a suitable 
transseries framework [63] for which one may develop resurgent analysis in a fashion similar 
to what we have worked out in section 4 (although the literature on this class of equations is 
considerably smaller than the one on non-linear differential equations). However, the difference 
equation we address in this problem, the string equation, arises from a matrix model set-up and, 
in particular, has very sharp physical requirements on what concerns double-scaling limits. In 
other words, our difference equation must relate to a differential equation, in a prescribed way, 
also at the level of resurgence. This will introduce some new features as we shall now see. 

Indeed, and as discussed previously, the transseries resurgent structure of the string equation 
is highly dependent upon a judicious choice of variables (the ones which properly implement the 
Painleve I double-scaling limit). As we discussed above, the naive choice of variables for TZ{x), 
i.e., the choice of variables that one would consider natural from a purely finite-difference string 
equation point of view, is not the one that leads to the correct double-scaling limit — for this, one 
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further needs to make the (7^ -dependent change of variables, from (cti, (T2) to {ai,a2), defined in 
(6.79-6.80). As a result, we get a new transseries representation for TZ{x); schematically 

^ a^a^ gl R^J"^""^ = Yl ^r^a^^r""""^ log^72A5s)«("l")<'-^ (6.84) 

n,m,g n,m,g,k 

In this new representation, different powers of Qs and entirely new powers of logg's appear^'^. 
Applying the resurgent formalism using the standard expressions for the alien derivatives can 
only give correct large-order formulae in one of these cases. We shall thus make the obvious 
assumption: we will assume that the correct representation is the one on the right-hand side 
above, which is the one leading directly to u{z) in the double-scaling limit. In the following, we 
shall find ample evidence supporting this assumption. 

Tests of Resurgence: Perturbative Sector 

As a first test of resurgence, let us study the large-order behavior of Rf'^^\ Since 

jj(n|0)(0> ^ ^(n|0)^ i^(0|m,)(0> ^ ^(0|m)^ (g_85) 

the result takes essentially the same form for either hatted or unhatted components. Applying 
our resurgent formalism to the i?-transseries, and making the above substitution, one finds the 
large-order prediction 

4oio)(,) ^ y (^lI ng - k/2) ^ n9-h-k/2) (,,0) , 

This result is valid for even values of 5, so that R^g^'^^ is defined. In here, our change to the 
hatted components has still played a role: if we had applied the resurgent formalism directly to 
the i2-transseries, we would not have found the terms of k/2 in the gamma function and in the 
power of A. Notice that this issue was already present in [12], albeit implicitly: in there, this was 
solved by leaving a 5^ -dependent factor in the Rg , leading to the somewhat counterintuitive 
result (equation (3.50) in that paper) of a ^g-dependent sf'^ Stokes factor. To the contrary, 
our present formalism leads to large-order formulae which are (^^"iiidependent — a more natural 
form for a quantity describing the coefficients in a ^^-expansion. Moreover, as we shall see, this 
procedure can be straightforwardly applied to all generalized instanton sectors, including the 
ones where the relation between the R and i?-coefficients is more complicated. 

We now wish to test the large-order formula (6.86). The first prediction we get from it is 
that the leading large-order behavior of R^^^^ is 

We have tested this behavior in a computer for a large range of x (or, equivalently, r) and A, 
and found that it was completely consistent (up to at least 10 decimal places in all cases) with 
a value of sf''^ equal to 

S<»' = i/A. (6.88) 



^"We have now labeled the coefficients of the logga terms with (k) to avoid confusion with the (still present) 
coefficients of the logz terms, which we are labeling with a [k] index. 
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This formula equals (3.50) in [12] if we take into account the removal of that was discussed 
above, as well as the definition of A in that paper which differs from ours by a factor of 2. 
To illustrate these tests, let us set A = 1/2 and plot the large-g' values of 

^-^^^R^Hx) (6.89) 

for a sequence of equally spaced values of r, defined as a function of x in (6.31), between and 
its double-scaling value r^g = 1/A = 2. As before, we obtain very precise large-5 values by 
calculating the above expression for values up to 5 = 50, and then applying a large number of 
Richardson transforms (10 in this case) to remove ^""-effects. The result is given by the blue 
dots in figure 6; the red line in that graph represents the expected result of 

d(i|o)/^\ - 



ivr 



° ^ ' 7r3/2(3-Ar)i/^(3-3Ar)i/^ ^ ' 

where we have inserted the explicit expression for given in (6.42). We see that the large- 

order results perfectly match the predicted values. At the smallest value of r, the error is 0.002%. 
This error is mainly due to the fact that, for small r, a very large amount of -Rg*^'"^ data is required 
to get good Richardson transforms. The error quickly decreases as r increases; from r = 0.18 
onward, it becomes stable at around 10~^^%. 

As a further test of the large-order formula (6.86), we could now study the next-to-leading 
order behavior in g~^, arising from R^^^\ and so on. However, as discussed earlier in section 5.4, 
for the Painleve I case, we can actually test all perturbative corrections at once by Borel-Pade 
resumming them and going straight to the corrections. That is, we calculate the quantity 

0(0) +00 -p / r l\ 
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Figure 7: The 2^^ corrections to tlie large-g beiiavior of the Rg , expressed in terms of the quantity Xg 
on the left-hand side of (6.92) (blue dots). The red line indicates the predicted value from the right-hand 
side of that equation. We have set A = 1/2; the variable along the horizontal axis is r. 



by Borel-Pade resumming the second term as an expansion in g ^, and then test the prediction 
that we get from (6.86): that the large-order behavior of this quantity is 

_.(2A(^ ~ -- (5f))'4^l°)(x) = ^, (6.92) 

r(<7-l) ' ttVW 27r2(3-Ar)i/2(3_3Ar)3/2' ^ ^ 

with R^Q^'^^ given in (6.46). Note that here we have also included a factor of — i to pick out the 
imaginary part of Xgi as in the Painleve I case, the correction in the large-order formula is 
purely imaginary, due to the fact that it comes from integrating around poles in the Borel plane 
with a given choice of itie-prescription. 

In figure 7, we plot the large-order quantity on the left-hand side of (6.92), calculated using 
the usual Richardson transform method, as well as the expected result on the right-hand side 
of that equation (the red line in the plot). We have once again set A = 1/2 and varied r. The 
large-order data starts at a value of r = 0.22; for smaller values, the amount of data required 
to get a good large-order approximation is too large to be calculated in a reasonable amount of 
time. The upper bound on r is again its double-scaling value rds = 1/A = 2. Akin to before, we 
find a very good match between the data and the prediction. For the smallest value of r, where 
the amount of data is barely sufficient, we find an error of 20%. The error reduces quickly as the 
value of Xg becomes larger: when r = 0.34 the error is already less than 1%, and it becomes as 
small as 0.007% near the double-scaling limit. 

As a final remark on the validity of the large-order formula (6.86), let us take its double- 
scaling limit using (6.73) and (6.75). After some straightforward algebra, one finds 

^(0|0)^. 1 VfgWc-i'l' ^(g-^/^) ^ T{g-h-k/2) MO) 

This formula agrees with the Painleve I large-order formula (5.58), provided that the Stokes 
constants for the quartic matrix model and for Painleve I are related by 

'S'i°QMM = ^ ^i%v (6.94) 
Inserting the values (6.88), (6.76), (5.84) for these constants, we see that this is indeed the case. 
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Tests of Resurgence: Instanton Sectors 

Having tested the large-order behavior of the perturbative part of TZ{x), we now want to switch 
to its (generahzed) instanton components, as this is where new Stokes constants and "back- 
wards/sideways resurgence" appear. The large-order behavior of the one-instanton coefficients 
still only depends on sf'^ (at least perturbatively in g^^), so the simplest coefficients to 
study for our purposes are the two-instantons coefficients, R^^^^\ 

Thus, our first task is to derive a large-order formula for these coefficients. For this, it turns 
out to be essential to use the hatted representation of the transseries given in (6.84). The reason 
is that the large-order behavior of the (2|0)-component of any transseries depends, through 
"sideways resurgence", on its (2|l)-components. The latter components contain logarithms, and 
so it is essential that we correctly include the logf^^ terms to get the correct large-order formula. 
After calculating the resulting large-order expression for the i?-transseries, we can then translate 
the result back to the i2-components using the relation (6.85), as well as the relation 

i?(2|i)(i>=_^^(i|0)^ (6.95) 
y g y 

that can be read off after expanding both sides of (6.84). Doing all of this carefully, one finds 
the following large-order expression 

^(2|0). ^ ^ 3^ • ^^S-^-'^) , 



{A{x)) 

2- (A(x))-'^-i 



4 

h=-~l 



(_1)»A5{°' ^, ,>/.„(l|0), ^ r(g-ft+l).B,2,,4(,)(9-ft+l) 



+ 



where Bs{a) is the shifted digamma function defined in (4.56), and we wrote the answer in terms 
of the purely imaginary combination 

5f ) = i5i\) + ^ 5f \ (6.97) 

which is also (compare against expressions such as (5.56)) the coefficient determining the large- 
order behavior of the "conjugate" coefficients B!"g^'^\ 

As we did several times before, (6.96) can now be tested on a computer. Doing this, we 
found that the above large-order formula holds and that, for a wide range of A and r, up to 8 
decimal places it is the case that 

"^l.QMM = (J ' (6.98) 

with the numerical value of pj given in table 2. As an illustrative example, we once again set 
A = 1/2 and evaluate the quantity 
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Figure 8: The large-g behavior of the Rg , expressed in terms of the quantity Xg on the left-hand 
side of (6.100) (blue dots). The red line indicates the predicted value from the right-hand side of that 
equation. We have set A = 1/2; the variable along the horizontal axis is r. 



where iJ^^'*^^^^^ ^^^^x) is the optimal truncation of the first three terms on the right-hand side 
of (6.96). To leading order, we expect this quantity to grow as 



1 



(6.100) 

In figure 8, we plot the large-order quantity on the left-hand side of this equation as blue dots 
and the prediction on the right-hand side as a red line, for values of r between r = 0.22 (where 
we can generate just enough data) and the double-scaling value r^s = 1/A = 2. We see that 
the results once again match the prediction very nicely. We have included explicit error bars 
(estimated by comparing the results for two consecutive values of g) to indicate that the results 
for the lowest values of r are still within the expectation. From r = 0.5 onwards, the error due 
to lack of data is negligible, and we get results which are correct up to 8 decimal places. 

As an extra check on the validity of the large-order formula (6.96), we can calculate its 
double-scaling limit using (6.73-6.75). It turns out that most of the logarithmic terms coming 
from -R^f' and S72AA cancel, leaving a single term proportional to logApi. All other terms 
reduce straightforwardly to terms involving the Painleve I coefficients, and in the end one finds 

(2|o)[o] _ v"„(3|o)[o] r (ff - /i - ^) {-iys[°l^ \^(_'i\h(2\im r (g - /t - ^) _ 

h=0 Apj h=0 Apj 

+00 



izl^!:?li£I V^r_n'^,,(i|o)M r (g - /i + ^) • Bap, [g-h+l] 



■'^)^ CSfj^yy^^ (^|Q)JO] r(g-/l+^) 



+ 



27ri ' .g~h+\ ^ ^ 

In this expression, everything is written in terms of Painleve I quantities, except for the com- 
bination CS\ ^yiy^ in the last term. If we now directly apply the resurgence formalism to the 
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2-instantons component of the Painleve I transseries, we find precisely the same large-order for- 

~(2) ~(2) 

mula, but with CS^ q^j^ replaced by pj. The two large-order formulae thus exactly coincide 
when (6.98) is valid, providing a good extra check on the validity of that equation. 

As a final test, we study the large-order behavior of the generalized (l|l)~instanton coeffi- 
cients, rIj'^^^ ■ Applying the same techniques as above, we find the large-order formula 



+ 



A5f|) - (^10) r(g-/.+ ^)-i372AA(.)(g-/^+^) 



h=0 
Jl) oo 



In this formula, a new Stokes constant appears, s[^\ We have checked by computer that, up to 
4 decimal places, it equals 

"^l.QMM = ~c~' (6.103) 

Furthermore, as we did before, one can also check that this result precisely leads to the correct 
Painleve I large-order formula in the double-scaling limit. 

For a graphical illustration of the S^^ tests, let us choose A = 1/2 as usual and calculate the 
quantity 

where Rg (x) is the optimal truncation of the first two terms on the right-hand side of 

(6.102). To leading order, we expect this quantity to grow as 

(-l).<i!^X,(.)~^4''°'W. (6.105) 

r (5 + 2) 1^ 

Figure 9 shows the large-order quantity on the left-hand side of the above equation as blue dots, 
and the prediction on the right-hand side as a red line. The variable r ranges between r = 0.10 
and the double-scaling value r = l/A = 2. In spite of the fact that the coincidence is not perfect 
(due to a lack of -Rg^'^^ data, the production of which consumes large amounts of computer time), 
the results still match the prediction within a few percent^^. 



Moduli Independence of Stokes Factors 

We have now explicitly calculated three Stokes factors for the quartic matrix model and we have 
seen that, up to high accuracy, they satisfy 

0(0) _^c.{0) qW _ ^-1 c.(i) c(2) _^-ic(2) 

'^l.QMM ~ ^ '-'l.PI' '-'l,QMM ~ '-'l.PI' '-'l.QMM ~ '-'l,PV {O.lUU) 

^^The reason that we can actually calculate s[^^ itself to higher precision is that, for that calculation, we can 
also take the optimal truncation of the third term in (6.102). 
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Figure 9: The large-g behavior of the Rg^^^\ expressed in terms of the quantity Xg on the left-hand 
side of (6.105) (blue dots). The red line indicates the predicted value from the right -hand side of that 
equation. We have set A = 1/2; the variable along the horizontal axis is r. 

In particular, since C ~ A~^/^, the above quartic matrix model Stokes factors depend on the 
parameter A. But this C-dependence is somewhat artificial: as we saw in (5.101-5.102), one can 
rescale Stokes factors by simply rescaling the parameters ai with some factor, c, resulting in 

5f ) ^ c2-2^-^5f , 5f ^ c^-^'^+^sf \ (6.107) 

Thus, by choosing c = C^^, we can actually make all three Stokes factors A~independent, and 
exactly equal to their Painleve I counterparts. Note that this is nothing but the scaling (6.82) 
that produces the Painleve I transseries solution u{z) out of the quartic matrix model transseries 
solution TZ{x), in the double-scaling limit^^. 

The statement (6.106) is much stronger than a statement just about the double-scaling limit: 
it says that, up to a trivial C~dependent rescaling, the quartic matrix model Stokes constants 
we have calculated are independent of the parameters of the model. That is, their value at any 
point in parameter space equals their value in the double-scaling limit — and hence the value of 
the Painleve I Stokes constants. That this is the case for S^^ alone is not too surprising: one can 
always choose a c in (6.107) in such a way that s[^^ becomes independent of the parameters. But 
that this is also the case for the other Stokes constants is indeed quite interesting. Nothing in the 
resurgence formalism seems to prevent these Stokes constants from depending on A, or — as we 
shall see in more detail in the next section — on some combination of A and the 't Hooft coupling 
t. The only consistency requirement is that the Painleve I Stokes constants are reproduced when 
taking the double-scaling limit, which is expressed in here by choosing appropriate resurgent 
variables that allow for the matching of transseries solutions off-criticality and at criticality, as 
in (6.84). Here, we find that this requirement is fulfilled in the simplest possible way: by having 
off-critical Stokes constants which are fully independent of the parameters. 

It would be very interesting to understand why the quartic matrix model Stokes constants 
that we have found are parameter-independent in the above sense. We have not been able to 

^^We could of course have chosen to absorb this scaling already in (6.79-6.80). The reason for not doing this 
was that it would have spoiled the simple relation (6.85) between the hatted and unhatted representation of our 
transseries. 
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find a compelling argument for this, but of course it is very natural to conjecture that the above 
is not a coincidence, but that it is actually true for all Stokes constants. That is, we conjecture 



This gives us conjectured values for many new quartic matrix model constants: up to the above 
A-dependent rescalings, they should be equal to the Painleve I values that we reported in table 
2. Combined with the further conjectures in footnote 36, this gives us for example the conjecture 
that the exact value of the index (0) Stokes constants is 

S'^U = " • (6.109) 

It would be very interesting to further test these conjectures, fully understand them from a 
physical point of view, and put them on a firm resurgent analysis mathematical footing. 

6.5 The Nonperturbative Free Energy of the Quartic Model 

Having fully constructed the two-parameters transseries solution for TZ{x), our final task is 
to translate this solution into an expression for the free energy T{t,gs)- In section 6.2, we 
already briefiy discussed how to do this. We saw that the Euler-MacLaurin formula leads to the 
expression (6.19), which we repeat in here for convenience: 



Ht,9s) = -(2log^- log — 



+ ^ / dx t - x log -^-^ + 
li Jo 

(6.110) 



2,-2 B2, d^^-^ 



(2^)! dx^s- 



x=oJ 9. 
{t - x) lo, 



x=0 



In applying this formula, we can of course choose any parametrization of TZ{x) we wish, and one 
will thus end up with the corresponding parametrization of the free energy J-{t, Qs)- To get a good 
double-scaling limit, in this section we shall once again work with the "hatted representation" 
that was introduced in (6.84). But do notice that, in order to avoid cluttering the notation too 
much, we will not put any hats on the corresponding coefficients of J-{t,gs)- 

As we saw in section 6.2, the above expression is valid for the full two-parameters transseries, 
meaning that we can apply it both in the perturbative sector and in the (generalized) instanton 
sectors. We shall next discuss its application in these different sectors. 

The Perturbative Sector 

In the zero-instanton sector, the above formula was already used in [25, 15] to compute the 
first few genus-^ free energies. There, it was found that the result takes its nicest form when 
expressed in terms of the 't Hooft coupling constant t = ggN and a variable denoted by a^. This 
variable was also introduced in (6.12); it determines the end-points of the eigenvalue cut and is 
defined as 

= ^ (l - Vl - 2Xtj . (6.111) 

Note that is also essentially equal (up to an exchange t -H- x) to the variable r introduced in 
(6.31). In terms of t and a^, it was conjectured in [25] and confirmed up to genus 10 in [15] that 



-82- 



the perturbative expansion coefficients of the free energy are of the form 



-^r^(t)= -^^(^)- (6-112) 

^ t9 {2t - a^y^'^ 

where is a homogeneous polynomial^^ of degree '^^^"^ ct^ ^'^^ ^• 

In the above expression, the factor of can be combined with the prefator gi to recover 
the original factor of iV^^ appearing in the large expansion of the matrix model free energy. 
Apart from this factor, we see that T^^'^^ (i) in (6.112) above, which could apparently be thought 
of as a natural function of t and a^, is actually just a function of tjc? . In other words, the 
perturbative free energy components do not depend on the two separate parameters o? and t 
(or, equivalently, A and t), but only on a single combination of the two. As we mentioned in 
the previous subsection, this result could have been anticipated: the coupling constant A in the 
quartic matrix model potential y(M) can be absorbed into t by rescaling the variable M . We 
shall see that this pattern naturally extends to the full transseries solution. 

Specifically, the first few T^^^^ {t^ are^^ 

^iTV) = ^ (at" - lOfa" + a' + \ie- log , (6,113) 

^r'(0 = 4log(^5^), (6.114) 



720*2 (2t - a^f 



{t - a^f (17260*5 - mMt^a^ - 2925*^04 + gSSt^a^ _ i35tQ,8 + 9q,io^ 

9072*4 (2t - a2 



We next want to investigate how these results extend to the (generalized) instanton sectors. 
The Nonperturbative (n|n) Sector 

To calculate the higher (generalized) instanton contributions to J-{t,gs), it is more convenient to 
use the result of the Euler-MacLaurin formula in the form (6.24), 

T{t + gs) - 2F{t) + F{t - gs) = log (6.117) 

A first consequence of this equation is that the instanton action A{t) of the free energy equals 
the instanton action for TZ{t) [12]. We constructed this action as a function of r and A in (6.36); 
expressed in terms of and t it takes the form 

A{t) = -tarccosh [^^^^) + j^^^ \/I^^^^3^. (6.118) 

Note that, once again, A(i)l gs is a function of the single combination of variables c? jt. 



^^As usual, there are exceptions at low genus, g = —2, in this case, where logarithmic contributions appear. 
®*In our present conventions, these results differ by an overall minus sign from those in [15]. 
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For the free energy, we therefore make the following two-parameters transseries ansatz^^ 

+00 +00 +00 
•^(*'5s) = J] J]cTrare-("-™)-^W/^» giJ^t^'^Kt), (6.119) 

n=Om=0 a=B^ 

where, as usual, is the lowest g for which a nonvanishing term is present. From the calcula- 
tions we shall present below, it is a straightforward exercise to calculate that 

n = m: = ^^±|^, (6.120) 

m = 0,n>0: /3;^^ = , (6.121) 

^-p n — m + 2 
m > 0,n > m : = , (6.122) 



and by symmetry = Now, it is a matter of plugging (6.119) into (6.117) and expanding 
in (qs, (Ti, (T2) to obtain equations for the J-g^^"^\t). When n = m, these are differential equations: 
one obtains 

Here, we have denoted the two-parameters transseries representation of the right-hand side of 
(6.117) by 

logM^£(t), (6.124) 

and C{t) has a transseries expansion completely analogous to (6.119). Expressing the 4"'"^(t) 
m terms of the R^J'^"'\t) is once again a straightforward exercise in Taylor expanding functions 
of transseries. The sum on the right-hand side of (6.123) is infinite, but only a finite number of 
terms contribute for any given choice of n, as g in £g"'"''(t) is bounded from below. 
Solving the above equations for n = 1, one obtains for the lowest two genera, 



^(ih),, V2t - aW2t + a2 t (t - a^) / 2^3 {it - a^) + 6V2t - a V2t + aA 



J'r'it) = ^ L^. 6.126 

6^/3a4(2t-a2)5/2(2i + a2)3/2 

Similarly, for n = 2 one finds 

^(^l^)(t) - i^log ( 1 (6 127) 

(2\2). . (t - a^Y (1696*6 - Siet^a^ _^ i896t''a^ - 5408*^06 + 2229t2a8 + 516ia^° + ISOa^^) 



486a8 (2t-a2)^(2t + Q2)3 



®^Our conventions diflFer from the usual "perturbative" ones, where J-g{t) denotes the function multiplying g^^ ^. 
When including instanton sectors, it becomes more convenient when the subscript of J-g"^^\t) simply indicates 
the power of Qs that it multiplies. Thus, our g should be thought of as an Euler number, not a genus. 
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In both cases, the g < results are exceptional, with logarithmic contributions. For all strictly 
positive g, one finds the following general structure of the solution: 

(t ^V" 

where ^^"'"^(t) is a homogeneous polynomial of degree 3^. Note that, formally, J-g^^^\t) in 
(6.112) is also of this form if we take the corresponding function Vg^^'^\t) (which is now no longer 
a polynomial) to be 

(t-a2)^'+ij2t + a2)33/2 

~t9' 



p(0|0)(t) = ^ ^-^^ 1 Sg{t). (6.130) 



In appendix C, we present some higher-genus examples of Vg^^^\t) for n = 1, 2. 
The Nonperturbative (n|m) Sector 

When n ^ m, the Euler-MacLaurin formula in the form (6.117), 

^(t + gs) - 2F{t) + F{t - gs) = log (6-131) 

gives, upon expansion in (5s, cri, 0-2); a set of algebraic equations^^ for Fg^^'^\t). For example, 
for the lowest two orders, one finds 

fP^^ = \ sinh-2 (^^'j £H-), (6.132) 



-n(n\m) 
0^ +1 



where i = n — m. Solving these equations is now straightforward (see also [12] where this was 
already done for the (l|0)-sector), and we find for example the one-instanton results 

^ ^ ~ 35/4a2 (2i _ «2)5/4 + „2)l/4 ' ^^'^^^^ 



^(110).. ^ (t -g^)^/^ (40^3 _i2t2c.2_2lto4_ 10,6) 

^ ' 6V2 33/4 a2 (2t - a2)i5/4 (2i + ^2)7/4 ' ^ ■ ) 



which agree with the results in [12], and the two-instantons results 



[t - a2)3 (736^4 _ io96t3a2 _^ 554^20,4 _ 25340^ + 22q8) 
162a4 {2t - a^f {2t + a^f 



By "algebraic", we mean that J-i"!"''^ itself occurs algebraically (and even linearly), so that no integrations 
are needed to solve the equation. Derivatives of lower J^^^ \t) still appear. 
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In appendix C, we present some higher-genus results, as well as some results for the generalized 
instanton sectors (2|1), (3|1), (3|2) and (4|2). Their logarithm-free part (we will discuss the 
logarithmic terms in a moment) satisfies the general structure formula 

' (a2)"+-(2t_a2)59/2(2t + a2)(39-5)/2 5 w, 

where Pi"l'")(t) is a homogeneous polynomial of degree (6(7 + 5 — 4) /2 (recall that here, as usual, 
S = {n + m) mod 2). This expression should be compared to the very similar result (6.51) for 
rIP^"^\x). Also note that, apart from the different degree of the polynomial, the result (6.129) 
for n = m is nothing but a specific case of the above equation. 

The Logarithmic Sectors 

As is familiar by now, whenever n > and m > 0, the J-g^^"^\t) contain logarithmic terms. 
Once again, these logarithmic sectors do not contain any new information: one finds that when 
n ^ m, J^g"'^"^\t) is of the form 

mm(n,m) , , 

_^(n|m)(^)^ J-(nMW(i).log^ (6.139) 



with 



JT{n\m)[k]^^^ ^ i_ ^ A (ra - m) y T-in-k\m-k)\0] _ 



(6.140) 



The function f{t) is essentially the same function as before (see (6.49)), but now conveniently 
written in the variables and t, 



81 ( a^-2t)'(a2 + 2tf 
16 aS {a^-tY 



fit) = o ,^4 ■ (6-141) 



For readability reasons, we have left some factors of A explicit in the above expressions, but in 
principle, they should also be rewritten in terms of these variables, that is 

A=^^, 6.142 

which is the inverse of (6.111). As before, one can also choose to sum all the logarithmic sectors 
resulting in the closed form 

+00 +00 +00 ^ N V j^(n-m)(Tio-2 

^it) = E E <a^e-(— Yl 5f -^i-i'-^'^kt) • Uiiw ' ^^•^'^^^ 

n=0 m=0 Q=0'J^ ^ 

for the two-parameters transseries. In here, we have introduced the shifted starting exponent 
Pnm = Pnm (n>m = 0), P'Zn = l^L + (n>m>0), (6.144) 

extended by symmetry to the cases where n < m. The reason for the shifted exponent in the 
cases where m > is that in these cases, J-g^^"^\t) starts off with a purely logarithmic term. 
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Double— Scaling Limit 

Because we started with the "hatted representation" for TZ{x), which gives the Painleve I solution 
u{z) in the double-scahng limit, it is very natural to expect that the corresponding free energy 
J-{t,gs) also gives the (2, 3) minimal model free energy F{z) in the double-scaling limit. Indeed, 
as we discussed in detail for TZ{x), the factor 



( fit) 
V5184A252 



A 



(n~m)aia2 

(6.145) 



nicely reproduces the structure of the log z terms in the Painleve I solution. Thus, all we need to 
check is that the coefficients J-p""'™^^*^' (t) have the correct double-scaling limit. Indeed, we have 
checked that in this limit, and for all of the examples presented in appendix C, 

oo 

Sf-^i"''"^'"'^ ^ F^"l"^(2), (6.146) 

with f("I"*)(z) given in (5.128-5.136). This once again underlines the fact that the hatted 
transseries representation is the correct representation to study when one is interested in the 
double-scaling limit. 

Stokes Constants for the Free Energy 

In the case of the (2,3) minimal string, we found simple proportionality relations (5.110-5.111) 
between the Stokes constants for the free energy F{z) and those for the solution, u{z), of the 
Painleve I equation. We were further able to derive these relations analytically, because the map 
between u{z) and F{z) (a double integration) is a very simple and linear map. 

Unfortunately, for the quartic matrix model, the situation is a whole lot more complicated. 
The Euler-MacLaurin formula (6.110) is very involved and it is difficult to deduce from it a 
direct relation between the large-order behavior of the R^J^^^^^^^'^x) and that of the J^g"''"^^^'^^ (t). 
Moreover, the computer generated data we have in this situation is insufficient to check or derive 
such a relation numerically, beyond the first Stokes constant. 

Nevertheless, one can still make an educated guess as to what the result could be. It was 
found in [15] (see equation (4.15) of that paper), both from a spectral curve analysis and using 
numerical results, that the large-order behavior of the perturbative series J-g^^^\t) is determined 
by the function 

(t - a^) 

Mt) = ^ ^77 777- (6.147) 

^' 33/4V^(2t-a2)5/4(2t + a2)V4 ^ ^ 



In our notation, this function corresponds to the combination 

wW = ^Sm--^iT(*)- (6.148) 
Thus, comparing (6.147) to (6.134), we find that 

We see from this that 'S'^'q^m exactly equal to the Stokes constant S'^'qmm 
transseries, presented in (6.88). This is very similar to what we found in the Painleve I case: 
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in the correct parametrization, the Stokes constants for the free energy F{z) are equal, up to a 
factor of to the Stokes constants for the corresponding solution u{z). Thus, we may make the 
natural guess that the same pattern holds for all Stokes constants of the quartic matrix model 
free energy, 

qiO)F _«2c(0)ii n(0)F _ n2qiO)R ((\^^^\\ 

Note that this guess can also be viewed as extending the parameter-independence of the quartic 
matrix model Stokes constants for TZ, to the corresponding Stokes constants for J^: it essentially 
states that, up to a trivial reparametrization, the quartic matrix model Stokes constants equal the 
Painleve I Stokes constants. It would be quite interesting to prove (or disprove) this statement. 



7. Conclusions and Outlook 

In this paper we have hopefully made a strong case for the existence of new, previously unno- 
ticed, nonperturbative sectors in string theory. The full structure we have uncovered was first 
anticipated in [13], by studying the asymptotics of instantons of the Painleve I equation, and 
first discussed, from a physical point of view, in [14]. But what exactly are these sectors? We 
hope to report on this question in upcoming work, but let us also make a few remarks herein. 

The physical instanton series is simple to understand: it corresponds to standard matrix 
model instantons [16, 17, 15, 23] which, in the double-scaling limit, become ZZ-brane amplitudes 
in Liouville gravity [20]. As we shift our attention to the remaining sectors the first thing 
one notices is that the structure of the transseries solutions we have addressed, where purely 
"generalized" instantons have an overall minus sign in front of the instanton action as compared to 
standard instantons^^, could seem to point towards understanding these new sectors as ghost D- 
branes [64] (or, in the matrix model context, their counterpart of topological anti-D-branes [65] 
as dictated by the correspondence in [3]). Indeed, these ghost D-brane sectors display this exact 
same feature as they have an overall minus sign in front of the Born-Infeld action [64] (also see 
the discussion in [14]). This is an appealing picture: for instance, in the examples we have studied 
the free energies FH^), with as many instantons as purely "generalized" instantons, were found 
to have a resulting perturbative series which is a series in the closed string coupling constant g^. 
However, both ghost D-branes or topological anti-D-branes have one further property [64, 65], 
which is that their free energies must satisfy 

^(n|m) ^^(n-m|0)^ n > m. (7.1) 

But this is a property we may explicitly check within our examples, and it is a property which 
is certainly not satisfied. To illustrate, let us recall in here the case of the Painleve I equation 
where we found 

^.211 w ^ 71 15 2999 25 25073507 35 2705576503 45 , 

F^^^^'iz) = z'^ H ^z~^ z"^ H ^z~^ ,(7.2) 

^ ' 864 18432^3 191102976 6794772480^/3 ^ ^ 

^^iinw ^ 1 5 37 15 6433 25 12741169 35 ^„ ^ 

F^^I^Uz) = z~8 H ^z~~8 z~~s H ^z~^ . (7.3) 

^ ' 12 768^/3 294912 283115520^/3 ^ ^ 

It is simple to see that these two sectors are not proportional to each other. Furthermore, one can 
also show that there is no reparametrization transformation that can achieve such proportional- 
ity. This is a straightforward consequence of (5.28) which states that, upon reparametrization, 

^^Of course this is only the case in our present setting of a resonant two-parameters transseries. When deal- 
ing with general multi-parameter transseries, required in the solution of matrix models with more complicated 
potentials, or in the solutions of the minimal series coupled to gravity, this simple scenario will no longer be true. 
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the only possible change of F^'^^^\z) is by a multiple of F^'^\^\z). Thus, if is not a 

multiple of in one representation, that statement is automatically true for any other 

reparametrization. Further notice that using the transseries structure of the free energy as in 
(5.123), where the transseries parameters also appear exponentiated, does not change this con- 
clusion. Indeed, the exponentiation (5.124) is just a convenient way to rearrange the logarithmic 
sectors, which can always be reversed (by expanding the exponential). In this case one would 
then apply the aforementioned argument to each separate logarithmic sector with the exact same 
conclusion. As such, although we cannot at this stage state what the new nonperturbative sectors 
are, it seems we can state what they are not. 

Another pertinent question is: why have we never seen these sectors before? The short 
answer is, of course, that two-parameters transseries were never addressed in a string theoretic 
context prior to [13]. Only by addressing the question of what controls the asymptotic behavior of 
rnu/ti-instanton sectors can one realize that indeed the familiar physical instanton series cannot 
be the full story. In fact, most large-order analyses have always been concentrated upon the 
leading asymptotics of the perturbative sector [31]. But, as we have shown at length in this 
paper, if we want to address harder questions than that, in the string theoretic nonperturbative 
realm, then the full multi-parameter transseries framework is indeed required. 

On the other hand there are examples of exactly solvable models, where full nonperturbative 
answers have been computed. Should any of these expressions have shown these new sectors? Of 
course in order to see them one would have to know what to look for. But when one rewrites one 
of these exact nonperturbative solutions in terms of semi-classical data, one usually does so only 
for real solutions around positive, real coupling, and in the one-parameter transseries framework! 
Let us briefly discuss the construction of real solutions, trivially generalizing a discussion in [12] 
to an arbitrary one-parameter transseries of the type (2.17), 



+ 00 

F{z,a) = ^a^e-^^'<^n{z). (7.4) 

n=0 

A real solution starts around positive real coupling z G M"*". But this is a Stokes line and we need 
to be careful in constructing such real solution. For instance, upon Borel resummation, either 
S+F or S-F, will display an ambiguous imaginary contribution to the solution which needs to 
be canceled, i.e., one needs to set^^ ImF(z,a) = 0. As it turns out [12], ImF{z,a) = if and 
only if Ima = ^Si. As such, and as long as the instanton action is real, a real solution can be 
constructed by considering [12] 

z, a) = S+F (z, a - ^Si] = S^F ( z, a + ^Si] , (7.5) 



2 J V 2 , 

where the transseries parameter in the expression above is now o" G M, and where the second 
equality follows trivially from the Stokes transition (2.41) 

S+F{z,a) =S-F{z,a + Si). (7.6) 

Expanding, it immediately follows 

Fk{z, a) = Me F^°^ (z) + a Me F^^^ (z) + (^a^ - ^Sf^ Me F^^) (z) + ■ ■ ■ . (7.7) 



'Notice that around the 6 = Stokes hne one has Imo = ^ (5+ — 5-) and Reo = | (5+ + 5- 
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Two things are to be noticed. The first is that indeed real solutions display instanton corrections 
(even if a = 0) . This is simply because the string equation (be it the Painleve I equation or the 
quartic string equation or any other) is non-linear and, although S+F or S-F may be solutions, 
their sum is, consequentially, not a solution. Indeed, their sum can only become a solution 
once we correct it appropriately by accounting for higher instanton corrections. The second 
point, however, is that this instanton expansion only includes information concerning Si, not 
about any of the other Stokes constants. This is to say, as long as we consider the expansion 
in semi-classical data around the (natural) 9 = Stokes lines, we shall find no indication of 
the multi-parameter transseries sectors. Searching for signs of these new generalized instanton 
sectors within nonperturbative answers must thus start by properly addressing what type of 
expansion one wants to do — as shown, the standard one will not do. 

In summary, we believe the most pressing question begging to be addressed is to fully un- 
derstand, from a physical string theoretic point of view, the generalized instanton series. As 
discussed, D-branes only yield information on a limited set of Stokes constants and, if one 
is to address nonperturbative questions where all Stokes constants play a role, some informa- 
tion is missing. Examples where all Stokes constants would be required involve general Stokes 
transitions — even if we are just addressing the perturbative series, Stokes transitions along 9 = n 
will require Stokes constants which, at this stage, have no first principles derivation. For instance, 
within the setting of the quartic model, one could imagine rotating the string coupling in the 
complex plane from the positive to the negative real axis. The saddle configuration would then 
change, from the one-cut spectral geometry we addressed in this paper to a two-cuts spectral 
curve. This change of background may be implemented within our framework — the transseries 
does provide the complete nonperturbative answer — via a Stokes transition, but in order to ex- 
plicitly construct the perturbative free energy around the new background, given the original 
one, we are still missing analytic expressions for the Stokes constants. This is a problem we 
hope to report upon soon. Furthermore, as one considers the two-cuts solution to the quartic 
matrix model, another double-scaling limit naturally appears: that of the Painleve II equation 
describing 2d supergravity. Given that our off-critical transseries construction was very much 
attached to implementing correct double-scaling limits, this is certainly an interesting problem 
to address. Finally, we have just started uncovering what we believe is a very general method 
towards the construction of explicit nonperturbative solutions in string theory. Still within the 
matrix model realm, addressing two-matrix models and their associated minimal series seems to 
be a direction of great interest. We hope to return to many of these ideas in the near future. 
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A. The Painleve I Equation: Structural Data 



The general two-parameters transseries solution of the Painleve I equation has the form 

+00 +00 



n=0 m=0 



with 



(n\m,) 



min{n,m) 
fc=0 



(A.l) 



(A.2) 



Table 4 shows up to which order in w we have calculated ^^!'\ n. The table is for A; = 0; as we 

(n\m) ' 

will see, the results for nonzero k are directly proportional to those. Moreover, we only list the 
entries for n > m; as we shall see in a moment, the coefficients for n < m can be easily obtained 
from those with n > m. It would go too far to reproduce all the data in this appendix — the 
interested reader may request a Mathematica notebook from the authors containing all calculated 
coefficients. Below, we reproduce part of the expansions for some small values of n and m. 



\n 
m\ 





1 


2 


3 


4 


5 


6 


7 8 


9 


10 





1000 


300 


300 


300 


300 


300 


300 


25 10 


10 


10 


1 




300 


300 


300 


300 


300 


23 


24 






2 






300 


300 


300 


300 


22 


23 






3 








300 


300 


20 


21 


22 






4 










20 


19 


20 


21 






5 












20 


19 


20 






6 














20 


19 






7 
















20 







Table 4: Order in w up to which we have calculated the ^j^'i^-). 
The first few ^f'\n\ are: 

(n|0) 



(0|0) 



1 



1 

48 



w 



49 
4608 



w 



1225 
55296 



w 



12 



64V3 



75 



8192 



w 



341329 
23592960\/3 



_ 1 

(2|0) - g 



55 



w 



576V3 



1325 
36864 



w 



3363653 
53084160\/3 



The first few \-,^ are: 

(n|l) 



[0] 

(1|1) 



—w 

11 



75 
512 



300713 



w 



1572864 



w 



10 



985 



597575 



(2|i) 72 



(3|1) 16 



w 



w 



= w H w ■ 

4608^/3 15925248 



3455 



1712825 



10368^/3 7962624 



w 



(A.3) 
(A.4) 

(A.5) 



(A.6) 
(A.7) 

(A.8) 
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In ^(ri|i), one sees the first logarithms appearing. One finds that *I'(i|i) has no logarithmic terms, 
and that 

^ fil 4 5 o 75 c / . 

^ 48 2048^/3 
, fil 4 9 55 /I 1325 fi . > X 

3^ 216 4608^/3 ^ ^ 

The reader may notice that these functions are very similar to the ^(^|o) fisted above: in fact, 
using (5.16), one can easily show that the recursion relations for the coefficients of the two power 
series are the same, and so they are equal up to an overal multiplicative constant. To be precise, 
one finds that 

$W , = - ^^""^^ (A.ll) 

(n|l) ^ ("-1|0) ^ ' 

This relation was first noted in [13], and all the formulae we have tabulated so far can in fact be 
derived from the formulae in that paper. However, with our methods one can easily go beyond 
the results of [13]. At the next level, *I'(n|2)) we find for example that 

^roi 5 4 54425 o 26442605 , 

^Si) = -6 - + 82944 - - 15925248 ^ + " " " ' ^^"^^^ 

^1, = + - + • • • , (A.13) 

24^/3 20736 1769472\/3 ^ ' 

47 4 54415 6750359 . 

72^/3 124416 5971968^3 

These functions, except for the diagonal one $(2|2)) &lso have parts proportional to logtf;. They 
are 

^fH 11 3 985 5 597575 . ^ , 

(^1^^ 18^/3 3456 3981312^/3 

^[11 3 4 3455 fi 1712825 o , 

$4„. = ^w^ + -w^ + ---. A.16 

2^3 3888 995328^/3 ^ ' 

The new phenomenon at this level is that we now also have log^ w contributions. These are found 
to be 

(^1^) 3 24\/3 1024 ' ^ ^ 

^ [21 16 9 55 4 1325 fi / . N 

Again, these functions have a close relation to the functions Iii fact, with a bit of work, 

one can show from the recursion relation (5.16) that terms with a given power of log w are always 
proportional to similar terms with lower n and m, as well as lower logarithmic power, 

^[", ,=^^^^$n, (A.19) 
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where in this expression we have assumed that n > m. Applying this formula k times, one can 
further express these coefficients in terms of log-free coefficients as 



^(n|m) - k\\ ^ ) '^{n-k\m-ky (^-20) 

This immediately implies that the logarithmic sectors are, from a certain point of view, artifacts 
of the resonant transseries solution — they do not contain any new physical content. Finally, we 
remark that we have only listed above with n > m. The expansions for n < m are very 

similar^^. In fact, one finds that 

^{n\m)[k] _ ^_Y-^{g-n-m)/2^(jn\n)[k] (A. 21) 

for n ^ m. This again generalizes a similar relation found in [13]. 
B. The Quartic Matrix Model: Structural Data 

In this appendix, we present some of the explicit polynomials that determine the full nonper- 
turbative solution (6.51) to the one-cut quartic matrix model. Recall from section 6.3 that this 
solution has the form 

^(^) = EE<^™^^"|"^^(^) (B.i) 

n=0 m=0 

with 

+00 

^(n|m)^^) ^g-(n-mM(x)/9, ^ g9RM^)(x), (B.2) 

Q — ^nm 

and that the expansion coefficients R^J^^^\x) can be expressed in terms of polynomials pj'^l'^^t'^] (x) 
as 

min(n,m) 

<'"^^(")= .P.r.-.Uir.-A.^^^ E log'=(/(x))-Pi"l™)W(x). (B.3) 



rP2 (3 - 3Ar)P3 (3 - Xrf 



The following table 5 shows to which order in we have calculated the Pg'^^^\x) polynomials: 



\ n 
m\ 





1 


2 


3 


4 


5 


6 


7 


8 9 


10 





100 


30 


30 


30 


10 


10 


10 


10 


10 10 


10 


1 




12 


4 


4 


4 


3 










2 






4 


2 


2 


2 











Table 5: Values for the highest g for which we have calculated Pg"'™-*. 

Note that the numbers in this table are actually smaller than the actual number of calculated 
polynomials. For example, at n = 5 and m = 2, g starts at /3nm. = —2. Therefore, the entry 
of 2 means that we have calculated the five leading orders. At each of these orders (except for 

^^Notice that the naive observation that all Painleve I coefficients with n < m are positive is, in fact, not true 
(even though the examples we have shown could seem to point in that way). This is only noticed for the first time 
when n = 3, m = 4 and at genus 11, so it is indeed an assumption which is hard to falsify! 
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the leading one), the expression contains three polynomials multiplying different powers of the 
logarithm. Therefore, this entry of 2 corresponds to a total of 13 polynomials. 

In the table, we have only mentioned the calculated polynomials for n > m. The ones with 
n < m differ from those only by a sign, 

p{n\m)[k] _ (B.4) 

For reasons of space, in this appendix we only reproduce a very small sample of the calculated 
polynomials. A Mathematica file containing all the calculated data is available from the authors. 

Let us begin with the perturbative results — that is, n = m = 0. At this order, the data is 
most easily reproduced in terms of the polynomials Sg introduced in (6.61). For the first three 
of those, we have 

27 

^2 = y, (B.5) 

5, = ™(5 + 2X), (B.6) 
o 

177147 

56 = (1925 + 2864X + UIX'^) , (B.7) 

16 

where we substituted X = Xr. These results exactly match the results that were found in [25, 15]. 

For the one-instanton contributions, appearing at n = 1 and m = 0, we list the first four of 
the polynomials Pg^^^\ 

Po^^l°^ = 1, (B.8) 





= _ ^ (6 + 3X - 6X2 + 2x^) , (B.9) 



pirn ^ _^ 3gj^ _^ 1665^2 - 2844X3 + ISOOX"^ - 536X^ + 68X^) , (B.IO) 
128 

= (30024 - 234900X + 608958X2 _ 3803895X3 + 6142554X^- 

3 5120 ^ 

-4634370X^ + 2034360X^ - 588060X'^ + 116520X^ - 12520X^) . (B.ll) 

These expressions agree with the one-instanton results presented in [12]. 

We now turn to some of the new results. For the two-instanton case, n = 2 and m = 0, 

P^' = -I, (B.12) 

pi2\o) ^ _ _^ ^^7^ _ -^Q2X2 + 22X3) ^ (^ -^^-^ 

8 

p{2\0) ^ _81 _^ ^gg^ ^ 5577X2 _ §204X3 + 4460X'^ - 1128X5 + UeX^) , (B.14) 

= (-20088 + 238140X + 989334X2 _^ 23247945X3 - 41702958X^+ 

+29306340X5 - 10628280X^ + 2188980X^ - 276120X^ + 20360X^) . (B.15) 

The main novelty of our method is that we can also calculate contributions with generalized 
instantons, having the "wrong sign" of the instanton action. For example, for n = m = 1, we 
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have, using the notation introduced in (6.59), 



Qi'^ = 3{2-X), (B.16) 

TOO 

Q^^^ = (72 + 220X - 380^2 + 207^^ - 48X^ + AX^) , (B.17) 



^^1) ^ 59049 ^272160 + 2748816X - 5760432X2 + 4023324X3- 
128 

-724722X'^ - 548049X^ + 380368X^ - 104016X^ + 14048X* - 784X^) . (B.18) 

These results for n = m = 1 do not yet show all the features of the "generalized instanton" 
expansions. As in the Painleve I case, we find that whenever n = m, there are no "open string" 
odd g contributions. Also, in these cases, there are no logarithmic contributions yet. Finally, the 
perturbative series start at g = 0. All three of these properties disappear when we go to cases 
where n ^ m. For example, when n = 2 and jn = 1, we find 

^^(2|i)[0] = ^ (54 - 45X - 6X2 ^ 8^3) ^ (B ^g) 

p{2\im ^ ^ ^224 _ 3132X + 1197X2 + 2052X3 - 2202X^ + 940X5 - 164X6) , (B.20) 

d(2|i)[o] ^ _9_ /_52488 + 317844X + 961794X2 + 7811559X3 - 22378842X^+ 
2 512 ^ 

+23547888X5 - 13285728X'^ + 4500468X'^ - 914760X* + 89840X^) , (B.21) 

for the logarithm-free contributions, and 

p(2|i)[i] = _L, (B.22) 

p^{2|i)[i] ^ _ A (6 + 3X - 6X2 + 2jj^3^ ^ ^g_23) 
32 

p^(2|i)[i] ^ _^ ^ 1665X2 - 2844X3 + 1800X^ - 536X5 ^ g8^6^ ^ ^3 24) 

for the one-logarithm contributions. Notice that the latter polynomials are essentially the same 
as the Pg^^'^^ reported starting in (B.8). In fact, we find in general that 

This relation can be easily derived from the recursion relations that follow from the string equa- 
tion. In fact, our expression above is simply the analogue of (A. 11), and it is the first in a 
sequence of equations that are analogous to the relations (A. 20) that we have found for the 
Painleve I transseries coefficients. For the general case, one obtains 



p{n\m)[k] ^ 1 {n-m) \^ p(n-fc|m-fc)[0] 

3 k\\ 12 ) 9+k 



(B.26) 



Thus, once again, the logarithmic contributions are simply related to the logarithm-free contri- 
butions, and do not seem to constitute new physical sectors. 
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a 





1 


2 


c 


3 

Q 


3 


27 




36 





-11664 




18 


23328 


122472 


x^ 


-38 


27432 


3170988 


x^ 


11 


-73476 


9125514 


X* 




49311 


25985394 


x^ 




-14442 


24283071 


x^ 




1667 


-11842992 


x^ 






3354462 


x^ 






-544628 


X9 






40996 



9 


1 

— 1 





1 




1 

12 


1 

16 


27 

128 




1 


18 


36 


x^ 




117 


468 


x^ 




-102 


5577 


x^ 




22 


-8204 


x^ 






4460 


x^ 






-1128 


X6 






116 



Table 6: Prefactor c and coefficients of the polynomials Pi 



(3|1)[0] 



(left) and I^^I^l (right). 



g 





2 


4 


c 


9 

4 


81 
16 


6561 
64 


X» 


-72 


326592 


-687802752 


X^ 


78 


255636 


-2925199980 


x^ 


-31 


-1268946 


9776740014 


x^ 


5 


1263654 


-10514590074 


X* 




603801 


4732494984 


x^ 




154827 


148363974 


x^ 




-20062 


-1271607633 


x^ 




950 


701712243 


x^ 






-203346798 








34993318 


^10 






-3454976 








156840 



Table 7: Prefactor c and coefficients of the polynomials Q 



(2) 



Using the formula (6.78), the reader can check that, in the double-scaling limit, the data we 
have presented so far exactly reproduces the expansions for 

cf)M (]>M <i)M and c&W rR27l 

*(0|0)' *(1|0)' *(2|0)' ^(1|1)' *(2|1) *{2|1)' 

that were listed in appendix A. For completeness, we also tabulate the coefficients of all other 
polynomials that are needed to reproduce the expansion coefficients we gave in that appendix. 



C. The Double— Scaling Limit: Structural Data 

In this appendix, and analogously to the previous one, we present some of the polynomials 
P^"l"'^(t) that deter mine the free energy (6.138) of the quartic matrix model. Table 11 shows 
to which index g we have calculated these polynomials. As will be clear when comparing this 
table to the analogous table in the previous appendix, the amount of available J-g^^"^\t) data 
is much smaller than the amount of R^g^^"^\x) data. The reason for this is that the procedure 
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128 
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X" 




72 


-7776 


443232 
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69984 
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101088 


-4000752 




36 


194400 


944784 






-90 


-137700 


24782112 




-90 - 


106272 


64513584 


x^ 




29 


44280 


-22509576 




29 - 


175392 


-20419776 








24687 


-5930982 


— ^ /I 




197802 


-135263034 


x^ 






-20094 


21534309 


X5 




-73836 


182249163 


x^ 






3941 


-17087760 


X6 




9937 


-108164682 


x^ 








7682442 


X^ 






35518077 


x^ 








-2022868 


x^ 






-6475770 


x^ 








240208 


x^ 






528388 


'able 8: Prefactor c and coefficients of the polynomials pj^l^^'^ 


(left) and Pg^^l^^I^l (right 


9 




-1 





1 
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48 


3 
128 


3 

2048 
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1 

16 


1 
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9 
512 


X 
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54 


-324 


-52488 


X^ 


36 





-11664 


X 


1 


-45 


3132 


317844 


xi 


18 


23328 


122472 


X 


2 


-6 


-1197 


961794 




-38 


27432 


3170988 


x^ 


8 


-2052 


7811559 




11 - 


73476 


9125514 


X 


4 




2202 


-22378842 






49311 


-25985394 


X 


5 




-940 


23547888 






14442 


24283071 


x^ 




164 


-13285728 






1667 


-11842992 


X 


7 






4500468 








3354462 


x^ 






-914760 


X8 






-544628 


X9 






89840 


X9 






40996 



Table 9: Prefactor c and coefficients of the polynomials p^'^'^^I^l (left) and pj**'^^!^] (right). 
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-2 


-1 
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-2 


-1 





c 


1 

288 


1 

256 


9 

4096 


c 


1 

144 


1 
192 


9 
512 




1 


6 


36 


X" 


1 


18 


36 


xi 




3 


36 


xi 




117 


468 


X2 




-6 


1665 


X2 




-102 


5577 


X3 




2 


-2844 


X3 




22 


-8204 


x^ 






1800 


X4 






4460 


X5 






-536 


X5 






-1128 


X6 






68 


X6 






116 



Table 10: Prefactor c and coefficients of the polynomials Pg ' ''■ ' (left) and Pg ' ' (right). 



used to calculate Fg (t) from Rg^'{x), using the Euler-MacLaurin formula, is rather time 
consuming. We have therefore chosen to do the tests of resurgence for the quartic matrix model 

directly at the level of R^g'^'^\x), where one can construct a sufficient amount of data much more 
easily. The J-'g"'™^(t) for which the data are presented in this appendix mainly serve the purpose 
of checking that the quartic matrix model free energy gives the (2, 3) minimal model free energy 
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m 
T 

1 

2 



12 3 4 

"25 7/2 4 

5 7/2 4 

4 7/2 4 



Table 11: Values of the highest g for which we have calculated 



g 


1 


3 


5 


c 


l_ 

6v^ 


1 


1 

378v^ 




_2 




—61908 






2835 


—574056 




8 


3642 


-1614616 






-16512 


1807479 






-5472 


8602998 






1950 


17467588 






51840 


-66986172 






-36000 


39683718 






-19200 


-60738324 






16640 


220690302 








-232460928 


til 






52828048 








14853888 








35051520 








-38348800 








9805824 



9 


2 


4 


c 


1 

486 


1 

43740 




130 


396710 




516 


3402120 




2229 


12327720 




-5408 


-20516720 




1896 


12385215 




-816 


-230785920 




1696 


536735424 


f 




-513929952 






490487040 






-569834240 






320398080 


til 




-6978048 






-34264576 



Table 12: Prefactor c and coefficients of the polynomials Vg^^^^t) (left) and Vg'^^^\t) (right). 



in the double-scaling limit. 

As usual, we only present results with n > m. The results for n < m are related to those by 

J:{n\m)^^^ = (C.l) 

Results for n = m = were already listed in (6.113-6.116) in the main text. We have also listed 
two exceptional results in there, (6.125) for n = m = 1 and (6.127) for n = m = 2. For all other 
(regular) results, we give the noninteger prefactors c alongside with the integer coefficients of the 
polynomials Vg'^^'^\t), defined in (6.138), in tables 12-15. 

In the first column of each table, we list the monomial that the coefficients in that row 
multiply. The corresponding power of a is easily derived from the fact that the whole polynomial 
is homogeneous in t and a^, with the highest power a pure power of t. Thus, if nmax is the index 
of the highest coefficient in a certain column, the coefficient in the row labeled of that column 
actually multipHes t"Q/2("max-n)^ 

D. Stokes Automorphism of Two— Parameters Instanton Series 

An expression for the general ordered product of k alien derivatives, of the form Y[i=i ^-^k+i-i^ ~ 
■^-4.A ■ ■ ■ ^iiA, acting on <^'(„|o), was presented in section 4.1, namely expression (4.36). In this 
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35/4 f 


)-33/4v^ 


144.31/4^ 
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517000 








-21 


2820 


3246300 
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2697 


5408118 








40 


-9224 


-10506063 










-2208 


-15792588 










3648 


-4743720 










1600 


44745600 












-12288960 












-17130240 












6540800 


9 


1 


2 


3 


4 


c 


]_ 


1 

162 


1 

648\/3 


1 

58320 




-1 


-22 


-316 


-22520 




4 


253 


10564 


3903200 






-564 


-41715 


-18769266 






1096 


168044 


125672865 






-736 


-341936 


-421619748 








393408 


941275296 








-281920 - 


-1561721280 








93952 


1764081600 




i 








-1258640640 










530946560 










-105113600 



Table 13: Prcfactor c and coefficients of the polynomials 'Pg^^'^\t) (top) and 'Pg'^^^\t) (bottom). 



appendix we shall outline an inductive proof of this result. First recall what this expression was, 

k k k 



i=l 



■4+l-,^^(n|0) 



E E n 

m=0 5ser{k,k-m+l) s=l 



(s + 1 - 6s) Sif + 



n 



i=l 



e(s + i-ds 



(D.l) 



i+m-J]jLi4|m)' 



Further recall that in section 4.1 we have explicitly shown that for the case of k = 2 (and 
analogously for the case of A; = 1) this closed form expression correctly reproduced the result we 
had earlier computed in (4.32). 

Assuming that the above result (D.l) holds true for a particular value of A; > 2, let us apply 
one more alien derivative A_£i^^^j^, with ik+i > 0, to this expression. Notice that this alien 
derivative, A_£, , will only act on v-* « 1^-,, and this action was already computed in 

(4.33). We thus find 



k+l 



fc+1 



4+l-,:A$(n|0) 



i=l 



m=0 



i=l 



q=0 



i=l 
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c 



3/2 5/2 7/2 

V2 1 1 



-14 
-21 
-96 
104 



184 



2665 
-3972 
-144 
-2848 

3200 



-638120 
-4532580 

-12820266 
13158375 
10689480 
70972776 

-114864000 
24333120 
11455488 
2252288 



g 


2 


3 


4 


c 


2 

81 


2 

729^/3 


1 

2916 




4 


341 


5032 




-49 


-6408 


-271696 




144 


25197 


1575868 




-310 


-113287 


-10610537 




184 


173664 


36078160 






-187692 


-86311034 






95792 


150067240 


f 






-173970320 








130130752 








-59481984 








12789248 



Table 14: Prefactor c and coefBcients of the polynomials Vg'^^^\t) (left) and Vg^^^\t) (right). The case 
V^'^^^\t) is exceptional, in the sense that it factorizes: the polynomial displayed here should be multiplied 
by (t-a^) to obtain T'f ^^(t)- 



Table 15: 
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27-33/4 


729-3V4v^ 


1944-33/4v^ 




-2 


896 


-39752 




6 
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24537 


-1595714 




29 


-30592 


1325412 






38919 


-3580714 






-79788 


17588671 






42836 


-19468680 


f 






8762744 








-5577248 








2872832 



9 


2 


3 


4 


c 


4 
243 


2 

2187v^ 


1 

39366 




-2 


-488 


-75448 




6 


9686 


2192188 




-42 


-54240 


-15529806 




29 


242728 


105787830 






-526987 


-368392458 






765960 


949167207 






-650210 


-1732001196 


f 




214280 


2088163092 








-1656852624 








806826880 








-179279104 



Prefactor c and coefficients of the polynomials 'Pg^^'^\t) (left) and 'Pg^'^'(t) (right). The case 



(t) is exceptional, in the sense that it factorizes: the polynomial displayed here should be multiplied 



by {At — a^) to obtain V2 (t). It is curious to see that the remaining factor is proportional to Vj^^l^ (t) 



<5ser(fc,fc-m-|-2) s=l 



+ ln-Y,^^ + s + l-Ss]s[ 



i=l 



(s + 1 - 5s) + (D.2) 



-q\m—qj 



In order to obtain the expression above, we have changed the variable in the first sum of (D.l) 
from X^^^o ~^ S?v=i' ^fter which one realizes that one may always add the term m' = as it 
is zero. The next steps include the change of variables X^^q ~ Sq^=m-g=o noticing that 
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one can further change the order of the sums as 



fc+l m k+1 k+1 

EE = EE- 

m=0 g'=0 q'=Qm=q' 



(D.3) 



In this process we thus obtain 

k 

^-^k+iAW^-lk+i-,A^{n\0) = 



i=l 



X E n 

6s&T{k,k-m.+2) s=l 



k+1 k+l / / k+1 \ \ 

E E "+o'-E'.Uitr-''' + «'-S5-f X 

q'=Om=q' \ \ i=l / / 



i=l 



_ ^k+2-q' 



The final step is to change variables yet one more time, as = ^,„'=fc+2-m=i' 

a new variable, jk+i = k + 2 — q' . Then 



(D.4) 



and introduce 



fc+i 



k+1 



^-4+1^11 ^-^fe+i-^^(«|o) = ^K+i,k+2-q' ( ( n + A: + 2-7fc+i 

i=l q'=0 



m'=l 



X ') + + 2 - jk+i) s^i;^"^ E n 



-44 



i=l 



Yli + s + l-5s 



i=l 



/ (5,er(fc,m') s=i 



m—q' 



(D.5) 



Finally recalling that 5s G r(A;, m') means that we are summing over all (5s : < < • • • < 5^ = 
m', and that now m' = 1, • • • , 7a;+i = A; + 2 — q'</c + 2, one can naturally rewrite the above 
expression as a sum over Young diagrams, of length /c + 1, obtaining 



h k+1 k+1 

^-4+iAn^-4+i-Aio) = E E n 

i=l q'=0 6ser{k+l,k+2-q') s=l 



(, + 1 _ s,) ^if + (D.6) 



+ ln-Y,^i + s + l-6s] S_ 



(is+dSs) 



i=l 



+q' m—q' ) ' 



This is the expected result for the ordered product of A; + 1 alien derivatives, acting on the 
instanton series ^(„|o)) as shown in (D.l). It thereby concludes our proof. 
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